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Abstract 

The article reviews recent developments in the theory of fluctuations and corre- 
lations of energy levels and eigenfunction amplitudes in diffusive mesoscopic samples. 
Various spatial geometries are considered, with emphasis on low-dimensional (quasi-lD 
and 2D) systems. Calculations are based on the supermatrix cr-model approach. The 
method reproduces, in so-called zero-mode approximation, the universal random ma- 
trix theory (RMT) results for the energy-level and eigenfunction fluctuations. Going 
beyond this approximation allows us to study system-specific deviations from univer- 
sality, which are determined by the diffusive classical dynamics in the system. These 
deviations are especially strong in the far "tails" of the distribution function of the 
eigenfunction amplitudes (as well as of some related quantities, such as local density of 
states, relaxation time, etc.). These asymptotic "tails" are governed by anomalously 
localized states which are formed in rare realizations of the random potential. The de- 
viations of the level and eigenfunction statistics from their RMT form strengthen with 
increasing disorder and become especially pronounced at the Anderson metal-insulator 
transition. In this regime, the wave functions are multifractal, while the level statis- 
tics acquires a scale-independent form with distinct critical features. Fluctuations of 
the conductance and of the local intensity of a classical wave radiated by a point-like 
source in the quasi-lD geometry are also studied within the cr-model approach. For a 
ballistic system with rough surface an appropriately modified ("ballistic") cr-model is 
used. Finally, the interplay of the fluctuations and the electron-electron interaction in 
small samples is discussed, with application to the Coulomb blockade spectra. 

PACS numbers: 05.45.Mt, 71.23.An, 71.30.+h, 72.15.Rn, 73.23.-b, 73.23.Ad, 73.23.Hk 
Keywords: Level correlations, wave function statistics, disordered mesoscopic systems, su- 
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1 Introduction 

Statistical properties of energy levels and eigenfunctions of complex quantum systems have 
been attracting a lot of interest of physicists since the work of Wigner [|lj , who formulated 
a statistical point of view on nuclear spectra. In order to describe excitation spectra of 
complex nuclei, Wigner proposed to replace a complicated and unknown Hamiltonian by a 
large N x N random matrix. This was a beginning of the random matrix theory (RMT), 
which was further developed by Dyson and Mehta in the early 1960's 0, |3|. This theory 
predicts a universal form of the spectral correlation functions determined solely by some 
global symmetries of the system (time- reversal invariance and value of the spin). 

Later it was realized that the random matrix theory is not restricted to strongly inter- 
acting many-body systems, but has a much broader range of applicability. In particular, 
Bohigas, Giannoni, and Schmit f§] put forward a conjecture (strongly supported by accu- 
mulated numerical evidence) that the RMT describes adequately statistical properties of 
spectra of quantum systems whose classical analogs are chaotic. 

Another class of systems to which the RMT applies and which is of special interest to us 
here is that of disordered systems. More specifically, we mean a quantum particle (an elec- 
tron) moving in a random potential created by some kind of impurities. It was conjectured 
by Gor'kov and Eliashberg || that statistical properties of the energy levels in such a disor- 
dered granule can be described by the random matrix theory. This statement had remained 
in the status of conjecture until 1982, when it was proven by Efetov |J. This became possible 
due to development by Efetov of a very powerful tool of treatment of the disordered systems 
under consideration — the supersymmetry method (see the review || and the recent book 
0). This method allows one to map the problem of the particle in a random potential onto 
a certain deterministic field-theoretical model (supermatrix a-model), which generates the 
disorder-averaged correlation functions of the original problem. As Efetov showed, under 
certain conditions one can neglect spatial variation of the cr-model supermatrix field (so- 
called zero- mode approximation), which allows one to calculate the correlation functions. 
The corresponding results for the two-level correlation function reproduced precisely the 
RMT results of Dyson. 

The supersymmetry method can be also applied to the problems of the RMT-type. In 
this connection, we refer the reader to the paper [EL where the technical aspects of the 
method are discussed in detail. 

More recently, focus of the research interest was shifted from the proof of the applicability 
of RMT to the study of system-specific deviations from the universal (RMT) behavior. For 
the problem of level correlations in a disordered system, this question was addressed for 
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the first time by Altshuler and Shklovskii |9j in the framework of the diffuson-cooperon 
diagrammatic perturbation theory. They showed that the diffusive motion of the particle 
leads to a high-frequency behavior of the level correlation function completely different from 
its RMT form. Their perturbative treatment was however restricted to frequencies much 
larger than the level spacing and was not able to reproduce the oscillatory contribution to 
the level correlation function. Inclusion of non-zero spatial modes (which means going beyond 
universality) within the a-model treatment of the level correlation function was performed 
in Ref. flOfl . The method developed in [1C] was later used for calculation of deviations from 
the RMT of various statistical characteristics of a disordered system. For the case of level 
statistics, the calculation of 
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valid for not too large frequencies (below the Thouless 
energy equal to the inverse time of diffusion through the system) was complemented by 
Andreev and Altshuler 
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whose saddle-point treatment was, in contrast, applicable for 
large frequencies. Level statistics in diffusive disordered samples is discussed in detail in 
Section |2] of the present article. 

Not only the energy levels statistics but also the statistical properties of wave functions 
are of considerable interest. In the case of nuclear spectra, they determine fluctuations of 
widths and heights of the resonances | 12| . In the case of disordered (or chaotic) electronic 



systems, eigenfunction fluctuations govern, in particular, statistics of the tunnel conductance 
in the Coulomb blockade regime [13|. Note also that the eigenfunction amplitude can be 

(though in this case one 



14, 15, 1 



directly measured in microwave cavity experiments 
considers the intensity of a classical wave rather than of a quantum particle, all the results 
are equally applicable; see also Section [5]). Within the random matrix theory, the distribution 
of eigenvector amplitudes is simply Gaussian, leading to \ 2 distribution of the "intensities" 
2 (Porter-Thomas distribution) 



0- 



A theoretical study of the eigenfunction statistics in a disordered system is again possible 
with use of the supersymmetry method. The corresponding formalism, which was developed 
in Refs. [|17], [18], [19], [2(J (see Section |37T|) , allows one to express various distribution functions 
characterizing the eigenfunction statistics through the cx-model correlators. As in the case 
of the level correlation function, the zero-mode approximation to the a-model reproduces 
the RMT results, in particular the Porter-Thomas distribution of eigenfunction amplitudes. 
However, one can go beyond this approximation. In particular, in the case of a quasi-one- 
dimensional geometry, considered in Section |3.2| , this cr-model has been solved exactly using 
the transfer-matrix method, yielding exact analytical results for the eigenfunction statistics 
for arbitrary length of the system, from weak to strong localization regime PTl IT? 



The case of a quasi- ID geometry is of great interest not only from the point of view of 
condensed matter theory (as a model of a disordered wire) but also for quantum chaos. 

In Section |3.3| we consider the case of arbitrary spatial dimensionality of the system. 
Since for d > 1 an exact solution of the problem cannot be found, one has to use some 

m 



approximate methods. In Refs. |23|, |25| the scheme of [|I(| was generalized to the case of 
the eigenfunction statistics. This allowed us to calculate the distribution of eigenfunction 
intensities and its deviation from the universal (Porter-Thomas) form. Fluctuations of the 
inverse participation ratio and long-range correlations of the eigenfunction amplitudes, which 
are determined by the diffusive dynamics in the corresponding classical system [^, E7[ ], 
are considered in Sec. [3.3.3 . 



Section [| is devoted to the asymptotic "tails" of the distribution functions of various 
fluctuating quantities (local amplitude of an eigenfunction, relaxation time, local density of 
states) characterizing a disordered system. It turns out that the asymptotics of all these 
distribution functions are determined by rare realizations of disorder leading to formation 
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of anomalously localized eigenstates. These states show some kind of localization while all 
"normal" states are ergodic; in the quasi-one- dimensional case they have an effective local- 
ization length much shorter than the "normal" one. Existence of such states was conjectured 
by Altshuler, Kravtsov, and Lerner [2S] who studied distributions of various quantities in 
2 + e dimensions via the renormalization group approach. More recently, Muzykantskii and 
Khmelnitskii [p9] suggested a new approach to the problem. Within this method, the asymp- 
totic "tails" of the distribution functions are obtained by finding a non-trivial saddle-point 
configuration of the supersymmetric a-model. Further development and generalization of 
the method allowed one to calculate the asymptotic behavior of the distribution functions 
of relaxation times |29|, [30], [H]], eigenfunction intensities [[32], |33|, local density of states f34|j . 
inverse participation ratio |35|, [36f| , level curvatures |$7], |38| etc. The saddle-point solution 
describes directly the spatial shape of the corresponding anomalously localized state |29|, . 

Section |5| deals with statistical properties of the energy levels and wave functions at 
the Anderson metal- insulator transition point. As is well known, in d > 2 dimensions a 
disordered system undergoes, with increasing strength of disorder, a transition from the 
phase of extended states to that of localized states (see e.g. 



35 for review). This transition 



changes drastically the statistics of energy levels and eigenfunctions. While in the delocalized 
phase the levels repel each other strongly and their statistics is described by RMT (up to 
the deviations discussed above and in Section 0), in the localized regime the level repulsion 
disappears (since states nearby in energy are located far from each other in real space). As 
a result, the levels form an ideal ID gas (on the energy axis) obeying the Poisson statistics. 
In particular, the variance of the number N of levels in an interval AE increases linearly, 
var( N) = (N) , in contrast to the slow logarithmic increase in the RMT case. What happens 
to the level statistics at the transition point? This question was addressed for the first 
time by Altshuler, Zharekeshev, Kotochigova, and Shklovskii in where a Poisson-like 
increase, var(iV) = x(N), was found numerically with a spectral compressibility x — 0.3. 
More recently, Shklovskii et al [41] put forward the conjecture that the nearest level spacing 
distribution P(s) has a universal form at the critical point, combining the RMT-like level 
repulsion at small s with the Poisson-like behavior at large s. However, these results were 
questioned by Kravtsov et al [4"2" who developed an analytical approach to the problem and 
found, in particular, a sublinear increase of var(iV). This controversy was resolved in [[13], [44 
where the consideration of JEJ was critically reconsidered and the level number variance was 
shown to have generally a linear behavior at the transition point. By now, this result has 



been confirmed by numerical simulations done by several groups 



|S). Recently, 



a connection between this behavior and multifractal properties of eigenfunctions has been 
conjectured . 

Multifractality is a formal way to characterize strong fluctuations of the wave function 
amplitude at the mobility edge. It follows from the renormalization group calculation of 
Wegner |5(| (though the term "multifractality" was not used there). Later the multifractality 
of the critical wave functions was discussed in [51] and confirmed by numerical simulations 
of the disordered tight-binding model in p2l R3|, |53, [551 156| . It implies, in very rough terms, 



that the eigenfunction is effectively located in a vanishingly small portion of the system 
volume. A natural question then arises: why do such extremely sparse eigenfunctions show 
the same strong level repulsion as the ergodic states in the RMT? This problem is addressed 
in Sec. 5T. It is shown there that the wavefunctions of nearby- in-energy states exhibit very 
strong correlations (they have essentially the same multifractal structure), which preserves 
the level repulsion despite the sparsity of the wave functions. 

In Sec. [5]2] we consider a "power-law random banded matrix ensemble" (PRBM) which 
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describes a kind of one-dimensional system with a long-range hopping whose amplitude 
decreases as r~ a with distance [j57|. Such a random matrix ensemble arises in various contexts 
in the theory of quantum chaos [58, 5S] and disordered systems [B(], 61, 62]. The problem 



can again be mapped onto a supersymmetric a-model. It is further shown that at a = 1 
the system is at a critical point of the localization-delocalization transition. More precisely, 
there exists a whole family of such critical points labeled by the coupling constant of the a- 
model (which can be in turn related to the parameters of the microscopic PRBM ensemble). 
Statistics of levels and eigenfunctions in this model are studied. At the critical point they 
show the critical features discussed above (such as the multifractality of eigenfunctions and 
a finite spectral compressibility < x < 1). 

The energy level and eigenfunction statistics characterize the spectrum of an isolated 
sample. For an open system (coupled to external conducting leads), different quantities 
become physically relevant. In particular, we have already mentioned the distributions of 
the local density of states and of the relaxation times discussed in Section |] in connection 
with anomalously localized states. In Section |6] we consider one of the most famous issues 
in the physics of mesoscopic systems, namely that of conductance fluctuations. We focus 
on the case of the quasi-one-dimensional geometry. The underlying microscopic model de- 
scribing a disordered wire coupled to freely propagating modes in the leads was proposed 
by Iida, Weidenmuller, and Zuk f63fl . It can be mapped onto a ID cx-model with boundary 
terms representing coupling to the leads. The conductance is given in this approach by 
the multichannel Landauer-Biittiker formula. The average conductance (g) of this system 
for arbitrary value of the ratio of its length L to the localization length £ was calculated 
by Zirnbauer pi ], who developed for this purpose the Fourier analysis on supersymmetric 
manifolds. The variance of the conductance was calculated in |)5| (in the case of a system 
with strong spin-orbit interaction there was a subtle error in the papers corrected by 

Brouwer and Frahm |66[ ). The analytical results which describe the whole range of L/£ from 
the weak localization (L <C £) to the strong localization (L ^> £) regime were confirmed by 
numerical simulations |67|, |68| . 

As has been already mentioned, the cx-model formalism is not restricted to quantum- 
mechanical particles, but is equally applicable to classical waves. Section [7] deals with a 
problem of intensity distribution in the optics of disordered media. In an optical experiment, 
a source and a detector of the radiation can be placed in the bulk of disordered media. 
The distribution of the detected intensity is then described in the leading approximation 
by the Rayleigh law ||6U| which follows from the assumption of a random superposition of 
independent traveling waves. This result can be also reproduced within the diagrammatic 
technique ]7D[ . Deviations from the Rayleigh distribution governed by the diffusive dynamics 
were studied in [71j for the quasi- ID geometry. When the source and the detector are 
moved toward the opposite edges of the sample, the intensity distribution transforms into 
the distribution of transmission coefficients ]7^, [73], [74| . 

Recently, it has been suggested by Muzykantskii and Khmelnitskii [75] that the supersym- 
metric cx-model approach developed previously for the diffusive systems is also applicable in 
the case of ballistic systems. Muzykantskii and Khmelnitskii derived the "ballistic cr-model" 
where the diffusion operator was replaced by the Liouville operator governing the ballistic 
dynamics of the corresponding classical system. This idea was further developed by An- 
dreev, Agam, Simons and Altshuler [76, 77] who derived the same action via the energy 
averaging for a chaotic ballistic system with no disorder. (There are some indications that 
one has to include in consideration certain amount of disorder to justify the derivation of 

m 



[77].) Andreev et al replaced, in this case, the Liouville operator by its regularization 



6 



known as Perron- Frobenius operator. However, this approach has failed to provide explicit 
analytical results for any particular chaotic billiard so far. This is because the eigenvalues 
of the Perron- Frobenius operator are usually not known, while its eigenfunctions are highly 
singular. To overcome these difficulties and to make a further analytical progress possible, 
a ballistic model with surface disorder was considered in [|7|, [79j . The corresponding results 
are reviewed in Sec. ||. It is assumed that roughness of the sample surface leads to the diffu- 
sive surface scattering, modelling a ballistic system with strongly chaotic classical dynamics. 
Considering the simplest (circular) shape of the system allows one to find the spectrum of 
the corresponding Liouville operator and to study statistical properties of energy levels and 
eigenfunctions. The results for the level statistics show important differences as compared 
to the case of a diffusive system and are in agreement with arguments of Berry [^, |3T] 
concerning the spectral statistics in a generic chaotic billiard. 

In Section |j| we discuss a combined effect of the level and eigenfunction fluctuations 
and the electron-electron interaction on thermodynamic properties of quantum dots. Sec- 
tion [9J] is devoted to statistics of the so-called addition spectrum of a quantum dot in the 
Coulomb blockade regime. The addition spectrum, which is determined by the positions of 
the Coulomb blockade conductance peaks with varying gate voltage, corresponds to a suc- 
cessive addition of electrons to the dot coupled very weakly to the outside world |82| . The 
two important energy scales characterizing such a dot are the charging energy e 2 /C and the 
electron level spacing A (the former being much larger than the latter for a dot with large 
number of electrons). Statistical properties of the addition spectrum were experimentally 
studied for the first time by Sivan et al [|83f| . It was conjectured in Ref. [83[ that fluctuations 
in the addition spectrum are of the order of e 2 /C and are thus of classical origin. However, it 



was found in Ref. p4| , jSq] that this is not the case and that the magnitude of fluctuations is 
set by the level spacing A, as in the non-interacting case. The interaction modifies, however, 
the shape of the distribution function. In particular, it is responsible for breaking the spin 
degeneracy of the quantum dot spectrum. These results have been confirmed recently by 
thorough experimental studies |86|, |87|] . 

The research activity in the field of disordered mesoscopic systems, random matrix theory, 
and quantum chaos has been growing enormously during the recent years, so that a review 
article clearly cannot give an account of the progress in the whole field. Many of the topics 
which are not covered here have been extensively discussed in the recent reviews by Beenakker 



|38fl and by Guhr, Miiller-Groeling, and Weidenmiiller 



2 Energy level statistics: Random matrix theory and 
beyond 

2.1 Supersymmetric cr-model formalism 

The problem of energy level correlations has been attracting a lot of research interest since 
the work of Wigner J]J . The random matrix theory (RMT) developed by Wigner, Dyson, 
and Mehta |2|, |3| was found to describe well the level statistics of various classes of complex 
systems. In particular, in 1965 Gor'kov and Eliashberg || put forward a conjecture that the 
RMT is applicable to the problem of energy level correlations of a quantum particle moving 
in a random potential. To prove this hypothesis, Efetov developed the supersymmetry 
approach to the problem || [7J . The quantity of primary interest is the two-level correlation 
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function?] 

R ( s ) = TZ»M E ~ w / 2 >(£ + w/2)> (2.1) 

where v(E) = V~ lr Ti5(E — H) is the density of states at the energy E, V is the system 
volume, H is the Hamiltonian, A = 1/ (v)V is the mean level spacing, s = uj/A, and (. . .) 
denote averaging over realizations of the random potential. As was shown by Efetov 0, 
the correlator ( |2.1| ) can be expressed in terms of a Green function of certain supermatrix 
cr-model. Depending on whether the time reversal and spin rotation symmetries are broken 
or not, one of three different cr-models is relevant, with unitary, orthogonal or symplectic 
symmetry group. We will consider first the technically simplest case of the unitary symmetry 
(corresponding to the broken time reversal invariance); the results for two other cases will 
be presented at the end. 

We give only a brief sketch of the derivation of the expression for R(s) in terms of the cr- 
model. One begins with representing the density of states in terms of the Green's functions, 

P{E) = 2^vj ddr ^ (r ' r) " ^ (r ' r) ] ' (2 ' 2) 

where 

Gl A (v 1 ,v 2 ) = {v 1 \(E-H±i V )- l \v 2 ) , r/^+0. (2.3) 
The Hamiltonian H consists of the free part Hq and the disorder potential U(r): 

H = H + U(r) ; H = ^p 2 ; (2.4) 

the latter being defined by the correlator 

(U(r)U(r')} = JL*(r - r') (2.5) 

A non-trivial part of the calculation is the averaging of the GrGa terms entering the corre- 
lation function {v(E + u/2)v{E — uj/2)). The following steps are: 

i) to write the product of the Green's functions in terms of the integral over a supervector 
field $ = (S , i,Xi,5 , 2 ,X2): 

G^ /2 (r 1 ,r 1 )G^ /2 (r 2 ,r 2 ) = J D$ D& S 1 (t 1 )S* 1 (t 1 )S 2 (t 2 )S* 2 (t 2 ) 

xexpjz J dv'^{v')K ll2 [E + (cu/2 + ^)A - i/]A 1/2 $(r') | , (2.6) 

where A = diag{l, 1, —1, —1}; 

ii) to average over the disorder; 

iii) to introduce a 4 x 4 supermatrix variable 1Z^ v {r) conjugate to the tensor product 
<Mr)$t(r); 

iv) to integrate out the $ fields ; 



x The two-level correlation function is conventionally denoted j3| || as i?2(s)- Since we will not consider 
higher-order correlation functions, we will omit the subscript "2" . 
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v) to use the saddle-point approximation which leads to the following equation for 1Z: 

1 



ft(r) = 7; — #( r > r ) ; 

Z71VT 

^(ri,r 2 ) = (r 1 |(E-Fo-^)- 1 |r 2 ) 
The relevant set of the solutions (the saddle-point manifold) has the form: 



n 



a- I Q 

2r 



(2.7) 
(2.8) 

(2.9) 



where / is the unity matrix, a is certain constant, and the 4x4 supermatrix Q = T~ l AT 
satisfies the condition Q 2 = 1, with T belonging to the coset space U(l, 1|2)/Z7(1|1) x Z7(l|l). 
The expression for the two-level correlation function R(s) then reads 

2 



R(s)- 
x exp 



1 

w 

T 



Re / DQ(r) 



d d rStiQAk 



J d d rStr [-DiyQf - 2iukQ\ 



(2.10) 



Here k = diagjl, —1, 1, —1}, Str denotes the supertrace, and D is the classical diffusion 
constant. We do not give here a detailed description of the model and mathematical entities 
involved, which can be found e.g. in Refs.|| [F], §, |9(J, and restrict ourselves to a qualitative 
discussion of the structure of the matrix Q. The size 4 of the matrix is due to (i) two types of 
the Green functions (advanced and retarded) entering the correlation function Q2.1| ), and (ii) 
necessity to introduce bosonic and fermionic degrees of freedom to represent these Green's 
function in terms of a functional integral. The matrix Q consists thus of four 2x2 blocks 
according to its advanced-retarded structure, each of them being a supermatrix in the boson- 
fermion space. 

To proceed further, Efetov || neglected spatial variation of the supermatrix field Q(r) and 
approximated the functional integral in Eq.( |2.10| ) by an integral over a single supermatrix 
Q (so-called zero-mode approximation). The resulting integral can be calculated yielding 
precisely the Wigner-Dyson distributionF]: 



R 



SUl 7TS , 



WD l 



(tts) 2 

the superscript U standing for the unitary ensemble, 
orthogonal (O) and the symplectic (Sp) ensemble are 



(2.11) 

The corresponding results for the 



R 



( 

WD I 



sin ( 7xs 



R 



Sp . 
WD I 



[ITS 
,2 



11 ■ f ^ 

-sign(s) 



1 - 



sin (27rs) 



smy 



(2vrs) 
dy . 



+ Si(27rs) 



cos 2tts 
2tts 



COS 71S 
71 S 

sin 2irs 
(2ns) 2 _ 



Sin 7TS 



7TS 



(2.12) 
(2.13) 



The aim of Sec. |2.2| will be to study the deviations of the level correlation function from the 
universal RMT results ( pUD , (gig) , (pl3|) . 



Strictly speaking, the level correlation functions (2.11 )— ( 2. 13 ) contain an addititonal term 5(s) corre- 



sponding to the "self-correlation" of an energy level. Furthermore, in the symplectic case all the levels are 



double degenerate (Kramers degeneracy). This degeneracy is disregarded in (2.13) which thus represents the 



correlation function of distinct levels only, normalized to the corresponding level spacing. 
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2.2 Deviations from universality 

The procedure we are using in order to calculate deviations from the universality is as follows 
Tcfl . We first decompose Q into the constant part Qo and the contribution Q of higher modes 



with non-zero momenta. Then we use the renormalization group ideas and integrate out 
all fast modes. This can be done perturbatively provided the dimensionless (measured in 
units of e 2 /h) conductance g = 2nE c / A = 2nuDL d ~ 2 > 1 (here E c = D/L 2 is the Thouless 
energy). As a result, we get an integral over the matrix Q only, which has to be calculated 
non-perturbatively. 

We begin with presenting the correlator R(s) in the form 

R(s) = tA?^T / DQexp{-S[Q]}\ u=0 ; 



(2ni) 2 du 2 

S[Q] = ~J Str(VQ) 2 + S J Str AQ + u J StiQAk (2.14) 

where 1/t — ttuD/4, s = irs/2iV, u = iru/2iV. Now we decompose Q in the following way 

Q(r) = T-'Q^Ta (2.15) 

where T is a spatially uniform matrix and Q describes all modes with non-zero momenta. 
When io E c , the matrix Q fluctuates only weakly near the origin A of the coset space. 
In the leading order, Q = A, thus reducing (|2.14|) to a zero-dimensional cr-model, which 
leads to the Wigner-Dyson distribution (2.11). To find the corrections, we should expand 
the matrix Q around the origin A: 

Q = A(l + W72)(l- W72)- 1 

/ W 2 W 3 \ 

= A[l + W + — + — + ... , (2.16) 



2 4 

where W is a supermatrix with the following block structure: 



*21 

Substituting this expansion into Eq. fl2.14p , we get 

S = So + S 1 + 0(W 3 ) 
'1 



W = ( ? ^ ) (2.17) 



S = I Str 



t iywy + SQoA + uQoAk 



Si = ^J Sti[SU AW 2 + uU 0k AW 2 ] (2.18) 

where Q = T _1 AT , U = T AT _1 , Uok = T AkT . Let us define 5 e ff[(5o] as a result of 
elimination of the fast modes: 

e -S cfi [Qo] = e -So[Q ]^ e -S 1 [Qo,W]+lnJ[W]^ w ^ (2.1$) 

where (. . .)w denote the integration over W and J[W] is the Jacobian of the transformation 
( p,15| ), ( p. 16] ) from the variable Q to {Qo, W} (the Jacobian does not contribute to the leading 



10 



order correction calculated here, but is important for higher-order calculations [F2^, pTf ) . 
Expanding up to the order W 4 , we get 

S cS = S + (St) - \(Sl) + ^(S x f + . . . (2.20) 

The integral over the fast modes can be calculated now using the Wick theorem and the 
contraction rules || ^8) : 

(StrW(r)PW(r')R) = II(r, r')(StrPStrP - StrPAStrPA) ; 
(Str[W(r)P]StT[W(r')R]) = n(r, r')Str(PP - PAPA), (2.21) 

where P and R are arbitrary supermatrices. The diffusion propagator II is the solution of 
the diffusion equation 

- PV 2 n( ri , r 2 ) = (ttz/)- 1 ^!-! - r 2 ) - V- 1 ] (2.22) 

with the Neumann boundary condition (normal derivative equal to zero at the sample bound- 
ary) and can be presented in the form 

U ^ = ^ E T Mr)Mr ' ] (2 - 23) 

where M (r) are the eigenfunctions of the diffusion operator —DV 2 corresponding to the 
eigenvalues e M (equal to Dq 2 for a rectangular geometry). As a result, we find: 

(Si) = ; 

(Si) = I J drdrTI 2 (r,r')(sStrQ A + uStiQ Ak) 2 (2.24) 
Substitution of Eq. flggg ) into Eq. ( p^0|) yields 

Se«[Qo) = JsStrQ A + ^uStiQ Ak + ^-^(sStrQoA + uStiQ Ak) 2 ; 
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ad 



2i 2i " 4g 

2 



AV 2 



drdrTI^r, r') = V ^ . (2.25) 

ni, ...,n d = 
+ . . . + > 

The value of the coefficient a<i depends on spatial dimensionality d and on the sample ge- 
ometry; in the last line of Eq.( |2.25| ) we assumed a cubic sample with hard- wall boundary 
conditions. Then for d = 1,2,3 we have a x = 1/90 ~ 0.0111, a 2 — 0.0266, and a 3 ~ 0.0527 
respectively. In the case of a cubic sample with periodic boundary conditions we get instead 

1 ~ 1 . . 

ad ~W7 ^ (n 2 + ... + n 2 y ' ^ 

n 1 ,...,n d = -oo 

n\ + . . . + n\ > 

so that ai = 1/720 ~ 0.00139, a 2 ~ 0.00387, and a 3 ~ 0.0106. Note that for d < 4 the sum 
in Eqs. ( |2.25|) , (|2.26|) converges, so that no ultraviolet cut-off is needed. 
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Using now Eq. (|2.14 ) and calculating the remaining integral over the supermatrix Qq, we 
finally get the following expression for the correlator to the 1/g 2 order: 



R(s) 



sin ( 7is ) 



7TS 



+ 



4a d 



sm ( 7fs 



(2.27) 



The last term in Eq. (|2.27 ) just represents the correction of order 1/g 2 to the Wigner-Dyson 
distribution. The formula (|2.27 ) is valid for s <^ g. Let us note that the smooth (non- 
oscillating) part of this correction in the region 1 <C s <C g can be found by using purely 
perturbative approach of Altshuler and Shklovskii |], [4pj . For s ^> 1 the leading perturbative 
contribution to R(s) is given by a two-diffuson diagram: 



Ras(s) - 1 



A* 

2tt 2 



Re 



E 



1 



qi = TTTli/L 

^ = 0,1,2,. 
1 ReV" 1 

) 7r 2 Z— / f — o'e _|_ ( <jr 



(Dcf - iu) 2 



2vr 2 ^ [-is + (n/2)gn^ 



(2.28) 



At s <C g this expression is dominated by the q 
of order 1/g 2 : 

Ras(s) = 1 - 



term, with other terms giving a correction 



2tt 2 s 2 



+ 



2a n 



(2.29) 



where was defined in Eq. fl2.25D . This formula is obtained in the region 1 -C s g 
and is perturbative in both 1/s and 1/g. It does not contain oscillations (which cannot be 
found perturbatively) and gives no information about actual small-s behavior of R(s). The 
result (|2.27|) is much stronger: it represents the exact (non-perturbative in 1/s) form of the 
correction in the whole region s ^ g. 

The important feature of Eq.( |2.27D is that it relates corrections to the smooth and os- 
cillatory parts of the level correlation function (represented by the contributions to the last 
term proportional to unity and to cos27rs respectively). While appearing naturally in the 
framework of the supersymmetric cr-model, this fact is highly non-trivial from the point of 
view of semiclassical theory [SO, ST], which represents the level structure factor K(t) (Fourier 
transform of R{s)) in terms of a sum over periodic orbits. The smooth part of R{s) cor- 
responds then to the small-r behavior of K(r), which is related to the properties of short 
periodic orbits. On the other hand, the oscillatory part of R{s) is related to the behavior of 
K{t) in the vicinity of the Heisenberg time r = 2n (t = 2n/A in dimensionful units), and 
thus to the properties of long periodic orbits. 

The calculation presented above can be straightforwardly generalized to the other sym- 
metry classes. The result can be presented in a form valid for all the three casesPl: 



2cid d 2 2 



(2.30) 



where (3 = 1(2,4) for the orthogonal (unitary, symplectic) symmetry; R^d denotes the 
corresponding Wigner-Dyson distribution ( p.ll|) -( p.l3|) . 

3 For all the ensembles, we denote by g the conductance per one spin projection: g = 2m/ DL d ~ 2 , without 
multiplication by factor 2 due to the spin. 
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For s — > the Wigner-Dyson distribution has the following behavior: 

R$ D ~ c^; s -> 

7T 2 7T 2 (2vr) 4 , . 

Cl = - ; c 2 = y ; c 4 = L^- (2.31) 

As is clear from Eq. fl2.30|) , the found correction does not change the exponent j3, but renor- 
malizes the prefactor eg: 

^»M=(l + 2(g + 2 f +1) g)c^ ; S ^0 (2.32) 

The correction to cp is positive, which means that the level repulsion becomes weaker. This 
is related to a tendency of eigenfunctions to localization with decreasing g. 

What is the behavior of the level correlation function in its high-frequency tail s ^> 
gl The non-oscillatory part of R(s) in this region follows from the Altshuler-Shklovskii 
perturbative formula ( j2.28| ). For s ^> g summation can be replaced by integration, yielding 

i? AS ( S )-loc<T d/ V/ 2 - 2 (2.33) 

(note that in 2D the coefficient of the term ( |2.33| ) vanishes, and the result for R^s is smaller by 
an additional factor 1/g, see [44[). What is the fate of the oscillations in R(s) in this regime? 



The answer to this question was given by Andreev and Altshuler |TI] who calculated R(s) 
using the stationary-point method for the a- model integral (|2.10|) . Their crucial observation 
was that on top of the trivial stationary point Q = A (expansion around which is just the 
usual perturbation theory), there exists another one, Q = kA, whose vicinity generates the 
oscillatory part of R(s) . (In the case of symplectic symmetry there exists an additional family 
of stationary points, see [fCT]]). The saddle-point approximation of Andreev and Altshuler is 
valid for s ^> 1; at 1 <C s -C g it reproduces the above results of Ref. flO| (we remind that 
the method of [ TO] works for all s <C g). The result of [ 1 L | has the following form 



R osc 


» 


COS27TS . . 

= 2vr 2 


R osc 


» 


COS27TS ~, . 
= 9^4 D ( S ) 


R Sp 


» 


_ cos2tt SjD1/2(s 



(2.34) 
(2.35) 

COS47TS ~, , , 

+ -^D\s) , (2.36) 



where D(s) is the spectral determinant 



D( 



6 m 



1 t— r / s z A z \ 1 



The product in Eq.( |2.37f) goes over the non-zero eigenvalues e M of the diffusion operator 
(which are equal to -Dq 2 for the cubic geometry). This demonstrates again the relation 
between R osc {s) and the perturbative part Q2.28D , which can be also expressed through D(s), 



Ras{s) - 1 -~2^ ds* (2 ' 38) 

In the high-frequency region s 3> g the spectral determinant is found to have the following 
behavior 

D(s) ~ exp { - — — — -. (—) } , (2.39) 

1 r(d/2)dsm(7rd/4) \g J 
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so that the amplitude of the oscillations vanishes exponentially with s in this region. 

Taken together, the results of 0] an d HJ] provide complete description of the deviations 
of the level correlation function from universality in the metallic regime g ^> 1. They 
show that in the whole region of frequencies these deviations are controlled by the classical 
(diffusion) operator governing the dynamics in the corresponding classical system. 



3 Statistics of eigenfunctions 



3.1 Eigenfunction statistics in terms of the supersymmetric a- 
model 

Within the RMT, the distribution of eigenfunction amplitudes is simply Gaussian, leading 
to the x 2 distribution of the "intensities" yi = iV|^>?| (we normalized yi in such a way that 
1^ ' 



(y) 
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V u (y) 

r°(y) 



-y/2 



(3.1) 
(3.2) 



Eq. ( p.2| ) is known as the Porter-Thomas distribution; it was originally introduced to describe 
fluctuations of widths and heights of resonances in nuclear spectra [12| . 



Recently, interest in properties of eigenfunctions in disordered and chaotic systems has 
started to grow. On the experimental side, it was motivated by the possibility of fabrication 
of small systems (quantum dots) with well resolved electron energy levels H, ||, H. Fluc- 



tuations in the tunneling conductance of such a dot measured in recent experiments [|94] , |95 
are related to statistical properties of wavefunction amplitudes |13|, 



97, m. When the 



electron-electron Coulomb interaction is taken into account, the eigenfunction fluctuations 
determine the statistics of matrix elements of the interaction, which is in turn important 
for understanding the properties of excitation and addition spectra of the dot IM [TOO , 
Furthermore, the microwave cavity technique allows one to observe experimentally spatial 
fluctuations of the wave amplitude in chaotic and disordered cavities JT^, [15], |To| . 

Theoretical study of the eigenfunction statistics in a rf-dimensional disordered system 
is again possible with use of the supersymmetry method [17, 19, 20]. The distribution 



function of the eigenfunction intensity u = \ip (i"o)| in a point ro is defined as 



V ^ = "7 (E <W°( r o)l 2 - u)5{E - E 



(3.3) 



The moments of V(u) can be written through the Green's functions in the following way 



■q-2 



<hKr )r> 



lim(2r ? )^ 1 (G| ? - 1 (r ,r )G' A (ro,ro)). 



2nuV v-+o 



(3.4) 



The product of the Green's functions can be expressed in terms of the integral over a super- 
vector field $ = (Si, xi, S 2 , X2), 



G q R \r , ro )G A (r ,r ) 



■2-q 



X 



(9-1) 

expjz J dr& (r)A 1/2 (E + ii]A - iT)A 1/2 $(r 



D^> D^(S 1 (r )Sl(r )y- 1 S 2 (T )S* 2 (r ) 



(3.5) 
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Proceeding now in the same way as in the case of the level correlation function (Sec. [2 . 1| ) , 
we represent the r.h.s. of Eq. (|3.5| ) in terms of a cr-model correlation function. As a result, 
we find[] 

(l^(ro)| 29 ) = ot/ lim^r 1 / DQQl-l b (r )Q 22 , bb (r )e- s ^\ (3.6) 
where S[Q] is the cr-model action, 



S[Q\ = ~ j d d r Str 



4 



(VQ) 2 - TTurjAQ 



(3.7) 



{(5 = 2 for the considered case of the unitary symmetry). Let us define now the function 
Y(Q )as 

Y(Q )= [ DQ(r)exp{-S[Q]} . (3.8) 

JQ(r )=Qo 

Here r is the spatial point, in which the statistics of eigenfunction amplitudes is studied. 
For the invariance reasons, the function Y(Q ) turns out to be dependent in the unitary 
symmetry case on the two scalar variables 1 < Ai < oo and — 1 < A2 < 1 only, which are the 
eigenvalues of the "retarded- retarded" block of the matrix Qq. Moreover, in the limit 77 — > 
(at a fixed value of the system volume) only the dependence on Ai persists: 

Y(Q ) = Y(X U A 2 ) -> Y a (27rur ] X 1 ) (3.9) 

With this definition, Eq. (|3.6|) takes the form of an integral over the single matrix Q$, 

(I^MI 29 ) = oT/ lim(27rz/r ? )"- 1 f DQ,Ql^ lbb Q^ bb Y{Q,) (3.10) 
ZV v-* J 

Evaluating this integral, we find 

(l^(ro)| 29 ) = ^q(q - 1) J duu^ 2 Y a {u). (3.11) 
Consequently, the distribution function of the eigenfunction intensity is given by [[17] 

V{u) = ^Y a (u) (U) , (3.12) 

where V is the sample volume. 

In the case of the orthogonal symmetry, Y(Qq) = F(Ai,A2,A), where 1 < Ai,A2 < 00 
and — 1 < A < 1. In the limit 77 — > 0, the relevant region of values is Ai ^> A2, A, where 

Y(Q ) ^ Y a (jrur]Xi) (3.13) 

The distribution of eigenfunction intensities is expressed in this case through the function 
Y a as follows |17| : 



1 r°° rP 



2v/2 d 2 f 00 dz 



itVu 1 / 2 du 2 







-Y a (z + u/2) (O) (3.14) 



z 

4 The first two indices of Q correspond to the advanced-retarded and the last two to the boson-fermion 
decomposition. 
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In the diffusive sample, typical configurations of the Q-field are nearly constant in space, 
so that one can approximate the functional integral (|3.8|) by an integral over a single su- 
permatrix Q. This procedure, which makes the problem effectively zero-dimensional and is 
known as zero-mode approximation (see Sec. |2.1| ), gives 



and consequently, 



Ya(z) 



V{u) 



-Vz 



(0,U) 



Ve 



-uV 



V c -uV/2 
2iru 



(U) 
(O) 



(3.15) 

(3.16) 
(3.17) 



which are just the RMT results for the Gaussian Unitary Ensemble (GUE) and Gaussian 
Orthogonal Ensemble (GOE) respectively, Eqs. (|3.1|) , (|3.2j ). 

Therefore, like in the case of the level correlations, the zero mode approximation yields 
the RMT results for the distribution of the eigenfunction amplitudes. To calculate deviations 
from RMT, one has to go beyond the zero-mode approximation and to evaluate the function 
Y a (z) determined by Eqs. (|3.8|), (|3T9| ) for a d- dimensional diffusive system. In the case of a 



quasi-lD geometry this can be done exactly via the transfer-matrix method, see Sec. |3.2j . For 
higher d, the exact solution is not possible, and one should rely on approximate methods. 
Corrections to the "main body" of the distribution can be found by treating the non-zero 
modes perturbatively, while the asymptotic "tail" can be found via a saddle-point method 
(see Sec. ^3] and Sec. £|). 

Let us note that the formula ( |3.8| ), ( |3.9| ) can be written in a slightly different, but com- 
pletely equivalent form |32|, [33| . Making in (|3.8| ) the transformation 



Q(r) 



Q(r) = V- 1 (r )Q(r)^(r ) 



where the matrix ^(r ) is defined from Q(r ) = V(r )AV 1 (r ), one gets in the unitary 

case 



(VQ) ; 



AQPi 



i,b 



(3.18) 



Y a {u) = DQexpl d d r Str 

jQ(r )=K [J 

where Pbb denotes the projector onto the boson-boson sector, and a similar formula in the 
orthogonal case. 

The above derivation can be extended to a more general correlation function representing 
a product of eigenfunction amplitudes in different points 



r {9} (ri, . . . ,r k ) = -L lC 9l (ri)||^ 2 (r 2 )| . . . \^(r k )\6(E - E a )\ 



(3.19) 



If all the points rj are separated by sufficiently large distances (much larger than the mean 
free path I), one finds for the unitary ensemble |T8 



ri,...,r k 



2V (qx + q 2 + ... + q k -l)U->o 
xJdQ gf i ; 66 1 (r 1 )Q 22)66 (r 1 )Qf 1)66 (r 2 ) . . . Qf 1)66 (r k )e 



lim(2vrz/r/) 9l+ - + ' 2 ^ 1 

-S[Q] 



(3.20) 



In the case of the quasi-lD system one can again evaluate Eq. (|3.20| ) via the transfer matrix 
method, while in higher d one has to use approximate schemes. The correlation functions 
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of the type ( |3.19| ) appear, in particular, when one calculates the distribution of the inverse 
participation ratio (IPR) P% = f d d r |^(r)|, the moments of which are given by Eq. (|3.19|) 
with qi = q 2 = . . . = qu = 2. We will discuss the IPR distribution function below in Sec. p. 2. 4 



and |3.3.3| . The case of k = 2 in Eq. (|3.19|) corresponds to the correlations of the amplitudes 



of an eigenfunction in two different points; we will discuss such correlations in Sec. |3.3.3j and 
in Sec. |4.1.1| (where they will describe the shape of an anomalously localized state). 

3.2 Quasi-one-dimensional geometry 
3.2.1 Exact solution of the cr-model 

In the case of quasi-lD geometry an exact solution of the cr-model is possible due to the 
transfer-matrix method. The idea of the method, quite general for the one- dimensional 
problems, is in reducing the functional integral ( |3.8| ) or ( |3.20|) to solution of a differential 



equation. This is completely analogous to constructing the Schrodinger equation from the 
quantum-mechanical Feynman path integral. In the present case, the role of the time is 
played by the coordinate along the wire, while the role of the particle coordinate is played 



by the supermatrix Q. In general, at finite value of the frequency rj in Eq. (377) (more 
precisely, rj plays a role of imaginary frequency), the corresponding differential equation is 
too complicated and cannot be solved analytically |J. However, a simplification appearing in 
the limit rj — > 0, when only the non-compact variable Ai survives, allows to find an analytical 
solution (1TH] of the ID a-moc 



There are several different microscopic models which can be mapped onto the ID super- 
matrix cr-model. First of all, this is a model of a particle in a random potential (discussed 
above) in the case of a quasi-lD sample geometry. Then one can neglect transverse variation 
of the Q-field in the cr-model action, thus reducing it to the ID form [|101| , [|. Secondly, 
this is the random banded matrix (RBM) model [ |102| , |T7| , |l~8|j which is relevant to various 
problems in the field of quantum chaos [|103j , |104|1 . In particular, the evolution operator of a 
kicked rotor (paradigmatic model of a periodically driven quantum system) has a structure 
of a quasi-random banded matrix, which makes this system to belong to the "quasi- ID uni- 
versality class" | IS, |105| | . Finally, the Iida- Weidenmiiller-Zuk random matrix model p3| of 



the transport in a disordered wire (see Section ^| for more detail) can be also mapped onto 
the ID cr-model. 

The result for the function Y a (u) determining the distribution of the eigenfunction inten- 
sity u = |^ 2 (r )| reads (for the unitary symmetry) 

Y a (u) = W (1) (uA£,t + )W {1) (uA^t_) . (3.21) 

Here A is the wire cross-section, £ = 2t[i>DA the localization length, r + = L + /£, r_ = 
with L_|_, L_ being the distances from the observation point ro to the right and left edges 
of the sample. For the orthogonal symmetry £ is replaced by £/2. The function W^'(z,t) 
satisfies the equation 

and the boundary condition 

W (l \z,0) = 1 . (3.23) 

5 Let us stress that we consider a sample with the hard-wall (not periodic) boundary conditions in the 
logitudinal direction, i.e a wire with two ends (not a ring). 
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The solution to Eqs.( |3.22| ), ( |3.23| ) can be found in terms of the expansion in eigenfunctions 



of the operator z 2 -^ — z. The functions 2z l l 2 Ki tl {2z l l 2 ), with K v {x) being the modified 
Bessel function (Macdonald function), form the proper basis for such an expansion [106], 
which is known as the Lebedev-Kontorovich expansion; the corresponding eigenvalues are 
-(1 + /i 2 )/4. The result is 

WV>{z,t) = 2Z 1 / 2 Ik^I 2 ) + - r^-^sinh^^(2. 1 / 2)e -^A . (3 .24) 



The formulas (|3.12j) , fl3.14|) , (|3.21|) , (|3.24j ) give therefore the exact solution for the eigen- 



function statistics for arbitrary value of the parameter X = L/£ (ratio of the total system 
length L = L + + L_ to the localization length). The form of the distribution function V(u) 
is essentially different in the metallic regime X <C 1 (in this case X = 1/g) and in the 
insulating one X ^> 1. We will discuss these two limiting cases below, in Sections [3.2.3| and 



3.2.4 respectively. 



3.2.2 Global statistics of eigenfunctions 

The multipoint correlation functions ( |3.19|) determining the global statistics of eigenfunctions 
can be also computed in a similar way. Let us first assume that the points lie sufficiently 
far from each other, |r$ — rj\ ^> I. We order the points according to their coordinates Xi 
along the wire, < x\ < . . . < Xk < L, and define ti — Ti = U — U-ii r i = t\- Then we 
find from Eq. ( |3.20| ) for the unitary symmetry 



vr«( ril ...,r k ) 



q x \ ...q h \ 



(q 1 + ...q k -2)\(ZA)<H+~+ik- 1 

k 

X 



poo 

/ dzz^W^&X- VrJfWfen,..,^) , (3.25) 

J s =l 

where the functions W^(z; r 1; . . . , t 8 ) are defined by the equation (identical to Eq. ( |3.22| )) 

dw^&n,...,^) ( 2 d 2 



dr. 



s 



(^z 2 ^-z^W^(z;n,...,T s ) (3.26) 



and the boundary conditions 

W {S \ Z] r 1; r 2 , . . . , T.-X, 0) = z^W^iz; n,r 2 ,..., r._i) . (3.27) 

Solving these equations consecutively via the Lebedev-Kontorovich transformation, one can 
find all the correlation functions VT* 9 ' (rj , . . . , r^) (in the form of multiple integrals over jj s ). 
This will be in particular used in Section [4 . 2 . 1| , where we will study the joint distribution 
function of the wave function intensities in two points (k = 2). 

We show now that the correlation functions ( 3.25| ) allow to represent the statistics of 
eigenfunction envelops in a very compact form. Making the substitution of the variable 
z = e e and defining the functions W^ s '(9; Ti, . . . , t s ) = z^ 1 ^ 2 W^ s \z; n, . . . , r s ), we can rewrite 
( P-26| ) in the form of the imaginary time Schrodinger equation, 

- ^_ = wW , h = -% + <? + j (3.28) 
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with the boundary conditions 

W^{9; n, . . . , r._i, 0) = e'- 1 "^'-^; n, . . . , r._i) , 

W^W^O) = e- e / 2 . (3.29) 

This allows to rewrite Eq. (|3.25|) as a matrix element, 

vr {q} (~ - " n^-Qk 



Ji, • • • ,r k 



(q 1 + ...q h -2)l(£A)*+-+*- 1 
x ( e -»/2| e -(^-E^ 1 r.)H eff .e e -r.^ eg ,_ 1 fl ^ le^ 2 ) (3.30) 

(these transformations are completely analogous to those performed by Kolokolov in ||107|1 , 
where the eigenfunction statistics in the strictly ID case was studied). Furthermore, the 
matrix element can be represented as a Feynman path integral, 

xAf [ D9(t) exp (- f dt (-9 2 + e e + 

Je(o)=e l ,e(x)=e r I J V 4 V J 

Xe 9if(ti)+926>(t2)+...+9fc6l(tfc) ^ (3.31) 

with A/" being the normalization constant, A/"^ 1 = J DO exp{ — | J dt6 }. The quantum 
mechanics defined by the Hamiltonian ( |3.28 ) (or, equivalently, by the path integral ( |3.31| )) 
is known as Liouville quantum mechanics [ 108| , p. 09)1 ; the corresponding spectral expansion 
is obviously equivalent to the Lebedev-Kontorovich expansion. 

Inserting here the decomposition of unity, 1 = J dwS(X~ 1 J dte e — w) and making a shift 
9 — > 9 + In w, we get 

X e X p|-i y a | e ffi«(*i)+-+9«.<'(**) < j (V 1 ^ dte e - lj . (3.32) 

According to ( |3.32j ), the eigenfunction intensity can be written as a product 

^(r) = $(r)*(t) , (3.33) 

where $(r) is a quickly fluctuating (in space) function, which has the Gaussian Ensemble 
statistics, (|$ 2<? |) = q^/V g , and fluctuates independently in the points separated by a distance 
larger than the mean free path. The function \l/(t) determines, in contrast, a smooth envelope 
of the wave function. Its fluctuations are long-range correlated and are described by the 
probability density 

V{6(t) = ln^ 2 (t)} =A/V x/ V M1 ^exp j-i J dt9 2 ^ . (3.34) 

The above calculation can be repeated for the case, when some of the points ri lie closer 
than / to each other. The result (|3.33|) , (|3.34j) is reproduced also in this case, with the function 
$(r) having the ideal metal statistics given by the zero- dimensional cx-model. This statistics 
||110| , |111| , |112j| is Gaussian and is determined by the (short-range) correlation function 

V(<S>*(rMr>)} = k 1 J 2 (\r-r'\) ] (3.35) 
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see Eq. ( |3.71| ) below. 

The physics of these results is as follows. The short-range fluctuations of the wave 
function (described by the function $(r)) have the same origin as in a strongly chaotic system, 
where superposition of plane waves with random amplitudes and phases leads to the Gaussian 
fluctuations of eigenfunctions with the correlation function ( |3.35| ) and, in particular, to the 
RMT statistics of the local amplitude, (|$ 2<7 |) = q\/V q . The second factor ^(t) in the 
decomposition (|3.33| ) describes the smooth envelope of the eigenfunction (changing on a 
scale ^> /), whose statistics is given by (|3.34|) and is determined by diffusion and localization 
effects. 

Let us note that in the metallic regime, X <C 1, the measure (|3.34|) can be approximated 

as 



V{9(x)} = exp 




(3.36) 



We will see in Section |J while studying the statistics of anomalously localized states in 
d > 1 dimensions, that the probability of appearance in a metallic sample of such a rare 
state with an envelope e e (r) is given (within the exponential accuracy) by the d- dimensional 
generalization of ( |3.36|) (see, in particular, Eqs. (^.10|) , fl4.77|) ). 

Finally, we compare the eigenfunction statistics in the quasi-lD case with that in a strictly 
ID disordered system. In the latter case, the eigenfunction can be written as 



i) 1D (x) 



— cos(kx + 5)^(x) 



(3.37) 



where \l/(x) is a smooth envelope function. The local statistics of iPid{x) (i-e. the moments 
( fj.4| )) was studied in ||113|| , while the global statistics (the correlation functions of the type 
( |3.19| )) in ||107| . Comparing the results for the quasi-lD and ID systems, we find that the 
statistics of the smooth envelopes \1/ is exactly the same in the two cases, for a given value 
of the ratio of the system length L to the localization length (equal to fiirvAD in quasi-lD 
and to the mean free path I in ID). In particular, the moments r^(r) = (|^ 2,? (r)|) are found 
to be related as 



A q T 



(?) 

Q1D 



qV 



ID ' 



(3.38) 



(2g-i) 

where the factor q\ 2 /(2q — 1)!! represents precisely the ratio of the GUE moments, (l^ 29 )) = 
q\/V q , to the plane wave moments, ((2/V) q cos 2q {kx + 5)) = (2g — l)\\/q\V q . For the case of 
the orthogonal symmetry of the quasi-lD system, this factor is replaced by g!. Equivalence 
of the statistics of the eigenfunction envelopes implies, in particular, that the distribution of 
the inverse participation ratio (IPR), 



Po 



dr 



(3.39) 



is identical in the ID [|114| , |115|| and quasi-lD [^, [18| cases (the form of this distribution 
in the localized limit L/C, 3> 1 is explicitly given in Sec. |3.2.4| below; for arbitrary L/£ the 
result is very cumbersome ||115|| ). 



3.2.3 Short wire 

In the case of a short wire, X = l/g < 1, Eq. (pl2|) , ( gig) , (|3~2TD , yield §§§ 

V (u) {y) = e -y [! + ax p _ Ay + y2) + . . j . y < x -i/2 (3 4Q) 
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-v/2 



2 ' 2 

p( l/ )( 2 /) = exp{- 2/ +f 2/ 2 X + ...} ; 
^ (0) (y) = 7feexp{i[-y + fy 2 X + ...]} 

p(v) 



x- 1 / 2 < ^ < X" 1 



cxp 



-2(5^/y~/X 



y>^ 



(3.41) 

(3.42) 
(3.43) 

(3.44) 



(a more accurate formula for the far "tail" ( 3.44j) can be found in Sec. |4.2. 1| , Eq. ( |4.74| )). 
Here the coefficient a is equal to a = 2[1 — 3L_L + /L 2 ]. We see that there exist three dif- 
ferent regimes of the behavior of the distribution function. For not too large amplitudes y, 
Eqs. (|3.40| ), (|3.41p are just the RMT results with relatively small corrections. In the inter- 
mediate range ( 3.42| ), ( |3.43| ) the correction in the exponent is small compared to the leading 
term but much larger than unity, so that V(y) 3> VnMriy) though lnV(y) — lnP^ MT (?/). 
Finally, in the large amplitude region, ( |3.44j ), the distribution function V(y) differs com- 
pletely from the RMT prediction. Note that Eq. (|3.44 ) is not valid when the observation 
point is located close to the sample boundary, in which case the exponent of (|3.44 ) becomes 
smaller by a factor of 2, see Sec. |4.2.3| . 

3.2.4 Long wire 

In the limit of a long sample, X = L/£ > 1, Eqs. ( gig ), ( PD , ( ET2ID , ( ggg) reduce to 



V {u \u) ~ 
V (0) (u) ~ 



X 1 2 (2 V / ^) + X 2 (2 v ^) 



L 

2i 2 AK l {2^Al) 



(3.45) 
(3.46) 



with £ = 2-nvAD as before. Note that in this case the natural variable is not y = uV, but 
rather uA£, since typical intensity of a localized wave function is u ~ 1/A£ in contrast to 
u ~ 1/V for a delocalized one. The asymptotic behavior of Eqs. (|3.45|) , (|3.46| ) at u ^> 1/A£ 
has precisely the same form, 



V{u) ~ exp(-2 / 3y / ttA£) 



(3.47) 



as in the region of very large amplitude in the metallic sample, Eq. ( 3.44 ). On this basis, 
it was conjectured in [18] that the asymptotic behavior ( |3.44j ) is controlled by the proba- 
bility to have a quasi-localized eigenstate with an effective spatial extent much less than £ 
("anomalously localized state"). This conjecture was proven rigorously in [[3(J where the 
shape of the anomalously localized state (ALS) responsible for the large-w asymptotics was 
calculated via the transfer-matrix method. We will discuss this in Section [| devoted to ALS 
and to asymptotics of different distribution functions. 

Distribution of the inverse participation ratio (IPR) is also found to have a simple form 
in the limit L > £ f2|, g5|: 



V(z) = 27r 2 ^(27r 2 ^ 4 -3A; 2 )e~ 7r2 

k=l 

4 d 



7T dz 



(3.48) 



k=i 
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Figure 1: Distribution function V(z) of the normalized (dimensionless) inverse participation 
ratio z = [/3 2 /(/3 + iy^vDA 2 P 2 in a long (L ^> £) quasi-lD sample. The average value is 
(z) = 1/3. From 



where z = irvDA 2 P2 in the unitary case and z = {irvDA 2 /3)P2 in the orthogonal case. (The 
second line in ( |3.48|) can be obtained from the first one by using the Poisson summation for- 
mula.) Therefore, the spatial extent of a localized eigenfunction measured by IPR fluctuates 
strongly (of order of 100%) from one eigenfunction to another. More precisely, the ratio of 
the r.m.s. deviation of IPR to its mean value is equal to 1/ Vo according to Eq. ( |3.48| ). The 
first form of Eq. fl3.48|) is more suitable for extracting the asymptotic behavior of V(z) at 
z ^> 1, whereas the second line gives us the leading behavior of V(z) at small z < 1: 



Viz) 



An 4 ze- n2z , z > 1 

47r -l/2 z -7/2 e -l/, z ^ 1 



(3.49) 



Therefore, the probability to have atypically large or atypically small IPR is exponentially 
suppressed. The function V[z) is presented in Fig. |l|. 

The above fluctuations of IPR are due to fluctuations in the "central bump" of a localized 
eigenfunction. They should be distinguished from the fluctuations in the rate of exponential 
decay of eigenfunctions (Lyapunov exponent). The latter can be extracted from another 
important physical quantity — the distribution function V x (v), where 

v = {2^DA 2 ) 2 \^ 2 a {v 1 )^ 2 a {v 2 )\ 

is the product of the eigenfunction intensity in the two points close to the opposite edges of 
the sample r\ —>■ 0, —>■ L. The result is |23| , [18| 



p ( _l n ^ ) = F[ _ ( ^l n ^ )/2 X] 2(2 ^ )i/2 exp|-A (2X/ ^ + i n , ;) 
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u 



u 



r 2 [(3-n)/2] 
T(u) 



F°{u) 



uY 2 [(l-u)/2] 
-kY(u) 



(3.50) 



Therefore, In v is asymptotically distributed according to the Gaussian law with the mean 
value (— hav) = (2j3)X = L/PtcuAD and the variance var(— lav) = 2(— lav). The same 
log-normal distribution is found for the conductance and for transmission coefficients of a 
quasi-lD sample from the Dorokhov-Mello-Pereyra-Kumar formalism ||1 16| , [7^| (see end of 
Sec. §3). 

Note that the formula ( |3.50|) is valid in the region of v <C 1 (i.e. negative lav) only, which 



contains almost all normalization of the distribution function. In the region of still higher 
values of v the log-normal form of V(v) changes into the much faster stretched-exponential 
fall-off oc exp{— 2\f2(3v 1//4 }, as can be easily found from the exact solution given in [^, [L8[]. 
The decay rate of all the moments (v k ), k > 1/2, is four times less than (—Inv) and does 
not depend on k: (v k ) oc e~ x ^ 2/3 . This is because the moments (f fc ), k > 1/2 are determined 
by the probability to find an "anomalously delocalized state" with v ~ 1. 



3.3 Arbitrary dimensionality: metallic regime 
3.3.1 Distribution of eigenfunction amplitudes 

In the case of arbitrary dimensionality d, deviations from the RMT distribution V(y) for not 
too large y can be calculated p4|, ^5[ via the method described in Section |2|. Applying this 
method to the moments (|3.6| ), one gets 



(mr)\ 2q ) = £ 
mr)\ 2q ) - 



Vi 

(2?-l 



1 + -Kq(q-1) + 



(U) 



.[l + Kq (q-l) 



(O) 



where k = Il(r,r). Correspondingly, the correction to the distribution function reads 



v{v) 

V(y) 



K 



l + -(2-4y + y 2 ) + 



k /3 y 2 



(U) 
(O) 



(3.51) 
(3.52) 

(3.53) 
(3.54) 



Deviations of the eigenfunction distribution function V(y) from its RMT form are illustrated 
for the orthogonal symmetry case in Fig. fj. Numerical studies of the statistics of eigenfunc- 



tion amplitudes in weak localization regime have been performed in Ref. ||117|| for the 2D 
and in Ref. [|1 18f for the 3D case. The found deviations from RMT are well described by the 
above theoretical results. Experimentally, statistical properties of the eigenfunction intensity 
have been studied for microwaves in a disordered cavity |15| . For a weak disorder the found 
deviations are in good agreement with (|3.54|) as well. 



In the quasi-one-dimensional case (with hard wall boundary conditions in the longitudinal 
direction), the one-diffuson loop II(r,r) is equal to 



k = Il(r, r) 



x 



1 



< x < L 



(3.55) 



so that Eqs. ( |3.53| ), (|3.54| ) agree with the results ( |3.40|) , (|3.41| ) obtained from the exact 
solution. For the periodic boundary conditions in the longitudinal direction (a ring) we have 
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Figure 2: Distribution V(y) of the normalized eigenfunction intensities y = V\ip 2 (r)\ in the 
orthogonal symmetry case. The dashed line shows the RMT result, Eq. ( |3.2| ), while the full 
line corresponds to Eq. fl3.54j ) with k = 0.4. 
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k = 1/Qg. In the case of 2D geometry, 



n< 



r,r) = — In j, (3.56) 
ng I 

with g = 2tivD. Finally, in the 3D case the sum over the momenta n(r, r) = (ttuV)^ 1 Ylq(Dq? 
diverges linearly at large q. The diffusion approximation is valid up to q ~ the corre- 
sponding cut-off gives Il(r, r) ~ Xjli^vDl = g~ l (L/l). This divergency indicates that more 
accurate evaluation of Il(r, r) requires taking into account also the contribution of the bal- 
listic region (q > Z -1 ) which depends on microscopic details of the random potential. We 
will return to this question in Sec. |3.3.4 . 



The formulas ( |3.53| ), ( |3.54|) are valid in the region of not too large amplitudes, where the 
perturbative correction is smaller than the RMT term, i.e. at y <C kT 1 ! 2 . In the region of 
large amplitudes, y > k~ x I 2 the distribution function was found by Fal'ko and Efetov [33|j , 
who applied to Eqs. (|3.12|) , (|3.14|) the saddle-point method suggested by Muzykantskii and 
Khmelnitskii P9[ . We relegate the discussion of the method to Section [| and only present 



the results here: 



W-exp^(- y + ^ + ...)}x{ ^ f Q ] 
V(y)~exp{-^-\n d (Ky) 



K- 1/2 < V < 



y ^ 



(3.57) 
(3.58) 



Again, as in the quasi-one-dimensional case, there is an intermediate range where a correction 
in the exponent is large compared to unity, but small compared to the leading RMT term 
[Eq. ( p.57| )1 and a far asymptotic region ( |3.58| ), where the decay of V(y) is much slower 
than in RMT. In the next section we will discuss the structure of anomalously localized 
eigenstates, which are responsible for the asymptotic behavior (|3.44|), (|3.58|). 



-i 



3.3.2 2D: Weak multifractality of eigenfunctions 

Since d = 2 is the lower critical dimension for the Anderson localization problem, metallic 
2D samples (with g ^> 1) share many common properties with systems at the critical point 
of the metal-insulator transition. Although the localization length £ in 2D is not infinite (as 
for truly critical systems), it is exponentially large, and the criticality takes place in the very 
broad range of the system size L £. 



Multifractality: basic definitions. The criticality of eigenfunctions shows up via their 
multifractality. Multifractal structures first introduced by Mandelbrot [|119|| are characterized 
by an infinite set of critical exponents describing the scaling of the moments of a distribution 
of some quantity. Since then, this feature has been observed in various objects, such as the 
energy dissipating set in turbulence 



120, 121 



strange attractors in chaotic dynamical 
systems ||123| , |124| , |125| , l26 |, and the growth probability distribution in diffusion-limited 
aggregation [127|, pH [L2!|; see Ref. |3| for a review. 

The fact that an eigenfunction at the mobility edge has the multifractal structure was 
noticed for the first time in [51j|, though the underlying renormalization group calculations 
were done by Wegner several years earlier [|5{J. For this problem, the probability distribution 
is just the eigenfunction intensity |V' 2 ( r )| an d the corresponding moments are the the inverse 
participation ratios, 

r d d r\ip 2q {r)\ . (3.59) 
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The multifractality is characterized by the anomalous scaling of P q with the system size L, 

P q otL- D «( q -V =L-< q \ (3.60) 

with D q different from the spatial dimensionality d and dependent on q. Equivalently, the 
eigenfunctions are characterized by the singularity spectrum f(a) describing the measure 
jjw of the set of those points r where the eigenfunction takes the value |V' 2 ( r )| ^ L~ a - The 
two sets of exponents r(q) and f{ct) are related via the Legendre transformation, 

T (q) = qa- f(a) ; f'(a) = q ; r'(q) = a . (3.61) 

For a recent review on multifractality of critical eigenfunctions the reader is referred to 



55, 131 



Multifractality in 2D. We note first that the formulas ( [37521) for the IPR's with 



< K.-V2 



can be rewritten in the 2D case (with (|3.56|) taken into account) as 



TO 

pRMT 
1 



1 



1 9(5-1) 



(3.62) 



where P,f MT is the RMT value of P q equal to glL" 2 ^" 1 ) for GUE and (2q 
GOE. We see that (|3.62| ) has precisely the form ( p.60|) with 



D n 



2 - 



(3ng 



1)11^-2(5-1) for 



(3.63) 



As was found in |32], |33|, the eigenfunction amplitude distribution ( |3.57D , (|3.58| ) leads to 



the same result ( p. 63 ) for all q <C Zftirg. Since deviation of D q from the normal dimension 
2 is proportional to the small parameter 1/irg, it can be termed "weak multifractality" (in 
analogy with weak localization). The result ( |3.63| ) was in fact obtained for the first time by 
Wegner [|50] via the renormalization group calculations. 

The limits of validity of Eq. ( |3.63| ) are not unambiguous and should be commented here. 
The singularity spectrum f(a) corresponding to ( |3.63|) has the form 



/(«) 



2-^(2 



— a 



so that f(a± = 0) for 



a± 



1 ± 



(3.64) 



(3.65) 



If a lies outside the interval (a_, a+), the corresponding f(a) < 0, which means that the 
most likely the singularity a will not be found for a given eigenfunction. However, if one 
considers the average (P q ) over a sufficiently large ensemble of eigenfunctions (corresponding 
to different realizations of disorder), a negative value of f{o) makes sense (see a related 
discussion in ||132| , |133|| ). This is the definiton which was assumed in |32|, |33| where Eq. ( |3.63| ) 
was obtained for all positive q 2/3tt g. 

In contrast, if one studies a typical value of P q , the regions a > a + and a < «_ will 
not contribute. In this case, Eq. (|3.63| ) is valid only within the interval g_ < q < q + with 
q± = ±(2/?7T(yf) 1//2 ; outside this region one finds ||134|, |135 



r(q) = D q {q - I) 



got- , 

qa + , 



q > qi, 
q<q- 



(3.66) 
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Therefore, within this definition the multifractal dimensions D q saturate at the values a + and 
a_ for q — » +00 and q — > —00 respectively. This is in agreement with results of numerical 



simulations [53 



3.3.3 Correlations of eigenfunction amplitudes and fluctuations of the inverse 
participation ratio 

In this subsection, we study correlations of eigenf unctions in the regime of a good con- 
ductor |25|, |27], |136| , |137|| . The correlation function of amplitudes of one and the same 
eigenfunction with energy E can be formally defined as follows: 



a 



;r 1; r 2 , E) = (\Mn)Mr2)\ 2 ) E = ^{J2 \Mn)Mr2)\ 2 5(E - e k ) 



(3.67) 



An analogous correlation function for two different eigenfunctions is defined as 
a(ri,r 2 ,£,w) = <|^ fc (ri)^(r 2 )| 2 )^ 

+ v - ti) 



A 2 R-\uo) /j2 \Mn)Mr2)\ 2 S(E - e k )5(E 

\ k^l 



(3.68) 



where R{uo) is the two-level correlation function (|2.1|). To evaluate a(r 1; r 2 , E) and cr(r 1; r 2 , E, ui), 
we employ an identity 



2tt 2 A _1 a(n, r 2 , E)5(uo) + A- 2 R(u)a(n, r 2 , E, uo) 
= Re [{G R { ri , n, E)G A (r 2 , r 2 , E + uo) - G R {r u r u E)G R (r 2 , r 2 , E + uo))] 



(3.69) 



where G R ' A (r, r', E) are retarded and advanced Green's functions and R(u) is the non- 
singular part of the level-level correlation function: R(oo) = R(uo) + 5(uo/A). A natural 
question, which arises at this point, is whether the r.h.s. of Eq. fl3.69|) cannot be simply found 
within the diffuson- Cooper on perturbation theory ||. Such a calculation would, however, 
be justified only for uo ^> A (more precisely, one has to introduce an imaginary part of 
frequency: uo — > uo + iT, and require that T ^> A). Therefore, it would only allow to find 
a smooth in uo part of a(ri, r 2 , E, u) for uo ^> A. Evaluation of a(ri, r 2 , E), as well as of 
a(ri,r 2 , E,uo) at uo ~ A, cannot be done within such a calculation. For this reason the 
non-perturbative supersymmetry approach is to be used. 

The r.h.s. of Eq. ( |3.69| ) can be expressed in terms of the supermatrix a-model, yielding: 



2tt 2 



A -1 a(ri, r 2 , E)6(u) + A"V(ri, r 2 , E, uo)R{uo) 



(3.70) 



[1 - Re(Q fe 1 fe 1 (r 1 )Q 22 (r 2 )) 5 - k d (v x - r 2 )Re(g^ 2 (r 1 )g 2 6 1 (r 1 )) 5 ] 



where (■ ■ -)s denotes the averaging with the sigma-model action and 
kd(r) = (ttu)" 2 (lm.G (r)) 2 is a short-range function explicitly given by 



k d {r) = exp(-r//) 



JoiPFr), 
(Pft)~ 2 sin 2 pft, 



2D 
3D 



(3.71) 



We consider the unitary ensemble first; results for the orthogonal symmetry will be presented 
in the end. Evaluating the a-model correlation functions in the r.h.s. of Eq. (|3.70|) and 
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separating the result into the singular the singular (proportional to 5(u)) and regular at 
uj = parts, one can obtain the correlation functions a(ri, r 2 , E) and cr(ri, r 2 , E, uj). The 
two-level correlation function R{uj) entering Eq.( |3.70 ) was studied in Section ^. We employ 
again the method of [[HJ described in Section ^| to calculate the sigma-model correlation 
functions (Qlb( r i)Qbb( r 2)} s and (Q 1 b 2 (r 1 )QH(r 2 )} s for relatively low frequencies uj < E c . 
First, we restrict ourselves to the terms of order g~ l . Then, the result for the first correlation 
function reads as 



(Qll{*i)Qli{*2)) 



2i 



exp(iirs) sin its 2i 



{7TS ) ITS 



-n( ri ,r 2 ) 



(3.72) 



where s = uj/A + iO. The first two terms in Eq. ( |3.72|) represent the result of the zero-mode 
approximation; the last term is the correction of order g~ x . An analogous calculation for the 
second correlator yields: 



(Q 6 t(ri)Q^(r 2 )> 



-2{- + 

ITS 



l+i 



exp(27rsj sm7rs 



7TS 



n(r 1; r 2 ; 



(3.73) 



Now, separating regular and singular parts in r.h.s. of Eq. ( |3.70| ), we obtain the following 
result for the autocorrelations of the same eigenfunction: 



v 2 (\Mn)Mr2)\ 2 ) E - 1 = k d (r)[i + n(n, n)] + n(n, r 2 ) 

and for the correlation of amplitudes of two different eigenf unctions 

v 2 (\m^i)M^)\ 2 }e^ - 1 = fcd(r)n(n, n), k ± i 

In particular, for = r 2 we have 

V 2 (\M^)\ 2 )e,u, - 1 = Ski + (1 + <yn(r, r). 



(3.74) 



(3.75) 



(3.76) 



Note that the result (|3.74j) for rj = r 2 is the inverse participation ratio calculated above 
(section |3.3.1| ); on the other hand, neglecting the terms with the diffusion propagator (i.e. 
making the zero-mode approximation), we reproduce the result of Refs. [ 110| , |1 1 1| , |112|| . 

Eqs. ( |3.75| ), ( |3.76| ) show that the correlations between different eigenfunctions are rela- 
tively small in the weak disorder regime. Indeed, they are proportional to the small parameter 
n(r, r). The correlations are enhanced by disorder; when the system approaches the strong 
localization regime, the relative magnitude of correlations, Il(r, r) ceases to be small. The 
correlations near the Anderson localization transition will be discussed in Section |5[ 

Another correlation function, generally used for the calculation of the linear response of 
the system, 




<^(ri)^(n)^(r 2 )^(r 2 )> 



r k (ri)M^)M^Kr 2 )S(E - t k )5{E 

can be calculated in a similar way. The result reads 

V^rMM^)M^)^2)) E ^ = fc d (r) + n(n,r 2 ), 



(3.77) 



(3.78) 



As is seen from Eqs. (|3.74j) , ( p.75| ), (|3.78| ), in the 1/g order the correlation functions 
a(ri,r 2 ,_E) and 7(ri, r 2 , E, u) survive for the large separation between the points, r ^> I, 
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while a(ri,T2, E,u) decays exponentially for the distances larger than the mean free path 
/. This is, however, an artifact of the g" 1 approximation, and the investigation of the 
corresponding tails requires the extension of the above calculation to the terms proportional 
to g~ 2 . The correlator (Q&&( r i)^?&&( r a)) S gets the following correction in the g~ 2 order: 

*<Q£(ri)Qg(r a ))* = -fi + 2/ 4 + exp(2 CTg )/ 3 - 2 » eXp(2Mra) (/ 2 - / 3 ) 

TTS 

_ exp^)-! /2 + 3A _ 

where we defined the functions 

/i(r ls r 2 ) = n 2 ( ri ,r 2 ), 

/ 2 (r 1; r 2 ) = (2V)- i y'rfr[n 2 (r,r 1 ) + n 2 (r,r 2 )], 

h = V- 2 J dvdv'Yl 2 (v,v'), 

/ 4 (n,r 2 ) = V- 1 J M(r, ri )n(r,r 2 ). (3.80) 

Consequently, we obtain the following results for the correlations of different (k 7^ I) eigen- 
functions at r > I: 

v'UM^M^e,. - 1 = \{h -h- 2/4) 

As it should be expected, the double integral over the both coordinates of this correlation 
function is equal to zero. This property is just the normalization condition and should hold 
in arbitrary order of expansion in g^ 1 . 

The quantities / 2 , ^3, and / 4 are proportional to g~ 2 , with some (geometry-dependent) 
prefactors of order unity. On the other hand, fx in 2D and 3D geometry depends essentially 
on the distance r = |ri — r 2 |. In particular, for I ^ r ^ L 

1 , oL 



ln z - , 2D 



/l (r 1 ,r 2 )=n 2 (r 1 ,r 2 )^<; ^? r 

{A7T 2 uDr) 2 ' 

Thus, for I < r <C L, the contributions proportional to fx dominate in Eq. (|3.81|) , yielding 

V 2 {\M*i)Mr2)\ 2 )E,u,-l = ~Il 2 (ri,r 2 ) , k^l. (3.82) 

On the other hand, for the case of the quasi-lD geometry (as well as in 2D and 3D for r ~ L), 
all quantities fx, / 2 , fy, and / 4 are of order of l/g 2 . Thus, the correlator cr(ri, r 2 , E, uj) 
acquires a non-trivial (oscillatory) frequency dependence on a scale uj ~ A described by the 
second term in the r.h.s. of Eq.( j3.81|) . In particular, in the quasi-lD case the function f 2 — f'3 
determining the spatial dependence of this term has the form 

A = (t)] ■ (3 ' 83) 
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where -84(2) = x A — 2x 3 + x 2 — 1/30 is the Bernoulli polynomial. 

Let us remind the reader that the above derivation is valid for uo <C E c . In the range 
uo > E c the o"-model correlation functions entering Eqs. (|3.70|) can be calculated by means 
of the perturbation theory ||, yielding 



v 2 {\u*im*2)?)E, 



1 

+ 2 



V J (rf.;.(r 1 )*(ri)'fc(r2)'</',*(r 2 )>£,„ 
where il a ,(r 1 ,r 2 ) is the finite-frequency diffusion propagator 

— <Mri)0 g (r 2 



1 + Re j^rUrx,^) 
k d (r) + ReL T a; (r 1 ,r2), 



n w (ri,r 2 ) 



7TZ/ 



,-1 



q 



Dq 2 



iuo 



(3.84) 



(3.85) 



and the summation in Eq. (|3.85|) now includes q = 0. As was mentioned, the perturbation 
theory should give correctly the non-oscillatory (in uo) part of the correlation functions at 
uo A. Indeed, it can be checked that Eqs. (|3.84j ) match the results (|3.75| ), (|3.78|) of the 
non-perturbative calculation in this regime. Furthermore, in the 1/g order [which means 
keeping only linear in ILj terms in ( |3.84| ) and neglecting — iuj in denominator of Eq. (|3.85|) 1 
Egs-d^D, (HID reproduce the exact results ( g^BD , fl3Jp even at small frequencies uo ~ A. 
We stress, however, that the perturbative calculation is not justified in this region and only 
the supersymmetry method provides a rigorous derivation of these results. 

Generalization to a system with unbroken time reversal symmetry (orthogonal ensemble) 
is straightforward [|138j| ; in the 1/g-order Eqs. (|3.74j ), ( |3.75|) , and (|3.78|) are modified as follows: 



^ 2 <|^(ri)^(r 2 )| 2 >i=; = [1 + 2k d (r)] [1 + 2n(n, r 2 )] , 

^ 2 (|^.(ri)^(r 2 )| 2 )^ - 1 = 2fc rf (r)n(n, r 2 ) , 
V 2 (r k (ri)M^i)M^(^)}E,u J = k d (r) + [1 + k d (r)} n(n, r 2 ), k^l. 



(3.86) 

(3.87) 
(3.88) 



IPR fluctuations. Using the supersymmetry method, one can calculate also higher order 
correlation functions of the eigenfunction amplitudes. In particular, the correlation function 
(\ip^(r 1 )\\ip^(r 2 )\) e determines fluctuations of the inverse participation ratio (IPR) Pi = 



j c?r|^ 4 (r)|. Details of the corresponding calculation can be found in Ref.[25|; the result for 
the relative variance of IPR, £(-82) = var(P 2 )/ (P2) 2 reads 



drdr'_ r[2/ . A 32a d 



(3.89) 



with a numerical coefficient a d defined in Section || (see Eqs. ( p.25| ), ( |2.26| )). The fluctua- 
tions (|3.89|) have the same relative magnitude (~ 1/g) as the famous universal conductance 
fluctuations. Note also that extrapolating Eq. (|3.89| ) to the Anderson transition point, where 
g ~ 1, we find 5(P 2 ) ~ 1, so that the magnitude of IPR fluctuations is of the order of its 
mean value (which is, in turn, much larger than in the metallic regime; see Section |5|). 
Eq. (|3.89| ) can be generalized onto higher IPR's P q with q > 2, 

~ l) 2 f ^n 2 (r,rO = * q \:] fa \ (3.90) 



(P q ) 2 ~ (3 2 ^ ^ ' J V 2 v ' ' (3 2 g 2 
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so that the relative magnitude of fluctuations of P q is ~ q(q — l)/g. Furthermore, the higher 
irreducible moments (cumulants) ((P q )), n = 2,3, . . ., have the form 



{p q ) n 



(n-l)! 
2 

(n-l)! 







q(q- 1) 



dr x . . . dr n 



Il(r 1 ,r 2 )...Il(r a ,r 1 ) 



Tr 







<?(<?- i)n 



(3.91) 



where II is the integral operator with the kernel Il(r, r') /V. This is valid provided q 2 n <C 
2/3irg. Prigodin and Altshuler [|137fl obtained Eq. fl3.91|) starting from the assumption that 
the eigenfunction statistics is described by the Liouville theory. According to ( |3.91| ), the 
distribution function V(P q ) of the IPR P q (with q 2 /f3irg <C 1) decays exponentially in the 
region q(q - l)/g < P q /{P q ) - K 1, 



V{P q 



exp 



2 A g(g-l) 



(3.92) 



where e± is the lowest non-zero eigenvalue of the diffusion operator —DV 2 . 

The perturbative calculations show that the cumulants of the IPR's are correctly repro- 



duced (in the leading order in 1/g) if one assumes ||137|| that the statistics of the eigenfunction 
envelopes |^ 2 (r)| smoot h = e 6 ^ is governed by the Liouville theory (see e.g. ||139| , |140|| ) defined 
by the functional integral 



DO 5 



d d r 



■e — 1 J exp 



d d r(V6f 



(3.93) 



We will return to this issue in Sec. ^ where the asymptotics of the IPR distribution function 
will be discussed. We will see that these "tails" governed by rare realizations of disorder 
are described by saddle-point solutions which can be also obtained from the Liouville the- 
ory description ( |3.93|) . The multifractal dimensions (|3.63|) can be found from the Liouville 
theory as well [|139] , |140|| . It should be stressed, however, that this agreement between the 
supermatrix cr-model governing the eigenf unctions statistics and the Liouville theory is not 
exact, but only holds in the leading order in 1/g. 

Let us note that the correlations of eigenfunction amplitudes determine also fluctua- 
tions of matrix elements of an operator of some (say, Coulomb) interaction computed on 
eigenf unctions ipk of the one-particle Hamiltonian in a random potential. Such a problem 
naturally arises, when one wishes to study the effect of interaction onto statistical properties 
of excitations in a mesoscopic sample (see Section PJ). 



3.3.4 Ballistic effects 

Ballistic systems. The above consideration can be generalized to a ballistic chaotic sys- 
tem, by applying a recently developed ballistic generalization of the cx-model [[75], [76], [FT]. 
The results are then expressed [J79] in terms of the (averaged over the direction of velocity) 
kernel g(i"i, nj.; 112) of the Liouville operator K = t>i?nV governing the classical dynamics 
in the system, 



n B (ri,r 2 ) = J dnidn 2 g[r 1 ,n 1 ;r2,n 2 y, 

Kg( ri , n i; r 2 , n 2 ) = (th/)" 1 [<f(n - r a )<$(m - n 2 ) - V' 1 ] . (3.94) 
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Here n is a unit vector determining the direction of momentum, and normalization J dn = 1 
is used. Equivalently, the function ns(ri,r 2 ) can be defined as 

dt J dn 1 g(r 1 ,n 1 ,t;r 2 ) , (3.95) 

where g is determined by the evolution equation 

(^- + v F n 1 V 1 )g(T 1 ,n 1 ,t;T 2 ) = , t >0 (3.96) 



\dt 

with the boundary condition 

~g\t=o = M- X [5( ri - r 2 ) - V- 1 ]. (3.97) 



Eq.( |3.94| ) is a natural "ballistic" counterpart of Eq. (|2.22| ). Note, however, that n£(ri,r 2 ) 
contains a contribution n^ri,^) of the straight line motion from r 2 to ri (equal to 
1/(tvpf\ti — r 2 |) in 2D and to l/2(pi?|ri — r 2 |) 2 in 3D), which is nothing else but the smoothed 
version of the function fcd(|ri — r 2 |). For this reason, n(ri,r 2 ) in Eqs. ( |3.86| )-( |3.88| ) should 
be replaced in the ballistic case by Il(ri,r 2 ) = n B (r 1 ,r 2 )-ng ) (ri,r 2 ). At large distances 
| i*i — r 2 | ^> Xp the (smoothed) correlation function takes in the leading approximation the 
form 

\Z 2 a( ri , r 2 , E) = 1 + |n B (n, r 2 ). (3.98) 

A formula for the variance of matrix elements closely related to Eq. ( [3.98| ) was obtained in 
the semiclassical approach in Ref. [|141||. In a recent paper ||142|| a similar generalization of 



the Berry formula for (V'fc( r i)V'ifc( r 2)) was proposed. 

Eq. ( |3.98j ) shows that correlations in eigenfunction amplitudes in remote points are de- 
termined by the classical dynamics in the system. It is closely related to the phenomenon 
of scarring of eigenfunctions by the classical orbits [ 143 , 144|1 . Indeed, if ri and r 2 belong to 



a short periodic orbit, the function n#(ri,r 2 ) is positive, so that the amplitudes \ipk(ji)\ 2 
and \ipkij2)\ 2 are positively correlated. This is a reflection of the "scars" associated with 
this periodic orbits and a quantitative characterization of their strength in the coordinate 
space. Note that this effect gets smaller with increasing energy E of eigenfunctions. Indeed, 
for a strongly chaotic system and for |r x — r 2 | ~ L (L being the system size), we have in the 
2D case ns(ri,r 2 ) ~ \f/L, so that the magnitude of correlations decreases as E^ 1 ^ 2 . The 
function ns(ri,r 2 ) was explicitly calculated in Ref. J73] for a circular billiard with diffuse 
surface scattering (see Section H). 

Ballistic effects in diffusive systems. We return now to the question of deviations of the 
eigenfunction amplitude distribution from the RMT in a diffusive 3D sample. As was shown 
in Sec. |3.3.1| , such deviations are controlled by the parameter k = n(r,r), see Eqs. ( |3.53| ) 



( |3.56| ). The physical meaning of the parameter k is the time-integrated return probability, 



see Eq. ( |3.95| ) generalizing definition of n(r!,r 2 ) to the ballistic case. The contribution to 
this return probability from the times larger than the momentum relaxation time, t > r, is 
given by 

n diff (r,r) = ( 7 r^)- 1 J2 i^r 1 ■ 

|q|<l/2 

The sum over the momenta diverges on the ultraviolet bound in d > 2, so that the cut-off at 
q ~ l/l is required. This results in Eq. (|3.56|) in 2D and in n dlfT (r,r) ~ 1/(/cf/) 2 in the 3D 
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case. There exists, however, an additional - ballistic - contribution to Il(r,r), which comes 
from the times t shorter than the mean free time r. Diagrammatically, it is determined by 
the first term of the diffuson ladder contributing to Il(r, r) (that with one impurity line), i.e. 
by the probability to hit an impurity and to be rejected back after a time i <C r. Contrary 
to the diffusive contribution, which has a universal form and is determined by the value 
of the diffusion constant D only, the ballistic one is strongly dependent on the microscopic 
structure of the disorder. In particular, in the case of the white-noise disorder we find 

(dq) 1,1 



In — , 2D 



n ball (r,r) = <J ^[(t) "l 9 ^ 3D (3 ' 99) 



Avv 2 f t J q 2 kpl ' 

Note that the integrals over the momenta are again divergent at large q - precisely in the 
same way as in the diffusive region, but with different numerical coefficients - and are now 
cut-off at q ~ kp. 

The total return probability is given by the sum of the short-scale (ballistic) and long- 
scale (diffusive) contributions. It is important to notice, however, that the single-scattering 
contribution (|3.99|) should be divided by 2 in the orthogonal symmetry case, because the 
corresponding trajectory is identical to its time reversal. Thus, k = Il dlff (r,r) + Il ball (r,r) 
for (3 = 2 and « = IT diff (r,r) + (l/2)n ball (r, r) for (3 = 1. We see that for the white-noise 
random potential in 3D the return probability is dominated by the ballistic contribution, 
yielding k ~ 1/kpl. In the 2D case, taking into account of the ballistic contribution modifies 
only the argument of the logarithm in ( p.56|) . Furthermore, even in the quasi- ID geometry 



the non- universal short-scale effects can be important. Indeed, if we consider a 3D sample of 
the quasi- ID geometry (L x , L y <C L z ), the diffusion contribution will be given by Eq. ( |3.55|) , 
Il dlfr (r,r) ~ 1/g, while the ballistic one will be Il ball (r,r) ~ 1/kfl. Therefore, the diffusion 
contribution is dominant only provided g <C kpl. 

On the other hand, let us consider the opposite case of a smooth random potential with 
correlation length d^> \p. Then the scattering is of small- angle nature and the probability 
for a particle to return back in a time t <C r is exponentially small, so that II ball (r,r) can 
be neglected. Therefore, the return probability k in Eqs. (|3.51|) - (|3.54 ) is correctly given by 



the diffusion contribution, see Eq. ( |3.56|) for 2D and the estimate below it for 3D. Thus the 
corrections to the "body" of the distribution function are properly given by the a-model in 
this case. 



4 Asymptotic behavior of distribution functions and 
anomalously localized states 

In this section, we discuss asymptotics of distribution functions of various quantities char- 
acterizing wavefunctions in a disordered system. Asymptotic behavior of these distribution 
functions is determined by rare realizations of the disorder producing the states, which show 
much stronger localization features than typical states in the system. We call such states 
"anomalously localized states" (ALS). 

It was found by Altshuler, Kravtsov and Lerner (AKL) [^] that distribution functions 
of conductance, density of states, local density of states, and relaxation times have slowly 
decaying logarithmically-normal (LN) asymptotics at large values of the arguments. These 
results were obtained within the renormalization group treatment of the cr-model. The 
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validity of this RG approach is restricted to 2D and 2 + e-dimensional systems, with e C 1. 
On the other hand, the conductance, LDOS and relaxation times fluctuations in strictly 
ID disordered chains, where all states are strongly localized, were studied with the use of 
Berezinski and Abrikosov-Ryzhkin techniques ||145| , |146| , |147j , |148| , [L13|| . The corresponding 
distributions were found to be of the LN form, too. It was conjectured on the basis of this 
similarity |28, 113 , |149|| that even in a metallic sample there is a finite probability to find 
"almost localized" eigenstates, and that these states govern the slow asymptotic decay of the 



distribution functions. Similar conclusion [25] is implied by the exact results for the statistics 
of the eigenfunction amplitude in the quasi-one-dimensional case, which shows the identical 
asymptotic behavior in the localized and metallic regimes, see Section 3, Eqs. ( |3.44|) , ( |3.47|) . 

A new boost to the activity in this direction was given by the paper of Muzykantskii and 
Khmelnitskii [p9|, who proposed to use the saddle-point method for the supersymmetric o- 
model in order to calculate the long-time dispersion of the average conductance G(t). Their 
idea was to reproduce the AKL result by means of a more direct calculation. However, 
they found a different, power-law decay of G(t) in an intermediate range of times t in 2D. 
As was shown by the author |3D| (and then reproduced in |TT| within the ballistic cr-model 
approach), the far asymptotic behavior is of log-normal form and is thus in agreement with 
AKL. The saddle-point method of Muzykantskii and Khmelnitskii allowed also to study the 
asymptotic behavior of distribution functions of other quantities: relaxation times |29| , |30[ 
3~i| , eigenfunction intensities [32, 33| , local density of states |j4], inverse participation ratio 
|35| , |3"E|1 , level curvatures J57| etc. The form of the saddle-point solution describes directly 



the spatial shape of the corresponding anomalously localized state [|^, [36|] . 

We will consider the unitary symmetry ((3 = 2) throughout this section; in the general 
case, the conductance g in the exponent of the distribution functions is replaced by {(3/2)g 
(we will sometimes do it explicitly in the end of the calculation). 



4.1 Long-time relaxation 



In this subsection we study (mainly following Refs. |2"5| , |31]| ) the asymptotic (long-time) be- 
havior of the relaxation processes in an open disordered conductor. One possible formulation 
of the problem is to consider the time-dependence of the average conductance G(t) defined 
by the non-local (in time) current-voltage relation 



I(t) 



I dt'G(t - t')V(t') 

J —oo 



(4.1) 



Alternatively, one can study the decay low, i.e. the survival probability P s {t) for a particle 
injected into the sample at t = to be found there after a time t. Classically, P s (t) decays 
according to the exponential law, P s (t) ~ e _<//<D , where is the lowest eigenvalue of the 
diffuson operator —DV 2 with the proper boundary conditions. The time to has the meaning 
of the time of diffusion through the sample, and is the Thouless energy (see Section |]). 
The same exponential decay holds for the conductance G(t), where it is induced by the 
weak-localization correction. The quantities of interest can be expressed in the form of the 
cr-model correlation function 



G(t), P s (t) ~ / ^e- 



-iujt 



DQ{r)A{Q}e 



S[Q] 



(4.2) 



where S[Q] is given by Eq. ( |3.7|) with rj — > — 2iu. The preexponential factor v4{Q} depends 
on specific formulation of the problem but is not important for the leading exponential 
behavior studied here. 
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Varying the exponent in Eq. ( |4.2|) with respect to Q and ou, one gets the equations |29| 

2DV(QVQ) + zw[A, Q] = (4.3) 
Y / rfr Str(AQ) = t (4.4) 

[We assume unitary symmetry (/3 = 2); in the orthogonal symmetry case the calculation is 
applicable with minor modifications and we will present the result in the end.] Note that in 
fact u plays in fl4.2|) a role of the Lagrange multiplier corresponding to the condition (|4.4j) . 
Therefore, it remains 

i) to find a solution of Eq.fl4.3j) (which will depend on u); 



ii) to substitute it into the self-consistency equation ([4.4|) and thus to fix a; as a function 
of t- 



iii) to substitute the found solution Q t into Eq. (|4.2|) , which will yield 

P s (t) ~ exp {^Str J (Vg 4 ) 2 } (4-5) 
Note that Eq. (|4.3| ) is to be supplemented by the boundary conditions 

Pleads = A ( 4 - 6 ) 

at the open part of the boundary, and 

V " ( 5 1 insulator = ( 4 - 7 ) 

at the insulating part of the boundary (if it exists); V n denotes here the normal derivative. 
It is not difficult to show [[^J that the solution of Eq. Q4.3| ) has in the standard parametriza- 



tion the only non-trivial variable - bosonic "non-compact angle" < Q\ < oo; all other 

coordinates being equal to zero. As a result, Eq. ( PO ) reduces to an equation for #i(r) (we 
drop the subscript "1" below) 

V 2 # + ^sinhfl = , (4.8) 
the self-consistency condition fl4.4j ) takes the form 

ixv J d d r (coshO - 1) = t , (4.9) 

and Eq. ( (4.5Q can be rewritten as 

\nP s (t) = / d d riyef . (4.10) 



For sufficiently small times, 6 is small according to ( |4.9| ), so that Eqs. ( |4.8| ), ( |4.9|) can be 
linearized 

\7 2 9 + 2~f9 = ; 2 7 = iu/D ; (4.11) 

^ ! d d rd 2 = t. (4.12) 



is related to the eigenvalue Ai used in Sec. |3.1| as Ai = cosh^i 
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This yields 

0=^y 2 Mr), (4.13) 

where <\>\ is the eigenfunction of the Laplace operator corresponding to the lowest eigenvalue 
= I /Dtp. The survival probability ( |4.10| ) reduces thus to 

lnP s (t) = ^ J d d r9V 2 9 = -ixvD lx J d d r9 2 = -t/t D , (4.14) 

as expected. Eq. ( |4.14j ) is valid (up to relatively small corrections) as long as 9 <C 1, i.e. for 
tA <C 1 (A = 1/uV being the mean level spacing). To find the behavior at t > A , as well 
corrections at t < A -1 , one should consider the exact (non-linear) equation (|4.8|), solution 
of which depends on the sample geometry. 



4.1.1 Quasi-ID geometry 



We consider a wire of a length L and a cross-section A with open boundary conditions at 
both edges, 9(-L/2) = 9(L/2) = 0. Eqs. Q, Q take the form 

(4.15) 
(4.16) 



9" + 2 7 sinh 9 = 

•L/2 



/ (io:(cosh 9 — 1) = t/ixvA 

J-L/2 



From the symmetry consideration 9(x) = 9{—x) and #'(0) = 0, so that it is sufficient to 
consider the region x > 0. The solution of Eq. fl4.15|) reads 



x 



00 



dd 



ie(x) 2 a/7 (cosh # — cosh??) 
where 6q is determined by the condition 9(L/2) = yielding 

d-d 



(4.17) 



a/7(cos1i# — coshi?) 



(4.18) 



In the large-t region (tA 3> 1) we will have ^ > 1, and Eqs. ( |4.17|) , (|4.18 ) can be simplified 
to give 



2x 



9{x) ~ 9 I 1 - 

Substituting this into the condition fl4.16|) allows to relate 9 to t, 



Jo 



71 



'-tA9 



(4.19) 
(4.20) 

(4.21) 



Finally, substitution of Eqs. (|4.19[) and ( 4.21|) into ( |4 . 1 0|) yields the log-normal asymptotic 
behavior of P a (t) p|: 



lnP s (t) ~ -g In 2 



tA 



ln(tA) 
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tA > 1 



(4.22) 



with g = 2ttvAD/L being the dimensionless conductance. We remind that Eq. ( |4.22| ) has 
been derived for the unitary ensemble {(3 = 2); in the general case, its r.h.s. should be 
multiplied by [3/2. 

Equation ( |4.22|) has essentially the same form as the asymptotic formula for G(t) found 
by Altshuler and Prigodin ||148|| for the strictly ID sample with a length much exceeding the 
localization length: 

G(t) ~ exp |-iln 2 (t/r)| (4.23) 

If we replace in Eq. ( J4.23|) the ID localization length £ = / by the quasi-lD localization 
length £ = ftm/AD, we reproduce the asymptotics ( 4.22|) (up to a normalization of t in the 



argument of In 2 , which does not affect the leading term in the exponent for t — > oo). This is 
one more manifestation of the equivalence of statistical properties of smooth envelopes of the 
wave functions in ID and quasi-lD samples [|l^] (see Sec. ||]). Furthermore, agreement of the 
results for the metallic and the insulating samples demonstrates clearly that the asymptotic 
"tail" ( 4.22|) in the metallic sample is indeed due to anomalously localized eigenstates. 



As another manifestation of this fact, Eq. ( |4.22j ) can be represented as a superposition 
of the simple relaxation processes with mesoscopically distributed relaxation times fl28| : 

P s (t) ~ J dt^e-^Vm) ( 4 -24) 

The distribution function V{t^) then behaves as follows: 

~ exp{-g ln 2 (#A^)} ; ^ > — = t D (4.25) 



This can be easily checked by substituting Eq. (|4.25|) into Eq. ([4.24 ) and calculating the 



integral via the stationary point method; the stationary point equation being 

2gt^\n(gAt^ = t (4.26) 

Note that Thouless energy t~^ determines the typical width of a level of an open system. 
Therefore, the formula ( |4.25| ) concerns indeed the states with anomalously small widths tT 1 



in the energy space. 

The saddle-point solution 9(r) provides a direct information on the spatial shape of 
the corresponding ALS. This was conjectured by Muzykantskii and Khmelnitskii |29| and 



was explicitly proven in for the states determining the distribution of eigenfunction 
amplitudes, see Sec. |4.2.1| below. Specifically, the smoothed (over a scale larger than the 
Fermi wave length) intensity of the ALS is |'?/' 2 (r)| = J\f~ 1 e e (r), where Af is the normalization 
factor determined by the requirement J d r\ip 2 (r) | = 1. We get thus from Eq. ( |4.19| ) 



I^WI smoo th=^- 2 * W/L (4-27) 



with 



9 ~ In 0tA ln(tA)J ~ In (^gAt^ \n 2 (gAt^ . (4.28) 
Thus, the ALS, which gives a minimum to the level width £T , has an exponential shape 

(pg), o. 
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The saddle-point method allows also to find the corrections to Eq. ( |4.14j ) in the interme- 
diate region to *C t <C A -1 , where 9q <C 1 ||150|| . For this purpose, we expand cosh^o and 
cosh-$ in Eq. ([4 .IS ) up to the 4th order terms, which leads to the following relation between 
7 and 9 : 

(4.29) 



7 



71 

21? 



-fa 



Further, we substitute Eq. (§37]) (|4-16|) 



dt?(coshi? - 1) 



/ 2 ^7 (cosh O - coshtf) 2nuA 
and expand cosh $ in ( |4.3U| ) up to the 4th order terms. This gives the relation 



(4.30) 



t 



el 



71 



1 



1 



2-nvA U 8V27V 96 
Using Eqs. ( 4.15 ), ( 4.16| ), we can rewrite the action in the form 



S = ^J^- I dx{9') 2 = 2nvADL 1 {cos\i9 - 1; 



2Dt~i 



Expressing now 9q and 7 through t according to Eqs. ( 4.30[) , (|4.31|) , we find 



-\nP s {t) =S 



t 



1 - 



1 t 



2n 2 gt 



+ 



D 



(4.31) 



(4.32) 



(4.33) 



with to = I? fi?D. In the general case, g is replaced by ((3/2)g here. Eq. ( [4.33| ) is completely 



analogous to the formula (|3.42| ), ( |3.43|) for the statistics of eigenfunction amplitudes. It shows 
that the correction to the leading term —t/tjj in lnPg becomes large compared to unity at 
t > y/gtoi though it remains small compared to the leading term up to t ~ gto ~ A -1 . 

The result ( |4.33| ) was also obtained by Frahm ||151|| from rather involved calculations 
based on the equivalence between the ID a-model and the Fokker-Plank approach and 
employing the approximate solution of the Dorokhov-Mello-Pereyra-Kumar equations in the 
metallic limit. The fact that the logarithm of the quantum decay probability, In Ps(t), 
starts to deviate strongly (compared to unity) from the classical law, lnP s cl (t) = —t/to at 
t ~ ^Dy/g was observed in numerical simulations by Casati, Maspero, and Shepelyansky 
| 152|| . For related results in the framework of the random matrix model see Sec. |4.1.3| . 



4.1.2 2D geometry 

We consider now a 2D disk-shaped sample of a radius R with an open boundary. If the 
problem is formulated in terms of the conductance, we can assume the two leads attached to 
the disk boundary to be of almost semicircular shape, with relatively narrow insulating in- 
tervals between them. Then we can approximate the boundary conditions by using Eq. ( |4.6| ) 
for all the boundary. In fact, in view of the logarithmic dependence of the saddle point 
action on R (see below), the result does not depend to the leading approximation on the 
specific shape of the sample and the leads attached. With the rotationally invariant form of 
the boundary condition, the minimal action corresponds to the function 9 depending on the 
radial coordinate r only. We get therefore the radial equation 

9" + 9'/r + 2 7 sinh = 0; < r <R (4.34) 
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(the prime denotes the derivative d/dr) with the boundary conditions: 

9(R) = , (4.35) 
6'(0) = (4.36) 



The condition (|4. 36|) follows from the requirement of analyticity of the field in the disk center. 



Assuming that characteristic values of 6 satisfy the condition 6 ^> 1 (which corresponds 
to tA ^> 1), one can replace sinh^ by e e /2. Eq.( f4.34| ) can be then easily integrated, and its 
general solution reads: 

e e(r) = 2C l C 2 r c ^ 2 

with two integration constants C\ and C%. To satisfy the boundary condition (|4.36|) , we have 
to choose C\ — 2. Furthermore, the above assumption 6(0) ^> 1 implies that 2C2/7 ^ !• 
Therefore, the second boundary condition ( |4.35| ) is satisfied if C 2 — 8/7-R 4 , and the solution 



can be written in the form 

e d{r) ~[(r/R) 2 + 1 R 2 /8}- 2 (4.38) 

Using now the self-consistency equation ( f4.9|) one finds 7 = kn 2 vjt. Finally, calculating the 
action on the saddle point ( |4.38| ), we find 



- In P s [t)=S~ 8n 2 uD ln(tA) , (4.39) 

The above treatment is valid provided 9'(r) < l~ x on the saddle point solution, which is the 
condition of the applicability of the diffusion approximation (here / is the mean free path). 
In combination with the assumption 8(0) ^> 1 this means that 1 -C £ A (R/l) 2 . 

Now let us consider the region of still longer times, t 3> A~ l (R/l) 2 . In order to support 
the applicability of the diffusion approximation, we should search for a function 9(r) mini- 
mizing the action with an additional restriction 9' < Al^ 1 . Here A is a parameter of the order 
of unity, which can not be fixed within the diffusion approximation. We will see however 
that the saddle-point action depends on I through \a(R/l) only, and thus does not depend 
on A in the leading order, so that we can set A = 1. Since the derivative has a tendency to 
increase in the vicinity of r = 0, the restriction can be implemented [0 via replacing the 
boundary conditions ( (4.36Q by 6'(r*) = 0, where the parameter r* will be specified below. 
The solution reads now: 

e s(r) = (r/R) c 2 , . 

The function 8(r) is meant as being constant within the vicinity |r| < r* of the disk center. 
The condition 6' < l^ 1 yields r* ~ IC. It is important to note that the result does not depend 
on details of the cut-off procedure. For example, one gets the same results if one chooses 
the boundary condition in the form 9'(r*) = l/l. The crucial point is that the maximum 
derivative 6' should not exceed l/l. The constant C is to be found from the self-consistency 
equation (|4.9|) which can be reduced to the following form: 



R\ C 2t C 2 



rj Tt 2 vR 2 C-2 
Neglecting corrections of the ln(ln •) form, we find 



(4.41) 



In(tA) ^ ln(tA)_ 
ln(R/n) \n(R/l) 1 ' } 
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The action ( |4.10| ) is then equal to 



S ~ ir 2 vD(C + 2) 2 ln(R/n) ~ t^ ^M^ 



(4.43) 



Combining Eqs.( f4.39| ) and ( |4.43| ) and introducing the factor (3/2 for generality, we get 
thus the following long-time asymptotics of P s (t) (or G(t)) [[301 : 



P.(t) ~ (tA)- 2 ^ 
P s (t) ~ exp 



n/3g In 2 (t/gr) ] 
4 ln(i?/Z) J 



1 « tA < (i?//) 2 
tA > (i?//) 2 



(4.44) 
(4.45) 



where g = 2irvD is the dimensionless conductance per square in 2D and r is the mean free 
time. 



The far asymptotic behavior ( f4.45| ) is of the log-normal form and very similar to that 
found by AKL (see Eq. (7.8) in Ref. |2£|). It differs only by the factor 1/g in the argument 
of In 2 . It is easy to see however that this difference disappears if one does the last step of the 
AKL calculation with a better accuracy. Let us consider for this purpose the intermediate 
expression of AKL (Ref. |28J, Eq. (7.11)): 



PC 

G(t) oc -- / 
r Jo 



1 

Au T 



trh 



(4.46) 



where u ~ ^ n 7 ^ n ^ ne wea ^ localization region in 2D, which we are considering. Eval- 
uating the integral ( 4.46 ) by the saddle point method, we find 



G(t) 



exp 



exp 



1 



In 2 



2ut 



Au T 

,2 



7rg \n 2 (t/gr) \ 
A ln(R/l) J 



(4.47) 



where we have kept only the leading term in the exponent. Eq. (^4.47| ) is in exact agreement 
with Eq. ( f4.45| ) for (3=1 (AKL assumed the orthogonal symmetry of the ensemble). There- 
fore, the supersymmetric treatment confirms the AKL result and also establishes the region 
of its validity. It is instructive to represent the obtained results in terms of the superposition 
of simple relaxation processes with mesoscopically distributed relaxation times ts'- 



G(t),P 8 (t) 



dt 



(4.48) 



Equations ( |4.44j) , ( |L45| ) lead then to the following result for the distribution function V{t^) 

(U/tp)-**^ , t D < ^ < t D (f ) 2 



exp 



Trf3g In 2 fa/r) 
4 \n(B./l) 



h > t D (fy 



(4.49) 



where to — R 2 /D is the time of diffusion through the sample. 

The smooth envelope of the ALS corresponding to the intermediate region to <C t^ <C 
to{R/l) 2 has according to Eq. (|4.37|) the following spatial structure 



I^WI sm0 oth = ^" lee 



r) 



1 



8(3^^ [(r/R) 2 + FP/WDtf] 



2 ' 



(4.50) 
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Thus, this ALS has an effective localization length £gf ~ -^(^d/^) 1 ^ 2 , with the intensity 
decreasing as 1/r 4 outside the region of the extent £ e /. As to the ultra-long-time region, 
Eqs. ( |L40|) , Q4.42| ) indicate that now the effective localization length is given by if 1 defined 



as 



I® = Itl ; 7, - ¥t^# • (4-51) 



\n(R/l) 



The ALS intensity decays in a power-law manner for r > I* , with an exponent depending 
on tf 

4.1.3 Random matrix model 

Here we mention briefly the results on the quantum decay law obtained by Savin and Sokolov 



[153 1 within the RMT model. This will allow us to see the similarities and the differences 
between the diffusive systems and the random matrix model. The model describes a Hamil- 
tonian of an open chaotic system by a Gaussian random matrix coupled to M external 
(decay) channels. The found decay law has the form 

P 3 (t) ~ (1 + T w t/M)- M , (4.53) 

where = MTA/2n is a typical width of the eigenstate, with T characterizing the channel 
coupling (T = 1 for ideal coupling, see also Sec. |^). In this case, the product MT plays a 
role of the dimensionless conductance g (in contrast to the diffusive case where g is governed 
by the bulk of the system, here it is determined by the number of decay channels and the 
strength of their coupling). For not too large t (tAT <C 1), Eq. ( |4.53| ) yields the classical 
decay law, P s (t) ~ e~ tFw , with the corrections of the form 

\nP s (t) = -tV w (l - T w t/2M + ...), (4.54) 

which is similar to the results found for the diffusive systems (see e.g. Eq. ( |4.33|) ). At large 
t 3> (AT) -1 , the decay take the power-law asymptotic form 



lnP s (t) ~ -M\n{Y w t/M) . (4.55) 
4.1.4 Distribution of total density of states 

Here we discuss the contribution of ALS to the asymptotic behavior of the distribution 
function V[v) of the total density of states (DOS), 

u(E) = —lm[ d d rG R (r,r;E) (4.56) 
nV J 

(in the present subsection we denote the average DOS as u to distinguish it from the 

fluctuating quantity v(E)). A resonance state with an energy E and width t^ 1 gives a 
following contribution to v{E): 

vals{E) = -% = -UAu (4.57) 

7T V 7T 
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We expect that the asymptotic behavior of V{y) at sufficiently large 5v = v—v§ is determined 
by a single anomalously localized state. This yields pB] 



V{y) ~p(i^ = <^~J ~ ^Pi - ^ 1112 ^^/^ )} ' 



(4.5* 



in the quasi- ID case and 



V(u) ~V[t 



7T 5v \ 



(g5u/u )-^ , 

irg In 2 (Su / i/qAt) 



cxp 



ln(L//) 



(4.59) 



in the 2D geometry. The far LN asymptotic tail in Eq. ( |4.59|) is in full agreement with the 
RG calculation by Altshuler, Kravtsov, and Lerner We find also an intermediate power- 
law behavior, which could not be obtained from the study of cumulants in Ref. [28]. We 



note, however, that this power-law form is fully consistent with the change of the behavior 
of cumulants (({Su) n )) at n ~ irg discovered in |28| . 



4.2 Distribution of eigenfunction amplitudes 
4.2.1 Quasi-ID geometry 

The spatial shape of the ALS determining the asymptotics of the distribution function of 
eigenfunction intensities can be found [36] via the exact solution of the a-model. We define 

Q(u, r) 



Wl>« 



V(u) 



(4.60) 



where 



Q(u, r) 



1 



J2\Mr)\ 2 5(\M0)\ 2 -^(E-E o 



(4.61) 



and V(u) is the distribution function of u = \ip 2 (0) \ defined formally by Eq. ( |3.3| ). According 
to Eq. ( |4.60| ), (|^ 2 (r)|) u is the average intensity of an eigenstate, which has in the point 
r = the intensity u (which will be assumed to be atypically large). The exact result for 
V(u) is given in the form of the Lebedev-Kontorovich expansion by Eq. (|3.12|) , ([3.21 ), ( 3.24j) . 
Calculating the moments (|?/>(0)| 2 |?/>(r)| 2 ' 3 ), and restoring the function Q(u,r), we find for 
r > / 



Q(u,r) 







V^Adu 

where the function W^ 2 \z; r 1; r 2 ) satisfies the same equation, as 

dWW(z;n,T 2 ) 



(* 2 IH w(2)( 



(4.62) 



(4.63) 



with the boundary condition 



W {2 \z;0,r 2 ) = zW (1 \ 



z,r 2 ) 



(4.64) 
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The solution of Eqs. ( gjgg ), QOjg ) is 



W (2) (*,tl,t 2 ) = 2^ / d//6(//, r 2 )^(2 v ^)e- 1 ^ 



r 2 ) 



/xsinh(7r//) 
2^2 



*^(t)fW(t74,r 2 ) 



(4.65) 



Substituting here the formula ( |3 . 24| ) for (z,t 2 ) and evaluating the integral over z, we 
can reduce Eq. (|4.65| ) for r 2 ) to the form 



x sinh 



Li sinh(7r/i) 



7T//1 



16tt 2 



r i + z 



r 



1 + 2/i 



:i+/^ 2 ) + 



/i sinh(7r/i) f dfii \x 



r i + 2 



Li - III 



2tt 3 

2 2 

e ^ 



(4.66) 



In the opposite case r < / we find 

1 



(4.67) 



where the function Y a (u) was defined in Eq. (|3.9|) . This formula is valid for any sample, 
which is locally <i-dimensional. In the case of the quasi-lD geometry we get 

Q(u,r) = I {*,(,) («£ + |) - 1} K'KA,)^"KA,)] , r « I 

(4.68) 



Insulating sample (L ^> £) The distribution P('u) is given by (|3.45|) , (|3.46|) . The "tail", 
Eq. (|3.47 ), at « > 1/A£ corresponds to atypically large local amplitudes. Analyzing the 
general formula for (\ip 2 (r)\) in this case, we find the following spatial structure of the 
ALS with |V 2 (0)| = u: 



(|^ 2 (r)|) u = 2^ M -l/2^-3/2 r -3/2 e -r./4 5 



1/2 



1 - / JVf 



2 > 



(\r(r)\) u = l(^) 1/2 [l + 2^lMk d (r)] 



r 

/ < r < f 
r < I 



(4.69) 
(4.70) 

(4.71) 



We see from Eqs. (|4.69Q , ( |4.70|) , ( |4.71]) that the eigenfunction normalization is dominated by 
the region r ~ £ e /, where £ e / ~ \f^JuA <C £ plays the role of an effective localization length. 
In the region £ e f C r C ( the wave intensity falls down as 1/r 2 , and crosses over to the 
conventional localization behavior at r ^> £. Therefore, the appearance of an anomalously 
high amplitude |V' 2 (0)| = u 3> 1/A£ is not just a local fluctuation, but rather a kind of 
a cooperative phenomenon corresponding to existence of a whole region r < £ e j with an 
unusually large amplitude |?/> 2 (r)| = \^/uj^A ~ 1/A£ e j. 

Eq. (|4.71|) describes a sharp drop of the amplitude from |V' 2 (0)| — u to (|^ 2 (r)|) = 
\^fuf^A at r ~ I. This "quasi-jump" happens on a short scale r ~ /^(wA^) 1 ^ 2 ^ -1 )]. 
To understand the reason for it, let us recall that the above formulas represent the average 
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intensity |V' 2 ( r )| (under the condition |V' 2 (0)| = u). One can also study the fluctuations of the 
intensity. It turns out |}6f that in the region r < r < ^ the fluctuations are of usual GUE 
type superimposed on the envelope (|4.70|) . It is not difficult to understand that the quasi- 



jump has the same origin as the GUE-like fluctuations at r > r . One can ask, of course, 
why this short-scale fluctuation happens exactly in the center of the smooth ALS "bump" 
with a probability close to unity. The answer is as follows. We are studying the states with 
an anomalously large local intensity u, which is an exponentially rare event. There are two 
sources which may favor the formation of such a high intensity: i) formation of an ALS with 
a spatially non-uniform smooth envelope, and ii) short-scale GUE-like fluctuations. Both 
these mechanisms have exponentially small probabilities to produce an enhancement of the 
intensity by a large factor. The found configuration of (|V ,2 ( r )l)u (short-scale quasi-jump 
( [4.71| ) on top of the smooth configuration ( [4.70| )) represents just the optimal combination of 



the two mechanisms: 



« = l^ 2 (°)lsmooth " J(u) = \ (^) V2 • 2MO V2 (4-72) 
For arbitrary geometry of the sample, the magnitude of the "quasi-jump" J(u) is given 



according to Eq. ( fL67|) by 



J(u) ~ -u-^- In YJu) ~ -u-^- In V(u) (4.73) 
du du 



(in the quasi-lD case this reduces to J(u) = 2{uA£ s ) 1 / 2 as stated above). The formula (|4.73p 
can be reproduced (within the saddle-point approximation) by writing the quantity u as a 
product u = u s J of the smooth part u s and the local fluctuating quantity J, with the latter 
distributed according to V(J) = e~ J . 

Metallic sample (L <C £) The asymptotic behavior of the intensity distribution function 
has the same stretched-exponential form, as in the localized regime, see Eq. ( |3.44j ). More 
accurately (with the subleading factors included), this formula reads p6i 



^, , 16 [A^^/LlLZ \ 4 r— ix 2 i ( 
V(u) = -J^^^ exp )- A ^A + ^l 

it V u L 4L + \ 



La- 



4L_ \ L. 



+ ...]> ■ (4-74) 



Calculation of Q(u, r) shows |36] that the ALS intensity has for I < r L the same form 



( 4.7G|) , provided the condition V A |^ 2 (0)| ^> g is fulfilled. This condition, under which the 



asymptotic behavior Q3.44 ), ( 4.74j ) is valid, acquires now a very transparent meaning. This 



is just the condition that the effective localization length of an ALS, £ e / = ^/£/uA is much 
less than the sample size L. Indeed, £ ef /L = ^J^/uAL 2 = ^g/uV. 

Near the sample edges, r~I> £ e j, the form of the ALS intensity is slightly modified 
by the boundary of the sample, see |j36"| . Finally, the "quasi-jump" of {\ip 2 {j-)\) u at r -C / has 
the same form Q4.71 ) as in the insulating regime. 



Saddle-point method The saddle-point method of Muzykantskii and Khmelnitskii can 
be also applied to the problem of the statistics of eigenfunction amplitudes, as was done by 
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Fal'ko and Efetov |32|, [HJ. In this case, one should look for the saddle-point of the functional 
integral (|3.8|) determining the function Y a (u) (which in turn determines the eigenfunction 
statistics, see Eqs. (|3.12| ), ( |3.14j) ). The saddle-point is again parametrized by the bosonic 



non-compact angle 9(v) only, and the corresponding saddle-point equation has the form 

7iuDV 2 6 - ue e = . (4.75) 



It is similar to the equation ( |4.8|) of the long-time relaxation problem, but with different sign 
of the second term. Also, the boundary conditions have now a different form: 

9(0) = (4.76) 

and the condition (|4.7| ) at the boundary (since we consider now a closed sample). Alterna- 
tively, one can write Eqs. ( ^4.75| ), ( f4.76| ) in a slightly different form by shifting the variable 
9 + ln(uV) — > 9. Then u is removed from the saddle-point equation and from the action, 
but appears in the boundary condition. 

The action determining the distribution function V(u) is given by 



\nV(u) = S = J d d r 



-—{V9) 2 + ue e 



(4.77) 



The formula ( [4.77] ) acquires a very transparent meaning if we take into account what was 
written in Section [4.2.1| concerning the two factors contributing to the large amplitude 



|-?/> 2 (0)| = u. Firstly, this is the non-uniform smooth envelope oc e 9(r ) yielding 

2 e m 1 

W (0) (smooth ~~ j dd re e{r) ~ J dd re e(r) ' 

the corresponding weight is represented by the first term in the action ([4.771) . Secondly, 
these are the local Gaussian (GUE-like) fluctuations of the wavefunction amplitude, which 
should provide the remaining factor ("quasi-jump") 

j = - = u I d d re 0{r) 

l^ 2 (0)| smO oth 



the corresponding probability V(J) = e J reproduces the second term in the action ([4.77|) . 



In the quasi-lD case and under the condition u 3> gV 1 , the solution of Eq. (|4.75|) reads 



33, 34 



e e(r) = L ^ 78 s 



Comparing Eq. ( |4.78| ) with Eqs. (|4.70 ), we see that the saddle-point solution nicely repro- 
duces the average intensity of the ALS, (|V' 2 ( r )l)u for r > I, up to an overall normalization 
factor. Also, the form of V(u) found in the quasi-lD case by the saddle-point method [[33 



is in a very good agreement with the exact results presented above. This agreement does 
credit to the saddle-point method and allows to use the saddle-point configuration for char- 
acterizing the shape of ALS in higher dimensions and for other distribution functions, where 
the exact solution is not available. 
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4.2.2 2D geometry 

For a 2D disk-shaped sample of a radius L with the high amplitude point r = in the center 
of the disk, the saddle-point solution of Eqs. (|4.75| ), (|4.76|) is found to have the form |32], [33 



l J J 8(1 - I2) 2 71UD \h 

for U < r < L , (4.79) 



where the exponent < \i < 1 depends on u and satisfies the equation 

L\ 2 " 2-/i L 2 w 



(4.80) 



We are interested in the asymptotic region uL 2 ^> 7tzaD ln _1 (L/Z), where the distribution of 
the eigenfunction intensity is given by Eq. ( |3.58j ), and an ALS is formed. Then the exponent 
H can be approximated as 

2HL/L) (4 ' 81) 
The lower cut-off scale Z* appears in Eq. (|4.79|) for the same reason as in the long-time 



relaxation problem (Sec. 4.1.2 ), i.e. because of the restriction of the diffusion approximation 
on the momenta q of the cr-model field: q < Z -1 . It is determined by the condition 9('r)\ r =i t ~ 
Z" 1 , which yields Z* ~ fil. The corresponding asymptotic behavior of Y a {u) (and consequently 
of the intensity distribution function), which was already quoted in Sec. [3.3.1| , is 




Y a (u),V(u) ~ exp { 7i vD v ln(L//) 7 } (4.82) 

Normalizing the expression ( [4.79| ), we find that the average ALS density for r > Z* is 
equal to 

/ \ -2/t f ,o / \ 2-2n^ ~ 2 

u I r \ I Izu / r x 



«^l>- = irt^UJ r~ 5(1^5 U 1 > (4 ' 83) 



The saddle-point calculation assumes that 9(r) is constant for r < Z*, so that Eq. (|4. 83| ) 



gives (|V' 2 ( r )l)M — ^vBp, * n ^ n * s re gi° n - However, for very small r < Z the average intensity 
{\i^ 2 {r)\) u changes sharply, as we have seen in the quasi-lD case. Using Eq. ( |4.67| ) for Q(u, r) 
and Eq. ( p.!2|) for V(u) , we get 



(|^ 2 (r)|) u = ^2 = [i _ k 2 (r) + J(u)k 2 (r)W\r = Q\) u , r < Z* (4.84) 



According to (|4.73| ), ( |4.82| ), the height of the quasi-jump is given by 

J(«) - In Uu) = In In - j K 4, „fl„ , (4.85) 
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which is precisely the factor by which the value of (\ip 2 (r = Z*)|)« found above differs from 
(|^ 2 (0)|) u = u. Combining Eqs. ([183]), ([184]), and ([185]), we get 

(\^ 2 (r)\) u = — -[1 - k 2 {r) + h(r)J(u)} , r < U (4.86) 
J(u) 

Therefore, in the 2D case the ALS determining the asymptotics of the amplitude distribution 
function has the power-law shape ( fl83|) with the short-scale bump ( |4.86| ). 



4.2.3 States localized near the boundary. 

We assumed in the above calculations that the center of an ALS is located far enough from 
the sample edge. For a quasi-lD sample, this means that £ e j <C L + , L_. In the 2D case this 
implies that the distance from the observation point to the boundary is of the same order 
of magnitude in all directions, so that ln(L//) is defined without ambiguity. Here, we will 
consider briefly the role of ALS situated close to the boundary, when these conditions are 
violated g|. 

We first consider the quasi-lD geometry. Let us calculate the distribution function V(u) 
in a point located very close to one of the sample edges. Formally, this means that L_ <C £ e f. 
Then the function W {1) (uA£,t_) in Eq. ( ^21| ) can be approximated by unity, and we get 



V(u) = V 4 ^ 1/4 £- 1/2 «- 3/4 exp |- 2v ^+ |f! ^1 - + . . J | (4.87) 

We see therefore that close to the boundary the distribution V(u) has the asymptotic 
decay V(u) ~ exp{— 2v / wA£}, which is slower than in the bulk of the sample, V(u) ~ 
exp{— 4y/uA!;}. This means that if we consider the distribution V(u) averaged over the po- 
sition of the observation point, its asymptotic tail will be always dominated by contribution 
of the points located close to the boundary, V{u) ~ exp{—2^/uA£}. This could be already 

anticipated from Eq. ( fL74| ), where the factor exp j 2 ^- + -^-j j increases strongly with 

approaching one of the sample edges. The same tendency, but in a weaker form, is observed 
in Eqs. (|OoD , (]3%I|) , fl3T42p , Calculating the average intensity (|^ 2 (r)|) u of the cor- 



responding ALS, we find that at r > I the ALS spatial shape retain the form fl4.70p , with an 
additional overall factor of 2. At small r, Eq. (|4.71| ) is slightly modified: 

mr)\)u = [l + V^MUr)} (4-88) 

In 2D, we can consider a sample of the semicircular shape, with the observation point 
located in the center of the diameter serving as a boundary. The saddle-point solution then 
has exactly the same form ( [4.79[ ), and the ALS intensity is still given by Eq. Q4.83j ), with 



an additional factor 2. The asymptotic form of the distribution function V{u) contains an 
extra factor 1/2 in the exponent: 



~ exp < -_ ) (4.89) 

This result is expected to be applicable to any 2D sample of a characteristic size L, with a 
smooth boundary and the observation point taken in the vicinity of the boundary. 
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We see therefore, that, very generally, the probability of formation of an ALS with the 
center in a given point is strongly enhanced (via an extra factor 1/2 in the exponent), if this 
point lies close to the sample edge. This leads to the additional factor 1/2 in the exponent 
in the asymptotic form of the distribution V(u) near the boundary. 



4.3 Distribution of local density of states 

We again assume the sample to be open, as in the problem of the distribution of relaxation 
times, Sec. |4.1| . Then it is meaningful to speak about the statistics of the local density of 
states (LDOS) p(E,r) = (— l/7r)ImG^(r, r; E). In a metallic sample, the LDOS is a weakly 
fluctuating quantity, whose distribution V(p) is mostly concentrated in a narrow Gaussian 
peak |28, |155|| with mean value (p) = v and the variance var(p/z/) ~ k C 1, where n is 



the usual parameter of the perturbation theory k = II(r, r) = Y2 q l/irvVDq , which already 



appeared in Sec. p.3.1| . This close-to-Gaussian shape of V(p) holds however in the region 



\p/u — l\ <C 1 only. We will consider in contrast the "tails" of the distribution, p/y<l and 
pjv 3> 1, where a much slower decay of V{p) will be found. 

We begin by discussing on a qualitative level a relation of the asymptotics of V(p) to 
the behavior of the distribution functions considered above. Typically, in an open metallic 
sample the LDOS p(E,r) is given by a superposition of ~ 1/tpA = g adjacent levels, since 
their widths are of order of l/tp. However, we can expect that for p much larger than its 
average value u, the asymptotic form of V(p) is determined by a probability to have a single 
narrow resonance, which gives this value of LDOS p(E,r). The most favorable situation 
happens when the resonance is located around the point r in the real space and around the 
energy E in energy space. The LDOS provided by such a resonance is: 

PALS = \^(r)\^t (4.90) 

7T 

where t^ 1 is the resonance width. Thus, the optimal fluctuation should provide now a 
maximum to the product of the local amplitude u = \ip 2 (r)\ and the inverse level width t^, 
and the asymptotics of the distribution V(p) should be related to that of V(u) and V{t^). 
In particular, in the quasi-lD case, where the distribution V(t^), Eq. (|4.25| ), decays much 
more slowly than V(u), Eq. ( |3.47 ), one should expect the asymptotic behavior of V(p) to be 



mainly determined by Vit^). We will see below that this is indeed the case. 

Now we turn to a formal calculation. The distribution function P p (p) of LDOS can be 
expressed through the function F(Ai, A2) introduced in Sec. [JTT] as follows WQ, |2(], |322| , |323 




P{p) = 5(p -!) + -—{ / o d\ 1 Y(\ 1 ) I ^ ) } , (4.91) 

where 

F(A0 = f dX 2 ^^ (4.92) 



-1 Ai — A2 

and p is normalized by its mean value: pjv — > p. Let us note the symmetry relation found 



in [gg 

Pip- 1 ) = p 3 P(p) . (4.93) 

It follows from Eq. Q4.91D and is completely independent of a particular form of the function 
^(Ai, A2). Obviously, Eq. ( |4.93|) relates the small-p asymptotic behavior of the distribution 
Tip) to its large-p asymptotics. 
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Asymptotic behavior of V(p) was studied in [53[] via the saddle-point method supple- 
mented in the quasi-lD case by the exact solution. The saddle-point equation has for this 
problem a very simple form, 

V 2 # = , 

with the boundary conditions 

Pleads = > W = P/ 2 ■ 
Quasi-lD geometry. In the quasi-lD case, the solution reads 



- ::..-„„-• : (4.94) 



( P /u) l - r / L + , r > o 

(p/4 Hr|/i ~ , r < 

where, as before, L + and L_ are the distances from the observation point r = to the 
sample edges. This yields the asymptotics of the distribution function 

V(p) ~ exp (4 (V- + 41 ln 2 (p/^)l • (4.95) 



4 £- , 

Let us note that this asymptotic behavior of V(p) in an open sample is strongly different 
from the asymptotics of V{u) in a closed sample, Eq. ( |3.47|) . As was explained above, 
the difference originates from the fact that V(p) is essentially determined by V(t^). To 
demonstrate this explicitly, we put the observation point in the middle of the sample, L + = 
L- = L/2. The configuration (f4.94|) acquires then precisely the same form as the optimal 



configuration Q4.19] ) for the relaxation time t^, and the asymptotics (|4.95|) of V(p) is identical 
to that of P(t ), Eq. ( ggg ). 

The corresponding values of t$ and p are related as follows: 

-gAtf In 2 (g At ) = p/u . (4.96) 

7T 

Now we calculate the value of the local amplitude |^ 2 (0)| for an ALS corresponding to the 
configuration ( |4.94| ). First, its smoothed intensity is given by 



n mooth = X'^ = l ^(^ 2WL . (4.97) 

Second, the quasi-jump induced by the GUE-type fluctuations gives an additional factor, 
which is found according to Eq. ( |4.73| ) to be 

Jifi) = -p§- p ^V(p) = 2gln(p/u) . (4.98) 



Combining Eqs. ( |4.96| ) , Q4.97D , and (|4.98j) , we can compute the LDOS ( |4.90| ) determined by 
this resonance state: 

PALS (E,0) = |^ 2 (0)| smooth • J{ P ) ■ ^ = " 2gHpM ■ 2g ^ (pM = P (4-99) 

We have explicitly checked therefore that the LDOS p is indeed determined by a single 
ALS, smoothed intensity of which is given by Eq. ( |4.97|) . There are three sources of the 
enhancement of LDOS: i) amplitude of the smooth envelope of the wave function, ii) the 
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short-scale GUE "bump", and iii) the inverse resonance width. They are represented by the 
three factors in Eq. ( |4.99| ), respectively. 

The result ( 4.95| ) can be also obtained [34] from the exact solution of the a- model. It 
was also shown in |34[ that in the long wire limit (L>() the whole distribution function 
takes a log-normal form analogous to that found in ||113|| for the case of a strictly ID sample 
by the Berezinskii technique. 



,0(r) 



p / IX 



V 



2D geometry. In the 2D case, we introduce again (as in Sec. |4.1.2| , |4.2.2j ) the small 
cutoff, Z* . The saddle-point solution reads 



(4.100) 



where tf 1 



j p l, and 



7p 



ln(p/z/) 



ln(L//i p) ) 

The distribution V(p) has the following asymptotics 



V(p) ~ exp 



ti 2 uD In 2 p 



ln(L/Zi' 



(4.101) 



So, the LDOS distribution has (as in the quasi-lD case) the log-normal form, which is now 



similar to the distribution fl4.82|) of the eigenf unction intensities (while the distribution of 
relaxation times has an intermediate power-law regime). This result is in perfect agreement 



with the asymptotic behavior of V(p) found by the renormalization group method in [28 . 
As for all the other distribution functions studied, the relevant ALS decay in a power-law 



fashion, |V' 2 ( r )| oc 

(r//i P V 7p - Lik e in the quasi- ID case, one can explicitly verify |36| that 



a single state with the spatial shape determined by ( |4.100| ) indeed provides, by virtue of 
Eq. ( [OOP , the value of LDOS equal to p. 



4.4 Distribution of inverse participation ratio 

In this subsection, we study the asymptotics of the distribution function of the IPR P 2 , 
Eq. ( |3.39|) . We have already considered the fluctuations of P 2 in Sec. |3.3.3| . As was explained 



there, the relative magnitude of the fluctuations is [r.m.s.(P2)]/(-P2) ~ 1/fi 1 - At 1/g <C 
P2I '(P2) -1<1 the distribution function is of the exponential form, 

P( P 2) . exp {-f I (<!-!)}. (4,02) 

Note that for negative deviations ^2/(^2) — 1 with ^2/(^2) _ 1| 3> 1/g the distribution 
function decays much faster ||137|| , so that the distribution is strongly asymmetric. As was 
mentioned in Sec. |3.3.2| , the "body" of the distribution V{Pq) is described properly (in the 
leading order in 1/g) by the Liouville theory. We will see below that this is also true for the 
asymptotic "tail" ofP(P 2 ). 

Our consideration of asymptotics of the IPR distribution is based on unpublished results 
p5| (partially announced in |3!J]). We derive first a relation between V(P2) and the distri- 
bution of level velocities V v (v) f25|. To this end, we consider a Hamiltonian 7i + aW, where 
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W is a random perturbation. Specifically, the matrix elements W rir2 are supposed to be 
independent Gaussian-distributed random variables with the mean value equal to zero and 
the variance 

(MC«%) = W o(ln - r a |)*(n - ri)*(r a - r' 2 ). 

We will assume that Wo(r) is a short-ranged function with some characteristic scale £. The 
level velocity t>„ corresponding to an energy level E n is defined as v n = dE n (a)/da (where 
E n (a) is the level of the perturbed Hamiltonian Ji + aW) and can be found within the 
conventional perturbation theory as 



d d r I d d r'Wrr'ip* n (r)ipn( r ') ■ 



Using ( |4.103| ), we find 



[v 2 n ) = w P 2 ; w = cj AWo(r) 



(4.103) 



(4.104) 



where c = 1 (1/2) if k f ( < 1 (resp. h f ( > 1). 

This consideration can be extended to higher moments of the level velocity (vn q ) as well. 
This leads to the following relation between the two distributions: 



VJv) 



(5(v - v n )) 
dP 2 



o [2nw P 2 



1/2 



exp 



2w P 2 



(4.105) 



On the other hand, the level velocity distribution can be expressed through the cx-model 
correlation function in the following way [|1 56|| . According to the definition, 



VJv) 



lim ^ (^2 6 ( E ~ E n(0))S(E + av- E n (a))J 



(4.106) 



so that V v {v) is determined by the a — > limit of the parametric level correlation function. 
The latter can be represented [as a generalization of Eq. ( 2.10|) 1 in terms of a correlation 
function of the cx-model ||157| , |158| , |159|| . This yields 



V v {v) = lim -ztj^-^j- 
v^o 8V 2 vVv 



DQ ( / d d rStrQ u k) ( / d d rStrQ 22 k ) e~ Sv[Q] ; 



S V [Q] 



cfrStr 



^ D <V Q r-^QA + ^^QAQA 



4 



4v< 



(4.107) 



Combining Eqs. ( [4.1 05| ) and ( [4.107] ), we find the expression for the IPR distribution function 
in terms of the cx-model, 



Vi{P 2 



-1 



1 



8V 2 2zv / 2^R 3/2 



lim 7]uV 



c+ioo 



du 



x 



J DQ (^J d d rStxQ u kj (^J d d r$txQ 22 kj e" 5u[Q] ; (4.108) 
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where S U [Q] = — u/2P 2 + S v [Q]\ v 2/ m=u . The saddle point configuration is again parametrized 
by the bosonic "angle" 9(r) only; the action on such a configuration is 



S u [9) 



a 



2P, 



+ / d d r 



4u 



(4.109) 



The corresponding saddle-point equations can be readily obtained by varying S u [9] with 
respect to 6*(r) and u: 



. D ^0 _ % JL e e + 



2u 6 







{7TU7])2 ' d d re 2e + 1 



. 



(4.110) 



4m 2 J ~ ' " ' 2P 2 

Making a shift = + In^iu/irvr]) and dropping the tilde, we can reduce them to the form 

V 2 9 - 7 (e e - e 20 ) = ; (4.111) 
P> = 2 / / d d re 2d , (4.112) 



where 7 = uj2i:vD. Integration of Eq. ( |4. 1 1 1| ) with the Neumann boundary condition yields 



re — e 



29 \ 



0, so that Eq. ( [4.1 12|) can be rewritten in the following form (invariant 



with respect to a shift of the variable 



Pi 



r e 



(J d d re e Y 



(4.113) 



The meaning of Eq. fl4.113Q is completely transparent if we recall that |V' 2 ( r )l smoo th 00 ■ 
The factor 2 comes from the GUE-like short-scale fluctuations. 

Taking into account Eqs. ( |4.111| ), ([4 . 1 1 2| ) , we find the action (|4.109| ) on the saddle-point 
configuration to be equal to 



Su[9] = ^- Jd d r(W) 2 . 



(4.114) 



The problem that we are solving is easily seen to be equivalent to searching for the minimum 



of ( |4.1 14 ) under the conditions J e 6 d d r = 1, J e 26 d d r = P 2 (or, equivalently, under the 



condition (|4.113| ) invariant with respect to normalization of e e ). This is nothing else but the 
optimum fluctuation problem for V(P 2 ) within the Liouville theory ( |3.93j ). Therefore, the 
Liouville theory (|3.93|) describes properly the asymptotics of the IPR distribution. 



Quasi-ID geometry. Equation 



has the following general solution 

„-0 7 



nil 8 1 

9 — 7e + 7e 



26 



Ci 



1 + v 7 ! - C1/7 sin( ^/C[x - C 2 ) 

The constants C\ and C 2 should be found from the boundary conditions 

9'{-L/2) = 9\L/2) = , 



(4.115) 



(4.116) 
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yielding 

cos ± C^j = . 

The solution providing the minimum to the action corresponds to G% = 0, \J~C~\ = tc/L and 
gives the wave function intensity 

^Mlsmooth = J-e^ = 7 1/2 ^ 2 1 + v^^VZ^sin(vrx/L) (4 ' H7) 
(with A being the sample transverse cross-section) and the IPR value 

P 2 = ^r- = P 2 GUE -^— -. (4.118) 
yi7r 7T 



Calculating the action ( |4.114j ) we find the following asymptotic behavior of the IPR distri- 
bution function 

_ m) = s. = ^ (ft - A) = £ (j*. - l) . (4,19) 

Eq. Q4.119Q is valid for l/g < P 2 /A GUE ~ 1 < V^- Therefore, in the quasi-lD case the 
exponential behavior (|4.102|) is not restricted to the region of small deviations from the 
average value; there is no change of the behavior of V(P 2 ) at P2/P 2 GUE — 1 ~ 1 (we will find 
such a change below in the 2D geometry case). 

The far asymptotics at P 2 S> P 2 GUE * s determined by the states, which have an effective 
localization length £ e f much smaller than the sample length L. The shape of these states is 
according to ( [4. 117p 

S^(r)\ smooth = 2 - [x + L ^ +ei ; (4-120) 

Let us note that the density of such a state is concentrated near the edge (x = —L/2) of 
the sample (of course, there exists an equivalent solution located near the opposite edge, 
x = L/2). Since the shape of a state with £ e f <C L and the corresponding action S u ~ 
7i 3 i/DA 2 P 2 /4: are only weakly affected by the sample length L, we may expect such states to 
determine the asymptotic behavior of Vi(P 2 ) also in the localized regime, I>(. However, 
if we compare the action found above with the asymptotics of ViiP-i) at P 2 ^> 1/C.A in an 
infinitely long sample [first line of Eq. fl3.49p 1, we find that the absolute value of the exponent 
in the latter case is larger by factor of 4. The explanation is as follows. The states determining 
the large- P 2 asymptotics in the limit L —>■ oo can be obtained from Eq. (|4.116|) by taking a 
solution which has a maximum in the middle of the sample (C 2 = tt/2, \fC\ = 2n/L); the 
result is 

^Mlsmooth = \JfQ f : = ^ ■ (4421) 



The only difference between (|4.120| ) and ( 4.121|) is that ALS is now located in the bulk of the 



sample. This leads to an extra factor 4 in the action (see similar discussion in Sec. [4.2.3| for 
the case of the statistics of eigenfunction amplitudes). However, in the limit of a long sample, 
1 the contribution of the states located near the boundary is additionally suppressed 
by a factor ~ 1/L as compared to that of the bulk ALS (which may be located everywhere 
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in the sample). Therefore, if one fixes P 2 and considers 7 ? j(P 2 ) in the limit L — > 00, only 
the contribution of the bulk ALS survives, despite the fact that it has the exponent 4 times 
larger than that of the ALS located near the boundary. In other words, the contribution of 
the states located near the boundary to the first line of Eq. (|3.49| ) (large- P 2 asymptotics of 
P 7 (P 2 ) at X = > !) is oc ^e-^/ 4 . 

On the other hand, if a sample with periodic boundary conditions in the longitudinal 
direction (a ring) is considered, only the bulk solution (C 2 — tt/2, \fC[ = 2n/L) will survive. 
Consequently, the asymptotic form of mP(P 2 ) will be different from ( 4.11 9|) by an extra 
factor of 4. 



2D geometry. Now we calculate the far asymptotics of V(P2) at P2 3> P 2 GUE . We assume 
first the periodic boundary conditions (so that the sample has no boundary) and search for 
a rotationally-invariant solution, which should satisfy the differential equation 

6" + -& - 1 (e d -e 2e ) = . (4.122) 
r 

(In fact, for the hard-wall boundary conditions the asymptotics is determined by the states 
located near the boundary; however, such states can be obtained from the symmetric solution 
by putting the center at the boundary and restricting the solution to the interior of the 
sample; see below.) From our experience in the quasi- ID case, we expect the solution to 
have a form of a bump concentrated in a region r < l p and decreasing with r outside this 
region. For r ^> l p the term e 2d is thus expected to become irrelevant, so that the equation 
takes the form fl4.75|) . Its solution is given by ([4.79|) and can be approximated at l p <C r <C L 



as e 9 ^ ~ A/r a with some coefficient A and exponent a. In order to have large P 2 , we require 
a > 1. Using the condition ( |4.113| ), we find scale l p to be given by 



P 9 1/2 , a > 2 



lp ~ ^ (£^2ap 2 )-l/[2(a-l)] ; Q < 2 [ ( 4123 ) 



up to a numerical coefficient of order unity. Therefore, the action (|4.114j) is equal to 

S u = a 2 n 2 isD \n(L/l p ) = T{u)k 2 vD ln(L 2 P 2 ) , (4.124) 



where 



, a 2 12 , a > 2 ; . 



2(o-l) ' 

Thus, the minimum of S u corresponds to a = 2, yielding 



P, 



-P*g/2 



P(P 2 ) ~ ( ) (4.126) 



The upper border of validity of d4T26|) is P 2 ~ P 2 RMT (L//) 2 ~ l/l 



]2 



For hard wall boundary conditions, the asymptotic behavior of P(P 2 ) will be, however, 
determined by configurations with a maximum located near the sample boundary. Assuming 
that the sample has a smooth boundary with the single characteristic scale L (for example, 
it is of the circular form), we can get such a state from the rotationally invariant bulk state 
by putting its center on the boundary and removing that half of the state which is outside 
of the sample. Such a truncated state will have the twice larger IPR and twice smaller 
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action compared to its parent bulk state. Consequently, the asymptotics of the distribution 
function for a sample with a boundary will be different from ( |4.126| ) by an extra factor 1/2 
in the exponent. 

Using the Liouville theory description ( |3.93|) , one can generalize the above consideration 
to the distribution V(P q ) of higher IPR's P q ( |3.90| ) with q > 2 ||160|| . We will assume that q 
is not too large, q 2 < 2Png, so that the average value (P q ) is at the same time the typical 
value of P q (see Sec. p.3.2| ). Then in the region q 2 /f3ng < P q /(P q ) — 1 < 1 the distribution 
has the exponential form (|3.92 ). At larger P q the optimal configuration is again of the form 
| (r) | smooth = = A/r a for r > Z p ; minimizing the action, we find a = g/4 and the 



smooth — 6 

distribution function 



Pn 



P, \{P q ) 



(4.127) 



This is valid for P q < P RMT (L//) 2 ; for still larger P q the corresponding optimal fluctuation 
would violate the condition of the applicability of the diffusion approximation 9' < 1/7. 
Incorporating this restriction (cf. similar situation for the distribution of relaxation times, 
Sec. |X2|) leads to 



a 



ln(P g /P RMT ) 
ln(L/Z) 



- 2 



and to the log-normal far asymptotics 

1 



V(P 



qj ~ 75" eX P 

r 1 



Aq 1 



\n(L/l) 



(4.128) 



(4.129) 



for (L/l) 2 < P 9 /P RMT < (L/l) 2q 



-2 



4.5 3D systems 

As we will see below, in the 3D case the states determining the asymptotics of the distribution 
functions have just a local short-scale spike on top of a homogeneous background. In this 
sense, no ALS is formed, in contrast to the quasi-lD and 2D situations. As a closely related 
feature, we find that the results in 3D are strongly dependent on microscopic details of the 
random potential. We start from discussing what the cr-model calculation gives when applied 
to the 3D geometry. Then we compare this with the results of the direct optimal fluctuation 
method [|161 |. 

Let us consider e.g. the long-time relaxation problem. The solution of the cr-model 



saddle-point equation has the form |30 



6{r) 



l*<r<R, 



where R is the system size (which is, in fact, irrelevant here) and 



L ~ I In 



r(k F iy 



C ~ In 



r(k F iy 



(4.130) 



(4.131) 



The reason for introduction of the short-scale cut-off length is the same as in 2D (see 
Sections |4.1.2| , [4.2.2|) : the gradient 9' should not exceed ~ l/l. Calculating the action, we 
find 

-\nP s (t) ~ (M) 2 hi 3 ' 



r(k F iy 



(4.132) 
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with an uncertainty in numerical prefactor. This uncertainty originates from the fact that 
the action is dominated by the ultraviolet (short-distance) region r ~ £*. Similar result is 
obtained for the eigenfunction amplitude statistics |33|. |3C 



Vv 

\nV(u) ~(k F l) 2 In 3 — — (4.133) 

(Mr 



and for the LDOS distribution function f3~4}] 



-\nV{p) ~ (M) 2 ln 3 - . (4.134) 

As we see, all these distribution functions are found to have the exponential-log-cube asymp- 
totic form. The corresponding eigenstates have the shape 

I V^M (smooth ^ ^ eX P j^" 1112 ^ 

where Ci ~ 1, Z t — t/[r(kpl) 2 ], Z u = Vuj (kpl) 2 , and Z p = p/u, i.e. they consist of the 
usual homogeneous background \ip 2 (r)\ smooth = 1/V supplemented by a short-scale bump 
containing only a minor part of the eigenfunction normalization. 

Muzykantskii and Khmelnitskii |3]J applied the ballistic a-model approach to the problem 
of the long-time relaxation. They reproduced the result (|4.132|) and found the numerical 
prefactor there to be equal to 7r/9\/3. On the other hand, Smolyarenko and Altshuler ||161 



employed the direct optimal fluctuation method to study the asymptotics of the distribution 
function V(u). This approach is similar to that used many years ago |162. 163| to calculate 



the tails of the density of states of a disordered system. The authors of ||161|| did not use the 
ex-model description, but rather searched directly for the optimal realization of the random 
potential (hence their term "direct optimal fluctuation method"). While having confirmed 
the exponential-log-cube form of the asymptotics in the 3D case, they found a prefactor 
smaller by ~ kpl compared to the a-model result, i.e. they obtained 

V(u) ~ k F l\n 3 (uV) . (4.135) 

The physical reason for this difference lies in the ballistic effects which have been already 
discussed in Sec. |3.3.3| in connection with perturbative corrections to the eigenfunction am- 



plitude distribution. It was shown there that in the 3D case and for the white-noise random 
potential the parameter k governing these corrections is dominated by a non-universal (de- 
pending on the type of the disorder) ballistic contribution yielding k ~ 1/kpl, while the 
diffusive contribution is ~ l/(kpl) 2 . This is in direct correspondence with Eqs. ( |4. 1 33|) , 
( |4.135| ), which show precisely the same difference. Eq. ( f4.135|) is again non- universal; its 



derivation by Smolyarenko and Altshuler ||161|| relies on the white-noise disorder assump 



tion. The corresponding optimal configuration of the potential found in ||161|| is nothing 
else but a potential wall surrounding the observation point with the height several times 
larger than Ep and the thickness ~ Xpln(uV). Such configurations are not included in the 
cx-model consideration, which assumes that the absolute value of the particle velocity does 
not change appreciably in space. If one would consider a smooth random potential, whose 
magnitude is limited from above by some value t/ max Ep, such configurations would not 
be allowed. Whether in this case the cx-model result would hold remains to be seen. 
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4.6 Discussion 



In Section £|, we have studied the asymptotic behavior of various distribution functions 
characterizing the eigenfunction statistics in a disordered sample. For this purpose, we used 
two methods of treatment of the cr-model: exact solution (in the quasi-lD case) and the 
saddle-point method. Physically, the saddle-point solution describes the relevant optimum 
fluctuations of the wave function envelope; probability of formation of such a fluctuation is 
found to be governed by the Liouville theory (|3.93| ). In the quasi-lD case, the results of the 
saddle-point method are in agreement with those of the exact solution of the cr-model. 

The 2D geometry is of special interest, since the eigenfunctions show the features of 
criticality. In this full agreement between the saddle-point calculation and the 

renormalization-group (RG) treatment of Altshuler, Kravtsov and Lerner |28|] was found 
for all the distributions, where such a comparison was possible, namely for P(t^), V(p) and 
V{y). This agreement is highly non-trivial, for the following reason. The RG treatment 
is based on a resummation of the perturbation theory expansion and can be equally well 
performed within the replica (bosonic or fermionic) or supersymmetric formalism. At the 
same time, the present approach based on the supersymmetric formalism relies heavily on 
the topology of the saddle-point manifold combining non-compact (Ai) and compact (A2) 
degrees of freedom. The asymptotic behavior of the distribution functions considered is 
determined by the region Ai ^> 1 which is very far from the "perturbative" region of the 



manifold Q ~ A (i.e. Ai,A2 — !)• It is well known ||164|| that for the problem of energy 



level correlation, the replica approach fails, since it does not reflect properly the topology 



of the cr-model manifold. The success of the RG treatment of [gg] seems to be determined 
by the fact that for the present problem (in contrast to that of level correlation) only the 
non-compact sector of the supersymmetric cr-model is essential, with compact one playing 
an auxiliary role. Let us note that the same situation appears in the vicinity of the Anderson 
transition |9Cj, [19], where the function Y(\i,\2) acquires a role of the order parameter 
function and depends on the non-compact variable Ai only. The above agreement found for 
the "tails" of the distributions in the metallic region provides therefore support to the results 
concerning the Anderson transition obtained with making use of the renormalization group 
approach and 2 + e expansion. 

We have found that the spatial structure of ALS relevant to the asymptotic behavior 
of different distributions may be different. This is because an ALS constitutes an optimal 
fluctuation for one of the above quantities, and the form of this fluctuation depends on the 
specific characteristic, which is to be optimized. Finally, we have discussed interrelations 
between asymptotics of various distribution functions. In the quasi-lD and 2D cases, we thus 
presented a comprehensive picture which explains all the studied asymptotics as governed 
by exponentially rare events of formation of ALS. 



5 Statistics of energy levels and eigenfunctions at the 
Anderson transition 

In d > 2 dimensions a disordered system undergoes, with increasing strength of the disorder, 



a transition from the phase of extended states to that of localized states (see e.g. [39| for 
review). This transition changes drastically the statistics of energy levels and eigenfunc- 
tions. In particular, at the mobility edge these statistics acquire distinct features reflecting 



criticality of the theory. This is the subject of the present section. Sec. fO is devoted to 
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the level statistics and Sec. [x2] to the eigenfunctions correlations at the mobility edge. In 
Sec. |0| we study the level and eigenfunction statistics in a quasi-one-dimensional model 
with long-range (power-law) hopping which undergoes the Anderson transition and shows at 
criticality all the features characteristic for a conventional metal-insulator transition point 
in d > 2. 



5.1 Level statistics. Level number variance 

The problem of the energy level statistics at the mobility edge was addressed for the first time 
by Altshuler et al, Ref. |0], who considered the variance (5N 2 (E)) = (N 2 (E)) - (N(E)) 2 of 
the number of levels within a band of a width E. This quantity is related to the two-level 
correlator ( |2.1p via 

r (N(E)) 

(6N(E) 2 } = / ((N(E)) - \s\)RV(s)ds , (5.1) 



(N(E)) 



or, equivalently, 



■(5N(E) 2 ) = / R^{s)ds , (5.2) 



d(N(E)) v ' 1 J_ {N(E)) 

where (N(E)) — E/A and R^ c '(s) = R(s) — 1 is the connected part of the two-level correlation 
function. In RMT, the 1/s 2 behavior of R^ c '(s) leads to the logarithmic behavior of the 
variance {SN 2 ) WD ~ (2/n 2 f3) ln(iV) for (N) ^> 1. In the opposite situation characteristic for 
the phase of localized states, when all energy levels are completely uncorrelated (known as 
the Poisson statistics), one gets (SN 2 )p = (N). Supported by their numerical simulations, 
Altshuler et al. 0] put forward a conjecture that at the critical point 

(57V 2 ) ~ X {N) , (5.3) 

where < x < 1 is a numerical coefficient (which is conventionally called now "spectral 
compressibility"). More recently, Shklovskii et al. introduced the concept of new uni- 
versal statistics at the mobility edge (see also Ref. ||165|| ). In Ref. the correlator R(s) 
at the mobility edge was studied by means of perturbation theory combined with scaling 
assumptions about a form of the diffusion propagator. It was found that for s ^> 1, 

R {c) (s) oc^T 2+7 (5.4) 

where 7 < 1 is certain critical index. The consideration of Ref. led to the conclusion that 



7 = 1 — 1/vd (where v is the critical exponent of the localization length), which was however 
questioned later ||166|| in view of the oversimplified treatment of the diffusion propagator at 



the transition point in [fj2 |. At any rate, the behavior ( |5.4| ) with some 7 is what one expects 
to hold at the mobility edge; the condition 7 < 1 follows from the requirement of convergence 
of / R(s)ds. 

Using Eqs. (|5.1|) , (|5.4j ) and the sum rule 

R(s)ds = (5.5) 



(implied by the conservation of the number of energy levels), the authors of [42| concluded 
that 

(5N 2 ) oc iV 7 (5.6) 
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with 7 < 1, in contradiction with Ref. [EO]. This conclusion was critically reexamined in 



Refs. [43], §4|, where it was shown that the asymptotic behavior (|5.4j) of the correlator R(s 



at s ^> 1 does not imply the absence of the linear term (|5.3|). The flaw in reasoning of 



Ref. [42j was in the assumption that the universal part of the correlator R(s) (which is the 
one surviving in the limit E/A = (N) = const, L — > oo) satisfies the sum rule ( |5.5| ). It turns 
out, however, that the sum rule is fulfilled only if all contributions are taken into account, 
including the non-universal contribution of the "ballistic" region uj ~ 1/r, where r is the 



elastic mean free time. To demonstrate this, we estimate below (following Ref. [43|) the 
contributions to the sum rule from all regions of the variable s. 



In fact, for the conventional model of a particle in a random potential the definition (|2.1| ) 
of the correlator and the sum rule relation ( |5.4|) should be modified, when the vicinity of the 
critical point is considered. The reason is that the Anderson transition point corresponds 
to a strong disorder regime EpT ~ 1, so that the condition u ~ 1/r implies to ~ Ep. 
On the other hand, the density of states v{Ep) can be considered as a constant only for 
small variations of energy uj <C Ep. This means that variation of v{Ep) should be taken 
into account in fl2.1|), (|5.5|). Besides, the condition EpT ~ 1 leads to a breakdown of the 
perturbation theory, that complicates the analysis of the "ballistic" region contribution. To 
get rid of these problems, we consider a different microscopic model which has exactly the 
same universal part of R(s), but whose density of states does not change within the range 
of uj ~ 1/r. This is so-called n-orbital Wegner model [|167|| , which can be considered as a 



system of metallic granules forming a d-dimensional lattice, each granule being coupled to 
its nearest neighbors. In the limit n ^> 1 this model can be mapped onto a supersymmetric 
cr-model defined on a lattice. The action of this cr-model reads |5U| (we consider the unitary 
symmetry for definiteness) 

s{Q} = J -J2 StI QiQj + e £ StlQiA ' (5 - 7) 

(v) * 

where the supermatrices Qi are defined on sites i of a (^-dimensional lattice with a lattice 
spacing a. Summation in the first term of Eq.( |5.7| ) goes over the pairs (ij) of nearest 
neighbors. The parameters T and e are related to the classical diffusion constant D, the 
density of states z/, and the frequency to as follows: 

T = ixvDa d ~ 2 , e = -iujixva d /2 . (5.8) 

The two-level correlator can be expressed through a correlation function of this a-model 
via the discretized version of Eq. fl2.10|) : 



R{s) = Re J Y[ DQj 



4V 



StrQikA 



- 1 } exp(-S{Q}) , (5.9) 



where V is the system volume. It is not difficult to prove explicitly P3| that the two-level 
correlation function of this lattice a-model, Eq. (|5.9|) , satisfies exactly the sum rule (|5.5|) . 
As will become clear below, the continuum version of the cr-model does not possess this 
property: there is a deficiency of the sum rule for it which is related with the contribution 
of the range of lj close to the ultraviolet cut-off. 

Let us stress that in the region u -C D/a 2 = E c (a) the correlator R(s) is universal, i.e. 
does not depend on microscopical details of the model. The region uj ~ D/a 2 plays a role 
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analogous to that of the ballistic region, u ~ 1/r, in the case of the usual model of particle 
in a random potential. Despite the non-universality of the correlation function R nu (s) in this 
"ballistic" domain, the corresponding integral contribution I nu to the sum rule is universal, 
because it determines, according to Eq. (|5.5| ), the sum rule deficiency for the universal part 
Ru(s): 

4 = / R$(s)ds ; I nu = J R { n c l(s)ds (5.10) 

As will be seen below, when the system is close to the Anderson transition the two regions 
of the variable s dominating the integrals I u and I nu , respectively, are separated by a para- 
metrically broad range of s giving a negligible contribution to the sum rule. For energy band 
width E lying in this range, we get from Eqs. fl5.2p , ( |5.10| ) 

(6N(E) 2 } ~ I U (N(E)} = -I nu (N(E)} , (5.11) 



'nil' 



i.e. just the linear term ( |5.3| ) with x = l u = —I n 

We turn now to the analysis of the correlator R(s) and of the sum rule in various situation. 
For completeness, we start from the case of a good metal, then we consider the critical point 
and the critical region cases. 

(1) Good metal. Here the following three regions with different behavior of R(s) can be 
found: 

(A) "Wigner-Dyson" (WD) region: u <C E c . The correlator R(s) in this region was 
studied in Sec. pT2| , see Eqs. ( |2.27| ), ( |2.30| ). The corresponding contribution to the sum rule 
can be estimated as (we omit numerical factors of order unity) 



[ 9 R(s) oc +l/g . (5.12) 
Jo 



(B) "Altshuler-Shklovskii" (AS) region: E c < u < E c (a). Here E c (a) = D/a 2 is the 
Thouless energy at the scale of lattice spacing a, which plays a role of the ultraviolet cut-off 
for the diffusion theory. The level correlation function is given by Eq. ( |2.33| ), yielding 



IB 



D/a 2 A / n \ |-1 

R(s)ds oc g~ d ' 2 — oc +l/g(a) , (5.13) 



where g(a) = g(a/L) d 2 is the conductance at the scale a. Note that for the case of a particle 
in a random potential, the following substitutions should be done: a — > I; E c (a) —> E C {1) = 
1/r; l/g(a) -> l/g(l) oc (er)^ 1 . 

(C) "Ballistic" region: uo > E c (a). To find the correlator R(s) in this range, we can 
neglect in the leading approximation the first term in Eq. ( |5.7p , that gives 

L\ d 1 



R^(s)ec- - -, (5.14) 

\ CL / S 

poo 

R {c \s)ds oc -l/g{a) . (5.15) 

J D/a 2 A 



The contribution ( |5.12| ) to the sum rule is dominated by the region u ~ A, whereas the 
contributions fl5.13|) , ( |5.15|) are dominated by a vicinity of the ultraviolet cut-off u ~ D/a 2 . 
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In d — 3 g ^> g(a), and these non- universal contributions Ib, Ic are much larger (by absolute 
value) than I a, but they should cancel each other according to the sum rule. The situation 
is more interesting in d = 2, where the contribution ( j2.33| ), (|5.13| ) of the AS region is absent 
in view of special analytical properties of the diffusion propagator ||168| , |169| , j44fl . In this case 



g(a) = g, so that the contributions Ia and Ic are of the same order, in agreement with the 
sum rule which prescribes their sum to be zero. We can identify then Ia with I u and Ic with 
I nu in Eq. ( |5.10| ), so that I u = —I nu oc 1/g. According to Eq. ( |5.11| ), this leads to a linear 
behavior of the variance, 

(5N 2 ) oc (l/g)(N) ; E c <^E^D/a 2 , (5.16) 

but the corresponding coefficient is of order 1/g <C 1. The numerical coefficient in (|5.16|) 
can be calculated using the explicit form of R^ c \s), yielding 

X = 4 - 7T7^ , 9 > 1 ■ (5-17) 

2f3ng 



With an increase in disorder strength, the coupling constant T in Eq. (|5.7| ) decreases. 
When it approaches a critical value T c , corresponding to the metal-insulator transition, the 
correlation length £ becomes large: £ a. Depending on the relation between £ and the 
system size L, the following two situations are to be distinguished. 

(2) Metal in the critical region: a <C £ <C L. In this case g(a) ~ g*, but g = g(L) g*, 
where g* is the critical value of the conduct anceQ. We find then as much as four different 
regions of s for the correlator R(s): 

(A) WD region: uj <C E c . The correlator R(s) and the contribution to the sum rule for 
this region are given by the same Eqs. ( |2.27| ), ( |2.30| ) and ( |5.12| ), as for a good metal. 



(B) AS region: £ c C w < £ C (C) = 9*^, where E c (£) and Ag are Thouless energy and 
level spacing for a piece of the sample with linear size £. The correlator R( c \s) is given by 
Eq. ( |2.33|) , and we find 



9 



9*A e /A /A \5-l 

R^\s)ds oc g~ d ' 2 g, ^ oc +l/g, , (5.18) 



where we have used the relations A^/A = (L/£) d and g/g* oc (L/£) d ~ 2 . 

(C) "Kravtsov-Lerner-Altshuler-Aronov" (KLAA) region: E c (£) = g*A^ <C uj <C g*A a = 
E c {a). The correlator R^(s) in this range was studied in Ref. [ |170|| ; the result reads: 

R^(s) ~ -g? (A^/A) 1 ^ s- 2+ ^ (5.19) 

and, consequently, 

Ic 



rg*A a /A 

/ R®(s)ds<x-l/g, (5.20) 

(D) "Ballistic" region: uj > E c (a). Here the equations (|5.14j) , (|5.15| ) hold, yielding 



I D oc -1/g, (5.21) 

The contribution I a is dominated by uj ~ A, the contributions Ib, Ic by w ~ A^, and 
finally, the contribution Ip comes from the "non-universal" region uj ~ EJa). Therefore, 



7 In 3D the value of g* is, of course, of order unity. However, if one consider formally d — 2 + e with e <C 1, 
then g» ~ 1/e is parametrically large. It is thus instructive to keep g„ as a parameter in all estimates. 
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this last contribution determines I nu in Eq. ( |5.10| ), leading according to Eq. (|5.11| ) to the 
linear behavior of the variance 

(5N 2 ) oc (l/g*)(N) ; E c (0 < E < E c (a) (5.22) 

with a coefficient ~ (which is of order unity in 3D). 

(3) Finally, a system is at the critical point, when £ ^> L. In this case, g(£) = g(a) = g*- 
The AS region disappears, and we have, in full analogy with the previous estimates: 

(A) WD region: uj < E c = g»A. 

I A oc +l/g* (5.23) 

(B) KLAA region: g*A < to < A a g*. 

R(s) oc -^ 7 s- 2+7 ; I B oc -1/g, (5.24) 

(C) "Ballistic" region: u > A a g*. 

I c oc -1/g, (5.25) 

Again the same conclusion, as for the critical region case, can be drawn: the "ballistic" 
contribution ( |5.25|) can be identified as I nu in Eq. ( |5.10|) , yielding: 

(5N 2 ) ~ X (N) , £ c < £ < E c (a) , (5.26) 

with x ~ 1/fl 1 *- 

Recently, a relation between the spectral compressibility x an d the multifractal dimension 
D 2 (for discussion of multifractality and corresponding bibliographic references, see Sec. |3.3. 2 
and |5.2| ) was proposed [f|^, p.66|| : 

X = ■ (5.27) 

The central idea of the derivation ^| is to consider the motion of energy levels, when the 
system is subject to a random perturbation. This allows to link the spectral statistics with 
the wave function correlations. In 2D, one can check, by comparing Eqs. ( |3.63|) and ( |5.27|) , 
that Eq. ([5.271) is indeed satisfied in leading order in 1/g. The general derivation of Eq. (|5.27|) 
is, however, based on certain approximate decoupling of a higher order correlation function 



171| , |172| , |49| ) so that it is not completely clear whether this is indeed an exact formula as 
argued in |^| or only an approximation valid for g* ^> 1. 

The linear behavior (|5.3| ) of the level number variance at the Anderson transition has been 
confirmed by now in numerical simulations by several groups E3, EBJ, E7I, 03. The spectral 



compressibility x is an important universal parameter characterizing the critical point of 
Anderson transition. Its universality is of the same sort as that of the critical indices, i.e. 
X depends only on the spatial dimensionality and on the symmetry (universality) class. Let 
us note that the whole level correlation function R^(s) is not as universal, since it depends 
also on the shape of the sample and on the boundary conditions [f|l| |173| , |174| , |175|| . This 



can be expected already from the perturbative 1/g 2 correction, Eqs. (|2.27|) , fl2.30|) , where 



the coefficient does depend on the shape and on the boundary conditions. Also, the 
shape dependence becomes evident if one considers the limit of an elongated sample with a 
length considerably larger than the transverse sample size. Indeed, let us consider (in 3d) 



a rectangular sample with L y = L z = L x /a, where a is a numerical factor ||173 |. If we fix 



a and consider the limit L — > oo, the level correlation function has a limiting form, which, 
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depends, however on a. In particular, at a ^> 1 the sample is of quasi- ID geometry with a 
ratio of the sample length to the localization length ~ a, so that the level statistics will be 
close to Poissonnian. 

Asymptotic behavior of the nearest-neighbor level spacing distribution function P(s) for 
s ^> 1 at the mobility edge has been also a controversial issue. We remind that for the 
Poissonnian statistics P(s) = e~\ while in the RMT P{s) ~ e - consts2 . While Refs. [0, |T 



conjectured that P(s) ~ e consts ; the authors of Ref. [|176|| found P(s) ~ e consts with the 



same index 7 as in Eq. (5^). Recent numerical studies [177, 178| | support the former result, 

P(s) ~ g-consts_ 



5.2 Strong correlations of eigenf unctions near the Anderson tran- 
sition 

In this section, we discuss correlations between the amplitudes of different, but close in 
energy, eigenstates at and near the critical point of the Anderson transition. Let us recall that 
in the metallic phase far from the transition point a typical eigenfunction covers essentially 
uniformly the sample volume. This is reflected in the inverse participation ratio P2, as well as 
in the higher moments P q = J d d r \ip 2q (v)\, which differ only weakly from their RMT values, 
see Eqs. fl3.51| ), (|3.52|) . When the system approaches the point of the Anderson transition E c , 



these extended eigenfunctions become less and less homogeneous in space, showing regions 
with larger and smaller amplitudes and eventually forming a multifractal structure in the 
vicinity of E c , see Sec. |3.3.2| . This multifractal behavior is characterized by the following 
behavior of the moments P q at the critical point 

P q oc L- D ^ q - l) , L < £ , (5.28) 

as well as in the conducting phase in the vicinity of the critical point: 

P q oc ^-o q ){q~i) L -d{ q -i) ; L>£ . (5.29) 

Here £ oc \E — E c \~ u is the correlation length and D q is the set of multifractality exponents, 
d > D 2 > D 3 > . . .. As Eq. (|5.29|) indicates, the eigenfunctions become more and more 
sparse when the system approaches the critical point of the Anderson transition from the 
metallic phase (i.e. when £ increases). Just at the mobility edge the scaling of IPR with the 
system size L becomes different, see Eq. ( |5.28| ), so that the eigenfunction effectively occupies 



a vanishing fraction of the system volume. At last, in the insulating phase any eigenstate 
is localized in a domain of finite extension £, and IPR remains finite in the limit of infinite 
system size L — > 00. 

This transparent picture serves as a basis for qualitative understanding of spectral prop- 
erties of disordered conductors. Indeed, as long as eigenstates are well extended and cover 
the whole sample, they overlap substantially, and corresponding energy levels repel each 
other almost in the same way as in RMT. As a result, the Wigner-Dyson (WD) statistics 
describes well energy levels in a good metal, see Section 0. In contrast, in the insulating 
phase different eigenfunctions corresponding to levels close in energy are localized far apart 
from one another and their overlap is negligible. This is the reason for absence of correlations 
of energy levels in this regime, i.e. the Poisson statistics. 

A naive extrapolation of this argument to the vicinity of the transition point would lead to 
a wrong conclusion. Indeed, one might expect that sparse (multifractal in the critical point) 
eigenstates fail to overlap, that would result in essential weakening of level correlations close 
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to the mobility edge and vanishing level repulsion at E 
H79i 



E r . However, numerical simulations 



show 



180, 181 



that even at the mobility edge levels repel each other strongly 
(although the level statistics is different from RMT). As we are going to explain now, this 
apparent contradiction is resolved in the following way: The critical eigenstates for nearby 
levels are so strongly correlated that they overlap well in spite of their sparse structure. 

To demonstrate the strong correlations, we will consider the relation between the overlap 
function cr(r, r, E, uj), Eq. fl3.68| ), and the "self-overlap" a(r, r, E), Eq. fl3.67| ) (the latter de- 
termining the IPR P2 according to P2 = J dra(r, r, E). As was shownQ in ||182|| , everywhere 
in the metallic phase (including the vicinity of the critical point) 



cr(r, r, E, u>) 



(3 



(3 + 2 



a(r, r, E) 



(5.30) 



for 10 <C A^, where Ag ~ l/^C d is the level spacing in the correlation volume. Eq. (|5.30| ) 
implies the following structure of eigenf unctions within an energy interval to < A^. Each 
eigenstate can be represented as a product ^(r) = ipi{v)^Eiy)- Here the function $£(r) is an 
eigenfunction envelope of "humps and dips" . It is the same for all eigenstates around energy 
E, reflects the underlying gross (multifractal) spatial structure, and governs the divergence 
of the inverse participation ratios P q at the critical point. In contrast, ipi(r) shows RMT- 
like fluctuations on the scale of the wave length. It fills the envelope function $^( r ) hi an 
individual way for each eigenfunction, but is not critical, i.e. is not sensitive to the vicinity of 
the Anderson transition. These Gaussian fluctuations are responsible for the factor (3/ {(3 + 2) 
(which is the same as in the corresponding Gaussian Ensemble) in Eq. Q5.30 ). 

As was already mentioned, this picture is valid in the energy window u 
the energy E; the number of levels in this window being large as A^/A ~ (L/£) d 3> 1 
in the limit L ^> £. These states form a kind of Gaussian Ensemble on a spatially non- 
uniform (multifractal for E — > E c ) background $£(r). Since the eigenfunction correlations 
are described by the formula Q5.30|) , which has exactly the same form as in the Gaussian 
Ensemble, it is not surprising that the level statistics has the WD form everywhere in the 
extended phase ||183j , |184|| . For larger u, <r(r, r, E, u) is expected to decrease as u~ r, / d , where 
r] = d — D 2 , according to the scaling arguments | |185| , |186| , |187] , |188| , |189| , |131|| (see below), so 
that 

a(r, r, E, cu)/a(r, r, E) ~ {uj/A i )-^ d ; u > A ? , (5.31) 



A^ around 



up to a numerical coefficient of order unity. 

The above formulas are valid in the metallic regime, 
critical point (£ ^> L) they take the form 



i.e. for L ^> £. Exactly at the 



and 



a(r,r,E,u)/a(T,T,E) 



a{r,r,E,u)/a{r,r,E) ~ (u/A) 



u < A 



-v/d 



uj > A 



(5.32) 



(5.33) 



Of course, Eq.( [5.32[ ) is not sufficient to ensure the WD statistics at the critical point, since 
there is only of order of one level within its validity range u ~ A. Indeed, the numerical 
simulations show that the level statistics on the mobility edge is different from the WD one 
50], |1J, [179], |T80|,[T8Tj. 



8 Eq. (5.30) was explicitly derived in [ 182 1 for the case of the sparse random matrix model [183, 184 1 
corresponding to the limit of infinite dimensionality, d — 00. So, even at d = 00, where the wave function 
sparsity (multifractality) takes its extreme form, nearby in energy eigenfunctions are fully correlated. 
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However, Eq. ( p.32| ) allows us to make an important conclusion concerning the behavior 



of R(s) (or, which is essentially the same, of the nearest neighbor spacing distribution P(s)) 
at small s = uo/A. For this purpose, it is sufficient to consider only two neighboring levels. 
Let their energy difference be uoq ~ A. Let us now perturb the system by a random potential 
V(r) with (V(r)) = 0, (V(r)V(r')) = TS(r — r'). For the two-level system it reduces to a 2 x 2 
matrix {V^}, i,j = 1, 2, with elements = J d d r V (r)ty* (r)^ j(r) . The crucial point is that 
the variances of the diagonal and off-diagonal matrix elements are according to Eq. ( |5.32j ) 
equal to each other up to a factor of order of unity: 

(Vu)/(\V&\) = <r(r,T,E,u)/a(T,T,E) ~ 1 (5.34) 

The distance between the perturbed levels is given by uo = [(Vn — V22 + ^o) 2 + l^l 2 ] 1 ^ 2 - 
Choosing the amplitude of the potential in such a way that the typical energy shift Vn ~ A 
and using Eq. (|5.34j ), we find ( | V12 1 2 ) ~ A. As a result, the probability density for the level 
separation uo is for w < A of the form dP ~ (uo/A) 13 duo /A, with some prefactor of order of 
unity. We thus conclude that in the critical point 

P(s) ~ c p s p , s < 1 , (5.35) 



with a coefficient c@ of order of unity, in agreement with the numerical findings [|41], |1 79| , |18 1 



Spatial correlations at the mobility edge. Another question that one can ask con- 
cerning the wave function correlations at the mobility edge is how the correlations decay 
with distance in real space. Let us consider correlations a(r, r', E) of the eigenfunction am- 
plitudes in two different points. Matching the behavior of a(r, r', E) at r = r' with that at 
|r — r'| ~ min{L, £} (where no strong correlations is expected), we find in the metallic phase 
(L»f>Z) ' " 

a{r,r',E) ~ Zr 2d f|r - r'p ; I < |r - r'| < f (5.36) 
and at the critical point (£ ^> L) 

a(v,r',E) L- M+ "|r - r'r" ; I < |r - r'| < L . (5.37) 

The same consideration can be straightforwardly applied to the higher correlation functions 
(|'?/> 29l (r)||'?/> 2 ' ?2 (r / )|), leading to the conclusion that they decay as 

(|^ 2,?1 (r)||V' 2 ' 72 (r , )l) oc \ r - r'\- [d+T(qi)+T(q2) - T ^ qi+q2)] , 



where r(q) = (q — l)D q has been already introduced in Sec. |3.3.2| . These results were 
obtained for the first time by Wegner [|190|| from renormalization-group calculations. The 
above formulas for the spatial correlations at the mobility edge can be generalized to the 
correlator <r(r, r', E, uo) of two different eigenfunctions by using the scaling assumptions that 
this is the length L u = L(A/uj) l l d = {l/uju) l l d which will set the correlation range at finite 
uo. This give for the correlations at criticality (I < |r — r'| < < L, £) 

a(r, r', E, uo) ~ L- 2d {uou)-^ d \r - r'^ . (5.38) 

Chalker and Daniell ||185| , |186|| conjectured that the diffusion propagator 7(r, r', E, uo), Eq. ( p. 77]) , 



has the same behavior (|5.38|) at the mobility edge. Extensive numerical simulations ||185| 



187 , 188] , |189|| fully confirm this conjecture, which determines the form of the anomalous 



diffusion at the critical point of the Anderson transition. 
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5.3 Power- law random banded matrix ensemble: Anderson tran- 
sition in ID. 

The ensemble of random banded matrices (RBM) is defined as the set of matrices with 
elements 



Hi 



Gija(\i-j\) 



(5.39) 



where the matrix G runs over the Gaussian ensemble, and a(r) is some function satisfying 
the condition lim^oo a(r) = and determining the shape of the band. In the most fre- 
quently considered case of RBM the function a(r) is considered to be (exponentially) fast 
decaying when r exceeds some typical value b called the bandwidth. Matrices of this sort 
were first introduced as an attempt to describe an intermediate level statistics for Hamil- 
tonian systems in a transitional regime between complete integrability and fully developed 
chaos ||191| and then appeared in various contexts ranging from atomic physics (see [|192|1 and 
references therein) to solid state physics [18] and especially in the course of investigations 
of the quantum behavior of periodically driven Hamiltonian systems [193 , 194 1 - The most 



studied system of the latter type is the quantum kicked rotor ||104|| characterized by the 
Hamiltonian 

_ + v(9) ^ ~ mT ) > ( 5 - 4 °) 



H 



where I = —ihd/dO is the angular momentum operator conjugated to the angle 9. The 
constants T and I are the period of kicks and the moment of inertia, correspondingly, and 
V{9) is usually taken to be V(9) = kcos9. Classically, the kicked rotor exhibits an unbound 
diffusion in the angular momentum space when the strength of kicks k exceeds some critical 
value. It was found, however, that the quantum effects suppress the classical diffusion in 
close analogy with the localization of a quantum particle by a random potential [ |195| , |104 . 



It is natural to consider the evolution (Floquet) operator U that relates values of the 
wavefunction over one period of perturbation, ip{9,t + T) = Uip(9,t), in the "unperturbed" 
basis of eigenfunctions of the operator I: [/) = ( 27r )i/2 exp(znfl), n = ±0, ±1, . . .. The matrix 
elements (m\U\n) tend to zero when \m — n\ — > oo. In the case V(9) = kcos9 this decay 
is faster than exponential when \m — n\ exceeds b w k/h, whereas within the band of the 



effective width b matrix elements prove to be pseudorandom |10 1| . 

Let us note, however, that the exponentially rapid decay of (m\U\n) in the above men- 
tioned situation is due to the infinite differentiability of V(9) = cos9. If we take a function 
V(9) having a discontinuity in a derivative of some order, the corresponding matrix elements 
of the evolution operator would decay in a power-law fashion in the limit |n — m\ — > oo. 
In fact, there is an interesting example of a periodically driven system where the matrix 
elements of the evolution operator decay in a power-law way, namely the so-called Fermi ac- 
celerator [Q. The power-law (pseudo-)random banded matrices appear also in other models 
of the quantum chaos, such as a close-to-circular Bunimovich stadium ||196| , [197], |198|| or a 
Coulomb center inside an integrable billiard p9*[ . 

One may also consider the random matrix ( |5.39| ) as the Hamiltonian of a one-dimensional 
tight-binding model with long-ranged off-diagonal disorder (random hopping). A closely 
related problem with non-random long-range hopping and diagonal disorder was studied 

[199]. Effect of weak long-range hopping on the localized states in a 

TO 



numerically in Ref 

3D Anderson insulator was discussed by Levitov 
power-law hopping appear also in other physical contexts |Kj ^Ij, [202 



324], Similar models with a 
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As was shown in ||102| , |18|j , the conventional RBM model can be mapped onto a ID 
supermatrix non-linear a— model, which allows for an exact analytical solution. The same 
cr— model was derived initially ||101| , [| for a particle moving in a quasi-lD system (a wire) 
and being subject to a random potential. All states are found to be asymptotically localized, 
with the localization length 

5 = g ' 2<2 ^7 oc " 2 ; B k =± a\r)r" ■ (5.41) 

r=— oo 

However, for the case of a power-like shape of the band, 

a(r) (x r~ a for large r , (5-42) 

this derivation should be reconsidered. As explained below, this leads to a more general 
ID cr-model with long-range interaction, which is much richer than the conventional short- 
range one. In particular, it exhibits the Anderson localization transition at a = 1. We 
mainly follow Ref. |57[ in our presentation. 



5.3.1 Mapping onto the effective a-model 



Let us consider the ensemble of large N x N matrices (N — > oo) defined by Eq. ( |5.39|) with 
a function a(r) having the form 

a(r) * { \r/b)-°, r > b ■ (5 ' 43) 

The parameter b will serve to label the critical models with a = 1. We will consider b to 
be large: b ^> 1, in order to justify formally the derivation of the a— model. We will argue 
later on that our conclusions are qualitatively valid for arbitrary b as well. We will call 
the ensemble ( |5.39| ), ( |5.43| ) the power-law random banded matrix (PRBM) model. While 



considering localization properties, we will restrict ourselves by one-loop calculation, and for 
this reason, by the orthogonal ensemble. 

The PRBM model can be mapped onto the cr-model on a lattice, 

S{Q} = -\(™A ) 2 StrJ2 [(A-% - Aj%] Q(i)QU) ~ ™ ]TStrQ(z)A . (5.44) 

ij i 

Here Q(i) = T i ~ 1 ATj satisfies the constraint Q 2 (i) = 1, A is the matrix with elements 
Aij = a 2 (\i — j\), A is given byf] A = Am ~ J2^L-oo 0,2 ( r )' anc ^ v ^ s ^ ne density of states: 

(4A -£ 2 ) 1/2 . (5.45) 



2nA 

The standard next step is to restrict oneself to the long wavelength fluctuations of the 
Q— field. For usual RBM characterized by a function a(r) decreasing faster than any power 

9 The expression for Ao is valid for a > 1/2 and in the limit N — » oo. When a < 1/2, Aq starts t o hav e 
a dependence on N which can be removed by a proper rescaling of the matrix elements Hij in Eq. (5.3E). 



Then, the properties of the model with a < 1/2 turn out to be equivalent to those of the GOE, so we will 
not consider this case any longer. 
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of r as r — > oo, this is achieved by the momentum expansion of the first term in the action 
O): 

[(A-% - A^%] Q(i)Q(j) = N- 1 £ [V " Ao 1 ] Q q Q- q 

w ll^ 1 E <? 2 w-* = i^J dx {dxQ)2 ' (5 - 46) 

where Q q is the Fourier transform of Q(i), and i?2 = J2i Aki(k — I) 2 , as defined in Eq. ( |5.41| ). 
This immediately leads to the standard continuous version of the nonlinear a— model: 



S{Q} = -^Str J dx 



l -D (d x Q) 2 + tuQA 



(5.47) 



with the classical diffusion constant Dq = 7TUB2, which implies the exponential localization 
of eigenstates with the localization length £ = ttuD cx b 2 . 

Let us try to implement the same procedure for the present case of power-like bandshape 
( |5.43j ). Restricting ourselves to the lowest order term in the momentum expansion, we arrive 



again at Eq. (|5.47|) as long as a > 3/2. This suggests that for a > 3/2 the eigenstates 
of the present model should be localized in the spatial domain of the extension £ oc vDq. 
However, in contrast to the usual RBM model this localization is power-like rather than 
exponential: |V ; ( r )| 2 °c r~ 2a at r ^> £. This is quite evident due to the possibility of direct 
hopping with the power-law amplitude. On a more formal level the appearance of power-law 
tails of wavefunctions is a consequence of the breakdown of the momentum expansion for 
the function A^ 1 — Aq 1 in higher orders in q 2 . The presence of power- law "tails" of the wave 
functions, with an exponent a determined by the decay of hopping elements, was found in 
numerical simulations in Refs. [ |i99j , |197| , |198| | . 

The most interesting region 1/2 < a < 3/2 requires a separate consideration. In this 
case, Eq. ( |5.46| ) loses its validity in view of the divergence of the coefficient B 2 . We find 
instead 



AKA^-A- 1 ) « A -A q = 2 dra 2 (r)(l-cosqr) 



oc 



2 r°° dr 

= — 11 dx(l-cosx) + (b\q\) 2a ^P(l-cosar) 

\q\ [io Jb\q\ x a 

w c a b 2a \q\ 2a - 1 for 1/2 < a< 3/2 and |g| < 1/6 , (5.48) 

where c a = 2 J °° ^(1 — cosx) is a numerical constant, c a = [2T(2 — 2a) /{2a — 1)] sin7ra; 
for a^l and c\ = is for a = 1. 

The corresponding long-wavelength part of the action, 

MQ} = "^Str J (dq) \q^- l Q q Q_ q , (5.49) 

where J (dq) = J ^ = N^ 1 ^ , can not be reduced to the local-in-space form in the coor- 
dinate representation any longer. Here, 1/t — \(nv) 2 cja 2a oc b 2a ~ l ^> 1 plays the role of 
coupling constant, justifying the perturbative and renormalization group treatment of the 
model described below. 
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Let us mention that if we consider the RBM model as a tight-binding Hamiltonian, 
the corresponding classical motion described by the master equation on the ID lattice is 
superdiffusive forl/2<«<3/2 with a typical displacement in a time t being r oc 

f l/(2«-l) 

(Levy law of index 2a — 1; see Ref. |[203|| for a review). As will be discussed in Sec. |5.3.3| , this 



influences the asymptotic behavior of the spectral correlation function for the corresponding 
quantum system. 

Perturbative treatment of the non-linear a— model: General formulas. Here we 
derive one-loop perturbative corrections to the density-density correlation function and in- 
verse participation ratios. Analysis of these expressions for various values of the power-law 
parameter a will be presented in Sec. f>.3.2\ 



The density-density correlation function (diffusion propagator) 
K(r u r 2 ;uj) = G l ^ +u;/2 (r l ,r2)G^~ ul/2 (r 2 ,r l ) 
can be expressed in terms of the a— model as follows ||: 

K(r u r 2 ; u) = -(tiu) 2 [ DQ QnM^k^iMn)^^ • (5-50) 



Here the indices p,p' of the matrix Q pp ' !a /3 correspond to its advanced-retarded block struc- 
ture, whereas a, j3 discriminate between bosonic and fermionic degrees of freedom. The 
matrix kpp is equal to 1 for bosons and (-1) for fermions. To calculate the correlation 
function ( |5.50| ) perturbatively, we use here the following parametrization |J of the matrix 
Q: 

, / W 2 W 4 \ 

Q = A(W + Vl + W 2 ) = A f 1 + W + — §- + •••) ' ( 5 - 51 ) 

where W is block-off-diagonal in the advanced-retarded representation. To get the pertur- 
bative expansion for K(ri,r 2 ;uj), one has to substitute Eq. ( |5.51| ) into ( |5.50|) , to separate 



the part quadratic in W from the rest in the exponent and to apply the Wick theorem. In 



the usual case, when the action is given by Eq. ( |5.47| ), the leading order (tree level) result 
reads in the momentum space as follows: 

K (q,u)= - 2 Y . . (5.52) 

The perturbative quantum corrections do not modify the general form (|5.52f ), but change 
the value of the diffusion constant. In particular, in one- loop order one gets Eq. (|5.52|) with 



D Q replaced by fl204 | 



D = D { 1 - — - n 2 ■ } ■ ( 5 - 53 ) 

y qi=-Kni/Li J 

This induces the standard weak-localization correction to the conductivity. 

Now we implement an analogous procedure for the non-local a— model of the type of 
Eqs. O, O: 

S{Q} = -Str U ( r ~ r')Q(r)Q(r') - i^- ^ StrAQ(r) , (5.54) 
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with U(r) oc r 2a as r — > oo, so that the Fourier-transform of U(r) behaves at small momenta 

as 

(j( q ) = 1/2 < a < 2 . (5.55) 

The exponent a is related to the parameter a of the RBM model by a = 2a — 1. In leading 
order, we keep in the action the terms quadratic in W only, which yields: 

2ttp 



K (q,uj) 



(5.56) 



8(7ri/t)- 1 |g| fT - iu 

corresponding to a superdiffusive behavior. 

To calculate the one- loop correction to K (q) (we set uj = for simplicity) we expand 
the kinetic term in S{Q} up to fourth order in W: 

Str U{r - r')Q(r)Q(r') | 4th order = £ jUtr W 2 (r)W 2 (r') . (5.57) 



The contraction rules are given by Eqs. (8), (16) of Ref. |25j (reproduced below as Eq. (|2.21|) 
for the unitary symmetry), with the propagator 11(g) replaced by 

U(q) = -4-r • (5-58) 
w 8|g| CT v ' 

Evaluating the one-loop correction, we get the following expression for the density-density 
correlation function up to the one-loop order: 

(vrz/) 2 f, „s \q + k\ a - \k\ c 



K~\q)=K Q \q) 



(dky 



\k\ a 



(5.59) 



Now we calculate the perturbative correction to the inverse participation ratios P q . The 
results of Ref. presented in Sec. |3.3.1| are straightforwardly applicable to the present 
case of power-law RBM, provided the appropriate modification of the diffusion propagator 
entering the contraction rules is made, see the text preceding Eq. ( |5.58j ). One finds 



(P q ) 



ji + ^g(g-l)X>(r,r)j 



(2g-l)H 



-i 



(5.60) 



where H(r, r') = (1/N) J2 Q ^( ( l) exp[ig(r — r')] and 11(g) is given by Eq. ( p.58|) . 



Renormalization group treatment. The effective a— model, Eq. ( |5.54j) , is actually of 
one-dimensional nature. However, for the sake of generality, we find it convenient to con- 
sider it here to be defined in <i-dimensional space with arbitrary d. The form (|5.56|) of the 
generalized diffusion propagator implies that d = o plays the role of the logarithmic dimen- 
sion for the problem. In the vicinity of this critical value we can carry out a renormalization 
group (RG) treatment of the model, following the procedure developed for general non-linear 
a— models in |[205| , [206] , [20 7| , |208|| . We begin by expressing the action in terms of the renor- 
malized coupling constant t = Z{ 1 tBfi d ~° where is the bare coupling constant, Z\ is the 
renormalization constant, and /i -1 is the length scale governing the renormalization 



S 



U(\r - r'\) Str \-W(r)W(r') + ^1 + W 2 (r) y/l + W 2 (i 



2tZ l 

rr 

^StrVl + W 2 {r 



VKVUJ 



4 



(5.61) 



10 Note that the W^-field renormalization is absent due to the supersymmetric character of the problem, 
which is physically related to the particle number conservation |p[ . 
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Expanding the action in powers of W(r) and keeping terms up to 4-th order, we get 
S = So + Si + 0(W 6 ) , 

S = (\r - r' \) Six (W {r) - W (r')) 2 -'-^^StiW 2 (r) , 



Si 



AtZ l 

/j, d - a 
8tZ 1 



U(\r - r'|)Str W 2 {r)W 2 {r') + Str W\r) . 

32 



(5.62) 



We have restricted ourselves to 4-th order terms, since they are sufficient for obtaining the 
renormalized quadratic part of the action in one-loop order. The calculation yields, after the 
cancelation of an J (dk) oc 5(0) term with the contribution of the Jacobian: 



Vad = S + (SJ = | / (dq)StvW q W- q 



ft 



d—a 



tz 1 



■U(q) 



(dk) 



U(k) - U(k + q) 



(5.63) 

According to the renormalization group idea, one has to chose the constant Z\(t) = l+at+. . . 
so as to cancel the divergence in the coefficient in front of the leading \q\ a term. 



5.3.2 Localization and eigenfunction statistics. 

In this subsection we analyze the model in different regions of the exponent a = 2a — 1 using 
the general formulas derived in Sec. |5.3.2 . 



Localized regime: l<a<2 (1 < a < 3/2). To evaluate the one-loop correction ( 5.59|) 
to the diffusion propagator, we use the expansion 



qk a q 2 



|q+ k 


a 




k 


cr 



- 1 



a F + 2¥ + a {2- 1 ) 



a 



qk 



+ 



\k\ a 



q < k 
q~>k 



(5.64) 



Thus, the integral in Eq. ( 5.59|) can be estimated as 



(dk) 



|q 


+ 


k 


a 




k 


a 



k>q 



k<q 



\k\ a 



(5.65) 



and we are interested in the particular case of the dimension d — 1. For a > d = 1 the 



integral diverges at low k, and the second term in Eq. (|5.65 ) dominates. This gives 



^ _1 i?r 



(5.66) 



where L is the system size determining the infrared cut-off (in the original RBM formulation 
it is just the matrix size N). This leads to the following one-loop expression for the diffusion 
correlator 2 _ 

K(q) = ; r 1 = r 1 - const If' 1 . (5.67) 



Now we turn to the renormalization group analysis, as described in Sec. |5.3.1| . In the 
one-loop order, the expression for the renormalization constant Z\ is determined essentially 
by the same integral J, Eq. (|5.65|) , with the RG scale \i playing the role of the infrared 
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cut-off (analogous to that of the system size L in Eq. ( |5.67| )). This yields (in the minimal 
subtraction scheme): 

Z 1 (t) = l-±-^- + 0(t 2 ), (5.68) 

Z7T — 1 

leading to the following relation between the bare and the renormalized coupling constant 
analogous to Eq. G5.67Q : 

V - * = — = - - -—V 1 - 17 ■ (5-69) 
t ts ts 2a — 1 



From Eq. ( |5.68|) , we get the expression for the /3-function: 

at 



Pit) 



= , J* . % m = -(* -l)t-f + 0(e) . (5.70) 



Both Eqs. (|5.67|) and ( 5.70|) show that the coupling constant t increases with the system size 



L (resp. scale /i), which is analogous to the behavior found in the conventional scaling theory 



of localization in d < 2 dimensions ||209| , |167j| . The RG flow reaches the strong coupling regime 
t ~ 1 at the scale /i ~ t 1 ^ . Remembering the relation of the bare coupling constant ts 
and the index a to the parameters of the original PRBM model: t^ 1 oc 6 2a_1 ; a oc 2a — 1, 
we conclude that the length scale 



f ~ t B 1/{a - 1] ~ (5.71) 

plays the role of the localization length for the PRBM model. 

This conclusion is also supported by an inspection of the expression for the IPR, Eq. 
( 5.60Q . Evaluation of the one- loop perturbative correction in Eq. ( |5.60|) yields 



2a-l 



where ((a) is Riemann's zeta-function; £(cr) ~ l/(o" — 1) for o close to unity. The correction 
term becomes comparable to the leading (GOE) contribution for the system size N ~ t l /^^ l \ 
parametrically coinciding with the localization length £. For larger the perturbative 
expression ( |5.72| ) loses its validity, and the IPR is expected to saturate at a constant value 
(P q ) ~ ^ q for N > C- 

In concluding this subsection, let us stress once more that the localized eigenstates in 
the present model are expected to have integrable power-law tails: |V ;2 ( r )| oc r~~ 2a = r~ a ~ l 
at r > £. 

Delocalized regime: < a < 1 (1/2 < a < 1). We begin again by considering the 
perturbative corrections to the diffusion propagator fl5.59Q . The first term in the r.h.s. 



of Eq. (|5.65|) proportional to J dk/k 2 is determined by the vicinity of its lower cut-off (i.e. 
k ~ q), whereas the second one is proportional to J dk/k a and thus determined by the vicinity 
of its upper cut-off (i.e. again k ~ q). Therefore, the integral (|5.65|) is now dominated by 
the region k ~ q, and is proportional to \q\. We get, therefore: 



(nv) 2 K-\q) = ^-C a \q\ (5.73) 
with a numerical constant C a . 
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We see that the correction term is of higher order in \q\ as compared to the leading one. 
Thus, it does not lead to a renormalization of the coupling constant t. This is readily seen 
also in the framework of the RG scheme, where the one- loop integral in Eq. ( 5.63| ) does not 



give rise to terms of the form \q\ a . One can check that this feature is not specific to the 
one-loop RG calculation, but holds in higher orders as well. For the case of a vector model 
with long-ranged interaction this conclusion was reached by Brezin et al [(210|| . Thus, the 



renormalization constant Z\ is equal to unity and the (3— function is trivial: 

P(t) = (1 - a)t . (5.74) 

This means that the model does not possess a critical point, and, for a < 1, all states are 
delocalized, for any value of the bare coupling constant t. This property should be contrasted 
with the behavior of a d— dimensional conductor described by the conventional local non- 
linear cr— model and undergoing an Anderson transition at some critical coupling t = t c 



T67|. 



Though all states of the model are delocalized, their statistical properties are different 
from GOE. In particular, calculating the variance of the inverse participation ratio Pi (see 
Sec. 13X31 , Eq. (ggg) ) we get: 



^ ' rr' q=nn/N;n=l,2,.. 

At a > 1/2 the sum over q is convergent yielding 

Thus, in this regime the fluctuations of the IPR are much stronger than for the GOE where 
6{P 2 ) oc 1/N. Only for a < 1/2 (a < 3/4) the IPR fluctuations acquire the GOE character. 
Considering higher irreducible moments (cumulants) of the IPR, ((-Rj 1 )), one finds that the 
GOE behavior is restored at a < = 1/n. In this sense, the model is analogous to 
a (i-dimensional conductor at d = 2/a. Therefore, only when o — > (correspondingly, 
a — > 1/2 in the original PRBM formulation) all statistical properties become equivalent to 
those characteristic of GOE. 



Critical regime: a = 1 (a = 1). As we have seen, a = 1 separates the regions of localized 
(a > 1) and extended (cr < 1) states. It is then natural to expect some critical properties 
showing up just at cr = 1. Let us again consider the generalized diffusion propagator, Eq. 
( |5.59| ). At a = 1 the one- loop correction yields 



7TZ/ 



; K-\q)=A\q\ 



t-'-^HlqlL) 



(5.77) 



As expected at the critical point, the correction to the coupling constant is of logarith- 
mic nature. However, Eq. ( 5.77 ) differs essentially from that typical for a 2D disordered 
conductor: 



t 



-i 



t 



-i 



ln(L//) 



(5.78) 



where the bare coupling constant ts corresponds to scale /. Comparing the two formulas, 
we see that in Eq. (|5.77|) the mean free path / is replaced by the inverse momentum q" 1 . 
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Therefore, the correction to the bare coupling constant is small for low momenta q ~ 1/L, and 
the correlator K(q) is not renormalized. This implies the absence of eigenstate localization, 
in contrast to the 2D diffusive conductor case, where Eq. ( p. 78 ) results in an exponentially 
large localization length £ oc expt^ 1 . On a more formal level, the absence of essential 
corrections to the low- 9 behavior of K(q) is due to the fact that the region k > q does 
not give a logarithmic contribution. This is intimately connected with the absence of t 
renormalization at a < 1. 

To study in more details the structure of critical eigenfunctions, let us consider the set 



of IPR P q . The perturbative correction, Eq. (|5.60|) , is evaluated at a = 1 as 



q(q 



1) ln(N/b) 

07T 



(2q 



(5.79) 



where the microscopic scale b, Eq. Q5.43Q , enters as the ultraviolet cut-off for the a— model, 
the role usually played by the mean free path I. This formula is valid as long as the correction 
is small: q <C In (N/b)\ . For larger q the perturbation theory breaks down, and one 
has to use the renormalization group approach. This requires [|5"0| p8fl introduction of higher 
vertices of the type z q J Str q (QkA)dr into the action of the non-linear a— model and their 
subsequent renormalization. The resulting RG equations for the charges z q read, in the 
one-loop order, 

dz n . » t 

-- q{q- l )^ z q , 



d In /i 1 ' 87T 

where is the renormalization scale. Integrating Eq. (|5.80|) , we find: 



TO 



(2q 



Ni- 



9(9-1) - 



(5.80) 



(5.81) 



Note, that this formula is reduced to the perturbative expression, Eq. (|5.79| ), in the regime 
g« [^ln(iV/6)]- 1/2 . 

The behavior described by Eq. ( |5.81| ) is characteristic of a multifractal structure of wave 
functions (see Sec. |3.3.2| , p^2[) . Comparing ( p.81| ) with Q3.60| ), we find the multifractality 
dimensions D q = 1 — qt/87T. The general formula valid for both orthogonal (/3 = 1) and 
unitary (/3 = 2) universality classes is 

(5-82) 

This form of the fractal dimensions is reminiscent of that found in two and 2+e dimensions for 
the usual diffusive conductor, see Eq. ( |3.63| ). The one-loop result ( |5.82| ) holds for q 8n/t. 

Now we can understand the reason for the g-dependent logarithmic correction to the diffu- 
sion propag atorK(g), Eq. ( gTrfl) . As was mentioned in the end of Sec. |5.2| , the multifractality 
of eigenfunctions determines the momentum dependence of the diffusion propagator at high 
q in the critical point [|185| , |186|| 

K-\q)^\q\{\q\L)^ ; rj = d - D 2 (5.83) 

(at finite frequency u the system size L is replaced by L w ~ (D/cu) 1 ^ 2 ). On the other hand, 
the logarithmic correction in Eq. ( |5.77| ) is the first term of the expansion 

-1 2 



-,K-\q) 



1- ±\n(\q\L) + C 2 

87T 



t 



in 



H\q\L) 



+ 



(5.84) 
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Figure 3: Typical eigenfunctions for the matrix size N = 800 and four different values of a: 



0.375, 0.875, 1.250, and 1.625. From 57 . 



where F(x) is some parameterless function. Since rj = 1 — D 2 — t/Air, Eq. (|5.84 ) has precisely 
the form expected from (|5.83|) , assuming that F(x) ~ e~ 2x at x ^> 1. 

The set of fractal dimensions D q (as well as spectral properties at a = 1, see Sec. |5.3.3|) 
is parametrized by the coupling constant t. Strictly speaking, the above a— model derivation 
is justified for t <C 1 (i. e. b 1). However, the opposite limiting case can be also studied, 
following Levitov |[200| , |201| , |324|| . It corresponds to a d— dimensional Anderson insulator, 
perturbed by a weak long-range hopping with an amplitude decreasing with distance as 
r~ a . The arguments of Levitov imply that the states delocalize at a < d, carrying fractal 
properties at a — d. The PRBM model in the limit b <C 1 is just the ID version of this 
problem. This shows that the conclusion about localization (derealization) of eigenstates 
for a > 1 (resp. a < 1), with o = 1 being a critical point holds irrespective of the particular 
value of the parameter b. The 5 < 1 limit of the PRBM model was also studied in Ref. [|B"T| 
where the IPR P 2 was calculated, yielding the fractal dimension D 2 oc b. Alternatively, 
the regime of the Anderson insulator with weak power-law hopping can be described in the 
framework of the non-linear a— model, eq. (|5.54j) , by considering the limit t > 1. Formally, 
the non-linear a— model for arbitrary t can be derived from a microscopic tight-binding 
model by allowing ii>l "orbitals" per site [ |167|| . 

To summarize, the PRBM model ( |5.39| ), ( p.43j ) with < b < oo or the er— model ( |5.54| ) 
with a coupling constant < t < oo represent at a = 1 a continuous family of critical 
theories parametrized by the value of b (respectively, t). 
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Figure 4: Index v characterizing the dependence of the inverse participation ratio (P2) on 
the matrix size N via (P2) oc 1/N U , as a function of a. Points refer to the best-fit values 
obtained from matrix sizes between N = 100 and N = 1000 (squares) or N = 2400 (circles). 
The dashed line is the theoretical prediction for the transition from v — 1, at small a, to 
v = 0, at large a. From |57 |. 



Numerical simulations. Numerical simulations of the PRBM model were performed in 

G [0,2] and 6 = 1. In Fig. |] we present typical eigenf unctions 



Ref. 57 



for values of a 

for four different regions of a. In agreement with the theoretical picture presented above, 
the eigenstates corresponding to a = 0.375 and a = 0.875 are extended, whereas those 
corresponding to a = 1.25 and a = 1.625 are localized. At the same time, one can notice 
that the states with a = 0.875 and a = 1.25 exhibit a quite sparse structure, as opposed to 
the other two cases. This can be explained by the proximity of the former two values of a to 
the critical value a = 1.0 where eigenstates should show the multifractal behavior. To get a 
quantitative information about the properties of the eigenstates, the mean value of the IPR, 
(P2), and its relative variance, 5 = ((Pi) — (P2) 2 )/ (P2) 2 , were calculated in [fTfll . At any given 
a the dependences of the quantities (P2) and S on the matrix size N were approximated by 
the power-laws (P 2 ) oc 1/N U , 5 oc l/iV 7 * for iV ranging from 100 to 2400. In Figs. £| and 
|5] the values of the exponents v and \i obtained in this way are plotted versus the PRBM 
parameter a. The expected theoretical curves are presented as well. We see from Fig. ^ 
that the data show a crossover from the behavior typical for extended states (v = 1) to that 
typical for localized states (u — 0), centered approximately at the critical point a = 1. The 
deviations from the sharp step-like theoretical curve v(a) can be presumably attributed to 
the finite-size effects which are unusually pronounced in the PRBM model due to the long- 
range nature of the off-diagonal coupling. The data for the exponent fi (Fig. ||) also show a 
reasonable agreement with the expected linear crossover, = 4(1 — a) for 3/4 < a < 1, see 
Eq. (grp . 

While the presented data are in good agreement with the above theoretical picture, a 
more detailed numerical investigation of the structure of eigenstates and of spectral statistics 
is certainly desirable. In particular, it would be especially interesting to study the critical 
manifold, a = 1, where the multifractal properties of eigenstates and intermediate level 
statistics are predicted by the theory. 
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Figure 5: The same as Fig. f|, but for the index //, derived from the N dependence of the 
variance 8 of the inverse participation ratio: 5 = ((P|) - (i^) 2 )/^) 2 oc The dashed 

line corresponds to the predicted linear crossover from /x=lata<3/4to/i = 0ata>l. 
From [57|. 



5.3.3 Spectral properties. 

Let us consider now the spectral statistics of the PRBM model. In the "metallic" regime, 
the leading correction to the Wigner-Dyson form R\y D (s) of the level correlation function 
is given by the properly modified Eq. fl2.3U[ ), i.e. 



R(s) 



1 n ^ 2 

1 + 2^ S 



Rwd(s) 



(5.85) 



where 



C 



n 2 (g) 



rr' q=-nn/N; n=l,2,... 

a > 1/2 
-AT 1 , a < 1/2 



t 2 
64tt 2ct 
const 



2a~2 



i 1 



(5.86) 



647r 2 CT& l-2<7- 



At a < 1/2 the sum divergent at high momenta is cut off at q ~ tt/6; the procedure leaving 
undetermined a constant of order of unity. The correlation function R(s) is close to its RMT 
value if a < 1 (the region of delocalized states), or else if a > 1 and the system size iV is 
much less than the localization length £, Eq. ( |5.71| ). Under these conditions, Eq. (|5.85|) 
holds as long as the correction term is small compared to the leading one. This requirement 
produces the following restriction on the frequency s = u/A: 



s < s r 



t- l N l - a , 

t-'b^N 1 / 2 oc (Nb) 1 / 2 



a > 1/2 
a < 1/2 



(5.87) 



At larger frequencies (s > s c ), the form of the level correlation function changes from 
the 1/s 2 behavior typical for RMT to a different one, in full similarity with the Altshuler- 
Shklovskii regime (|2.33j ) in the case of a conventional d- dimensional conductor. Extending 



77 



Eq. ( [2.28D to the present case, we find 



K \ t^" 1 ^" 1 oc (JV&)" 1 , a<l/2 ' 

At last, let us consider the level statistics at criticality, a — 1. In this case the coefficient 
of proportionality in the asymptotic expression (|5.88|) vanishes in view of analyticity: 



M , At f°° dx . 
R {c) (s) — / ^ = . 5.89 

This is similar to what happens in the case of a 2D diffusive conductor (see Sec. |5.1p . A more 
accurate consideration requires taking into account the high-momentum cut-off at q ~ b^ 1 . 
In full analogy with the 2D situation mentioned we find then a linear term in the level 
number variance: 

(5N 2 {E)) ~ X (N(E)) ; X = J R {c \s)ds = ^- . (5.90) 

The presence of the linear term ( |5.90|) (as well as the multifractality of eigenfunctions, Sec. 
5.3. 2| ) makes the case a = 1 similar to the situation on the mobility edge of a disordered 



conductor in d > 2 discussed in Sec. pT| , |5T2| . Let us finally mention that the value of x-> Eq. 
( |5.90| ), is in agreement with the formula ( |5.27| ), where D 2 is given by Eq. fl5.82| ) with q = 2, 
and d = 1 in the present case. We stress, however, that both Eqs. (|5.82 ) and Q5.90 ) have 



been derived in the leading order in t/87r <C 1. As we discuss below, at t/8ff > 1 Eq. ( J5.271 ) 
appears to be violated. 

In fact, one can calculate the whole two-level correlation function using the results of 
Refs. [0, [ll[] presented in Sec. |2.2| . We consider the unitary symmetry (/3 = 2) here for 



simplicity. In this case, the level correlation function is described by a single formula R(s) = 
Rpert(s) + R osc (s) [with R pert (s) and R osc (s) given by Eqs. ( |2.38[ ), fl2.34| ) respectively] in the 



whole range of s. The spectral determinant D(s), Eq. fl2.37|) can be easily calculated [^T 
yielding 

1 (TTSt/16) 2 

sinh 2 (7rs£/16) 
in the case of periodic boundary conditions and 



D(s) = ^.7: (5.91) 



D(s) = \ (5.92) 
v ; s 2 sinh(7rst/8) v ; 



in the case of hard-wall boundary conditions. The resulting expressions for R(s) are 

R (c) (s) = sin 2 7TS (TTgt/16) 2 

{S> (tts) 2 sinh 2 (7rst/16) 



(5.93) 



for periodic boundary conditions and 
R^(s) 



An 2 s 2 



(nst/8) 2 . . TTst/i 

1 + V / ; , - 2 cos (tts) 



sinh 2 (7rst/8) sinh(7rst/8) 



(5.94) 
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for the hard- wall boundary conditions. Calculating the spectral compressibility \ = J R^(s)ds, 
we reproduce Eq. (|5.90| ). 

There exists a deep connection between the PRBM ensemble and the random matrix 



ensemble introduced by Moshe, Neuberger and Shapiro in Ref. |[212|| . This latter ensemble 
is determined by the probability density 



V(H) oc J dUexp{-TrH 2 -h 2 N 2 Tr([U,H)[U,H]*)} 



(5.95) 



where J dU denotes the integration over the group of unitary matrices with the Haar measure. 
The parameter h allows to interpolate between the limits of the Wigner-Dyson (h = 0) and 
the Poisson (h — > oo) statistics. The relation between the ensemble ( |5.95|) and the PRBM 
ensemble is established as follows ||211|| . For any given unitary matrix U the level statistics 
in the ensemble determined by the density ( |5.95| ) (without integration over U) is determined 
by the eigenvalues e tdk (k = 1, . . . , N) of U. For iV> la typical matrix U has an essentially 
uniform density of eigenvalues, so that we can consider as a typical representative the matrix 
with equidistant eigenvalues, 9k = 2nk/N. Indeed, weak fluctuations of 9k around these 
values do not affect essentially the behavior of a(r) in Eq. ( |5.97| ) below. On the other hand, 
the matrices with strongly nonuniform density of eigenvalues represent a vanishingly small 
fraction of the whole group volume and can be neglected. Making a transformation to the 
eigenbasis of U, we can transform V(H) to 



V{H) ex exp j- \ H -^\ 2 [l + ( 2Nh ) 2 sin2 ~ j)} | , 
so that the matrix H has a form of the PRBM in the critical regime a = 1 with 



(5.96) 



a 2 (r) 







4 -I | 1 sin 2 (7rr/Af) 
1 ~+~ *3 (tt/N) 2 







1?' 
1 



4 1 + (r/b)- 



for r <^ N 



(5.97) 



and b = l/2nh. Calculating (for h <?C 1) the coupling constant in the center of the band 
(E = 0) according to Sec. |5.3.1| , we get t = 4/6 = 8irh. 

A nice feature of the ensemble (|5.95|) is that at (3 — 2 the level correlation function 
can be calculated exactly for arbitrary value of h [|212j1 (we remind that the above cx-model 
study of the PRBM ensemble was restricted to the regime 1/27T& = h <C 1). After the 
integration over the unitary group with making use of the Itzykson-Zuber formula, the joint 
probability distribution of the eigenvalues of H is found to be equal to the probability 
density of coordinates of a system of non-interacting ID fermions in a harmonic confinement 
potential at finite temperature. The parameter h of the model ( |5.95| ) determines, in this 
formulation, the ratio of the temperature to the Fermi energy, h = T/Ep, or, in other words, 
the degree of degeneracy of the Fermi gas. Calculating the two-particle correlation function 
for non-interacting fermions, one finds the sought level correlation function | 212 |. In the 
center of the band the result reads 



(u{-w/2) V (u/2)) - M0)> 2 = - 



Nh 



dp 



a 2iphNu) 



(5.98) 



79 



where C = [e l l h — 1] 1 and the mean level spacing A is given by 

The spectral compressibility x is found from ( |5.98| ), ( |5.99| ) to be 

fdp(l + Ce p2 )- 2 

x-'-JWTc^ (5100) 

The formulas ( 5.98|) , ( |5.99|) can be simplified in the limits h <C 1 and h ^> 1. For h <C 1, 
when the gas is strongly degenerate, one gets (as usual, we introduce s = u/A) 

,s, , /sin7rs\ 2 / n 2 hs/2 \ 2 . 
#( c ) ( s ) = - 5.101 

v ; V ™ / V sinn7r W 2 / 

in precise agreement with the result ( |5.93|) obtained directly for the PRBM ensemble (with 
the above identification of the parameters of the two models, t = 87th). 

In the opposite limit, h > 1, the gas is almost classical and the correlations are weak, 

R( c \s) = - e - 2 ^ 2 * 2 . (5.102) 

Note that the corresponding spectral compressibility x ~ 1 — 1/ y2h approaches unity in the 
limt h — ► 00, whereas the formula ( 5.27]) would imply x < 1/2. Therefore, at least for the 
PRBM model, Eq. ( |5.27| ) is not an exact relation, but rather an approximation valid in the 
close-to-RMT regime 1/27T& = /i<Cl. 



6 Conductance fluctuations in quasi-one-dimensional 
wires 

This section is devoted to a study of the conductance fluctuations of a quasi-one-dimensional 
disordered system :HJ. As was already mentioned in Section |3|, there exist different 
microscopic models which can be mapped onto the same ID supersymmetric cx-model and 
thus belong to the same "quasi-one- dimensional universality class". Our treatment below 
will be based on the Iida-Weidenmuller-Zuk (IWZ) model |63| representing a wire as a 
sequence of coupled N x N random matrices, the first and the last of which are coupled 
to the states propagating in the leads. Using the multi-channel Biittiker-Landauer formula 
213| , |214| , |215| , |216| |, the mean conductance and its variance can be expressed in terms of end- 
to-end correlation functions of the cr-model. Similarly to what was done in Sec. |3.2j, where 



the eigenfunction statistics was studied, these correlation functions can be calculated exactly 
via the transfer-matrix method. This means that calculation of the functional integral can 
be reduced to solution of a "Schrodinger equation", which can be found in terms of an 
expansion in corresponding eigenfunctions. The results depend on a single parameter L/£ 
(where L is the sample length and £ the localization length) and do not depend on details 
of the underlying microscopic model. On top of the IWZ-model and two other models from 
the quasi- ID class already discussed in Sec. |3.2| (the thick wire model of Efetov and Larkin 
and the random banded matrix model), we mention here a system of weakly coupled 2d 
layers in strong magnetic field |l|, [2T|, |2(J |21|, |22|, |2||, |2|, |§. In th e quantum 
Hall regime the transport in each layer is due to edge states. Tunneling between the layers 
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leads to appearance of the transport in the transverse direction. The coupled edge states 
propagating on a surface of the cylinder form a 2d chiral metal, which can be described 
by a directed network model introduced by Chalker and Dohmen [217]. Mapping of this 
problem onto the supersymmetric spin chain was done in [ 218 , 221]; further mapping onto 



the supersymmetric a-model was presented in |[222|| . If the number of layers is sufficiently 
large, the system is of quasi-one-dimensional nature. Very recent numerical study of the 
directed network model |68|] showed perfect agreement with the analytical results presented 
below. First experimental realization of the multilayer quantum Hall system has been also 
reported recently ][224| . 



6.1 Modeling a disordered wire and mapping onto ID cr-model. 



The presentation below is based on Ref. ||65|| . We consider a quasi one-dimensional disordered 
wire of length L, decomposed into K boxes with linear dimension i. With each box i 
{% = 1, . . . , K) we associate N electronic states (/i = 1, . . . , N). The boxes 1 and K are 
coupled to ideal leads k on the left (k = L) and right (k = R) side of the disordered region. 
In each of these leads we have a number of modes \E, a, k) (a = 1,...) with transverse 
energy e a whose longitudinal momentum k is defined by their total energy E being equal to 
e a + h 2 k 2 / (2m*). Here m* is the effective mass of the electrons. We work at zero temperature. 
The Hamiltonian of the system reads 

H = / dE \E,a,K)E{E,a,K\+ J] \ifj,)H%{ju\ 

K,a ta i,j,fx,v 

+ / dE(\E,a,K)Wl(E,K)(i t x\+c.c.) . (6.1) 

i=l,K fi,a,K Jea 

The iV x iV-matrices {H™ u ) are taken to be members of the Gaussian unitary (GUE), orthog- 
onal (GOE) or symplectic (GSE) ensemble. For definiteness, we will consider the orthogonal 
symmetry case in the present section; the results for all the three symmetry classes will be 
presented in Sec. |B~2| . In the GOE case, the elements i?" are independent real Gaussian 
random variables with zero mean value and 



,2 



(H^H% U ) = — 2 . (6.3) 



A 2 

( H »V H ™>V>) = (<W*V + <W^')- ( 6 - 2 ) 

States in adjacent boxes are coupled by another set of Gaussian random variables with zero 
mean values and 

iV 2 ' 

All other matrix elements H l J v vanish. The coupling between channel and box states is 
effected by the matrix elements W*(E, k). We assume that they do not depend on E and 
k, that they obey the symmetry W* = = and that they fulfill the orthogonality 
relation 

«^Wl ll W* b = 8 0b x {i = l,K) (6.4) 

with x a normalization constant. Eq. ( |6.4| ) is convenient but does not result in any loss of 
generality p. We note that vanishes unless i — l,K. 
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The conductance of the system is given by the many-channel Landauer-Biittiker formula 

\ 2W{ , pTi m, m, 

The S-matrix of the IWZ-model reads ||225|| 

Sg = $ ab S KK ' - 2m £ ^(^p- 1 )^,^), (6.6) 

where 

= E5 ij 5^ v -H% + ^J^V 1 + in£,W* (6.7) 

iVrf, 

^ = ^EW (»* = !,*)■ (6-8) 

a=l 

The S'-matrix in Eq. ( |6.5| ) is taken at the Fermi energy E = Ep, and denotes the number 
of open channels at this energy. It follows that in order to obtain (g) and (g 2 ) we have to 
calculate the ensemble-average of a product of two resp.four propagators D^ 1 . Further 
manipulations are completely analogous to those outlined in Sec. |2.1| for the case of the 
level correlation function. We define a supervector ip = (S 1 , Sf, Xi, Xi-> £>h X2, X2) °f rea ^ 
commuting (S) and complex anticommuting (x) variables, with the first four components 
corresponding to the retarded and the last four to the advanced sector. Then 

{D~ l )% = J dty\ (S\)i{S\% exp {y^D^^j (6.9) 

where D = diag(L>, D, D, D, D\ D\ D\ D^), A = diag(l, 1, 1, 1, -1,-1,-1, -1). Products 
of two and four propagators can be expressed in similar fashion (see Eq. fl6.14|) below). After 
averaging of Eq^^) over the Gaussian distribution of random variables entering D, we 
perform the Hubbard-Stratonovitch transformation introducing 8x8 supermatrices TZi (i = 
1 . . . K) conjugate to the dyadic product tpiipl and then take the large- N limit. As a result, 
the integration over TZi is restricted to solutions of the saddle-point equation TZi = X 2 (E — 
TZi) -1 , which have the form 



E / E^ 

TZ= — -id X 2 - —TAT' 1 = a ■ I -id ■ Q , (6.10) 
with Q = TAT- 1 , Q 2 = 1. As a result, Eq. @ takes the form: 

((D- 1 )^ = J D[Q] expf^X; StrQ^ 

x J dty] (S{)i(S{% exp ^A^iE + iflA + iS&A 1 ^ 

= ((Sl)i(Sl)»)Q- (6-11) 

Here (E)% = ES^P* , = (fi 1 )^ 1 ^ 1 + [Q,*)^* V K and [Q)% = 5^Q^. The last 

line of Eq. ( |6.11| ) introduces a short-hand notation (. . .)q. 

According to Eqs. Q6.5|) , (|6.6|) , the first two moments of the conductance are given by 



<*> = 8 E tf^MD-i)${D-^)%) , (6.12) 
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(<7 2 > = 64 £ n^^^ n ?P<(^" 1 )^( i3 " lt )iv(^" 1 )^P" lt )^>- ( 6 - 13 ) 

The averages of products of the Green's functions relevant for (g) and (g 2 ) can be expressed 
similarly to Eq. fl6.11| ) as follows: 

(P-^P-^J*) = ((^)i(^)f(5S)J(5i)W gi 

((5 1 1 )i(5 1 1 )f(5 2 1 )J(5 2 1 )i,(5i 2 )^(5i 2 ) p X (5 2 2 )?(5 2 2 )^) Q . (6.14) 

For simplicity, we set = below. Evaluating the contractions in ( |6.14| ), performing the 
convolution with the projectors Q 1 , Q K in (|6.12|) , ( |6.13|) and using N c ^ ^> 1, one comes 



after some algebraic manipulations to the following a-model representation^] of (g) and (g 2 ) 

n 2 r 

(9) = -^J D[Q] QfQ 5 K eMS[Q]} , (6.15) 

(g 2 ) = D[Q] QfQfQxQx exp{-S[Q}} , (6.16) 

S[Q] = §- ( Yl StrQiQ^i + ^p StrA{Q 1 + Q K ) . (6.17) 

i 

Here = N C ^T , with T = 45x/(S + x) 2 being the so-called transmission coefficient, and 
£ = (4t> 2 /A 2 )£. In the weak-disorder limit, £ 3> the discrete sum — (£/8£) ^ StrQjQ i+ i 
can be replaced by an integral 

* fL dz Strf^j , (6.18) 



16 y o \ dz 

which is precisely the continuous version of the ID cr-model (see Sections |3p with £ being 
the localization length equal to £ = 2nuAD for a thick wire. The second term in Eq. (|6.17|) 
containing Qi and Qk describes the coupling of the wire to the leads. 

Let us note that the duplication of the number of field components in the supervector 
ip (which is forced by the supersymmetric formulation) gave us enough flexibility to write 
down the expressions Eq. ( |6.14| ) for products of two as well as of four propagators. As a 
result, we were able to express both (g) and (g 2 ) in terms of correlation functions of the 
usual "minimal" cx-model. We did not have to introduce Q-matrices of larger dimension. 
This fact no longer holds true when higher moments of the conductance are to be calculated. 
For that purpose we would have to introduce supervectors of larger size to express higher- 
order products of propagators in a form analogous to ( |6.14 ). The increase in the number 



of supervector components would force us to deal with Q-matrices of larger dimension. 



This enlargement, while posing no difficulty to perturbative calculations p8[ , gives rise to 
substantial complications in the case of the exact treatment. 

To introduce the transfer operator for the "partition sum" in Eq. ( |6.16| ), we define the 
function 



W(Q K , Qi; L/2£) = j D[Q 2 ] • • • D[Q K . X \ exp £ Str(Q m - Q,) 2 \ , (6. 



19) 



ii 



We use here a brief notation for the matrix elements of the Q-matrices. Each of the two upper indices 



runs from 1 to 8 according to ordering of components of a supervector introduced before Eq. (6.9) 
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which has the property 

W(Q K+1 , Qr, (L + l)/2£) = J D[Q K ] W(Q K+1 , Q K ; £/2§W(Qk, Qi, L/2£) 

= j D\Q K ] exp (J^ Str(Q K+1 - Q K ) 2 ^j W{Q K , Q l5 L/2£). (6.20) 

In the continuum limit £ 3> £, when Q« and Qi+i are close to each other, one can reduce the 
integral equation (|6.20|) to a differential one: 



d t W = AqW, 

\imW(Q,Q';t) = S(Q,Q') 



(6.21) 
(6.22) 



where Aq is the Laplace operator on the nonlinear space of Q-matrices. Eq.( |6.21| ) has a 
form of the heat (or diffusion) equation on this space, with W being the corresponding heat 
kernel. 



6.2 Conductance fluctuations 



Expressions Q6.15 ), Q6.16Q contain the effective number of open channels, iV c h, as a parameter. 
As was shown in |]63fl , the existence of this parameter sets an additional length scale L$ = 
£,/N c h. For a very short system (L <C L ), the conductance (g) is controlled by the sample- 
leads coupling and is of order of In this case, the zero-mode (i.e, spatially constant Q) 
approximation can be used. We are interested here in the opposite limit L ^> L , where the 
resistance of the sample is dominated by the bulk of the system. In this case, N c ^ becomes 
an irrelevant parameter and the heat kernel itself already contains all information about (g) 
and (g 2 ). Since we assume that N ch ^> 1, the crossover length scale L is much less than £, 
so that the ratio of the sample length to the localization length L/£ in the bulk-dominated 
regime can be arbitrary. 

Using the generalized Fourier expansion of the heat kernel in the eigenf unctions of the 
Laplace operator Aq, one can find after very lengthy calculations |64|, |65j (there was a subtle 
error in the GSE case in these papers corrected in [p^D the following exact results for (g) 
and (g 2 ) for an arbitrary value of the parameter 



GOE : 

(9 n )(L) 



2^ 



roc 

- / c/Atanh 2 (7rA/2)(A 2 + l)"Vn(l,A,A)exp 
2 Jo 

/oo 
dX l dX 2 l{l 2 - l)A 1 tanh(7rA 1 /2)A 2 tanh(7rA 2 /2 

x p n (l, Ai, A 2 ) ] [ (-1 + al + ia 1 X 1 + ia 2 X 2 )~ 1 exp 



+ 2 



f,(Tl ,CT2=±1 



-^(/ 2 + A 2 + A 2 + l) 



GUE : 

(g n )(L) 

GSE : 

(g n )(L) 



2G2N-1 



dAAtanh(7rA/2)/(A 2 + l 2 )~ 2 p n (l, A) exp 



L 



(I 2 + A 2 ) 



2 5 E 

h,l 2 e 2N- l 
k + h e 4N - 2 



dXX(X 2 + 1) tanh(7rA/2)/ 1 / 2 p n (A, l 2 , h) 
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Figure 6: Average conductance (g) normalized to its Ohm's law value as a function of the 
sample length (measured in units of £ = 2ttuAD). The full, dashed and dot-dashed lines 
correspond to the orthogonal, unitary, and symplectic symmetry, respectively. 



n 



-1 + iaX + a J i + a 2 h) 1 



cxp 



a ,01 ,f2=±l 

The polynomials p n in the above expressions are given by 



-^? + Z 2 2 + A 2 



GOE 



GUE 



GSE 



Pi(Z,Ai,A 2 ) = Z 2 + A 2 + A 2 + 1, 

Pi{l, Ai, A 2 ) = h\i + \ A 2 + 2l A + 3Z 2 (A 2 + \\) + 2l 2 - \j - \j - 2) 
Pl (l,\)=l 2 + \ 2 , 
p 2 (/,A) = ^(/ 2 + A 2 ) 2 , 



Pl (x,h,h) = x + ii + i*- 1, 

p 2 (A, h, l 2 ) = hlf + l* + 2A 4 + 3A 2 (/ 2 + / 2 ) - 2A 2 + l\ + / 2 - 2). 



(6.23) 



(6.24) 



The results (|6.23| ) are presented in Figs. || |^ for all the three symmetry classes. The 
above GUE formula is written for the case of broken time-reversal symmetry, but preserved 
spin-rotation invariance. The result for systems with broken time-reversal symmetry and 
with strong spin interactions (e.g. systems with magnetic impurities, or systems with strong 
spin-orbit interaction in magnetic field; to be denoted as GUE' below) is 



(«f) GUE '(L) = (l/2)«<<f) UUi W2). 



n/ mGUE/ 



(6.25) 
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Figure 7: The same as Fig. ^, but for the conductance variance var(g). 



The results ( |6.23| ) have been confirmed by numerical simulations of the IWZ-model |67j and 
of the directed network model [Q . In particular, we present in Figs. |] and |9] the numerical 
results of |)7| for the symplectic symmetry class. 

Let us discuss the behavior of Eqs. ( |6.23| ) in the limits of short (L <C £) and long (L £) 
wire. Condition L <C £ corresponds to the metallic (weak localization) region, (g) ^> 1. We 
find in this case the following perturbative (in L/^) expansion for (g), (g 2 ) and var(g): 



(g)(L) 
(g 2 )(L) 

and var(g(L)) 



2£ 2 2 L 4 /L . 

— 1 1 — + O 

L 3 45 £ 945 V £ 1 



2A_8f 52 
T J ~ 3L + 45 



136 L 
945 £ 



+ ^ 



8 32 L 
15 ~ 315 f 



(6.26) 



The perturbative results for the symplectic case are related to those for the orthogonal class 
via the symmetry relations p2e| |27|, [228] (g n ) Sp (L) = (-l/2) n (g n )°(-L/2) and have the 
form 



(<7>(£) 

(<? 2 >(£) 
and var(g(L)) 



2f 1 1 L 

— H 1 

L 3 90 £ 



1 



1890 



2£Y 4£ 13 17 L 

— + -— H 1 

L / 3 L 45 945 £ 



2 

15 + 315 £ 



(6.27) 
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Figure 8: Numerical data |J7| for the average conductance in the case of symplectic sym- 
metry. The box size N and the number of channels iV c h are equal to 10 (squares), 20 
(diamonds), 60 (triangles), and 100 (stars). Each data point corresponds to an average over 
100 realizations of disorder. The full line is the theoretical prediction. 
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Figure 9: The same as Fig. Kj, but for the conductance variance. 



Finally, for unitary symmetry we have 

(g)(L) 

and var(g(L)) 





2 L 




45£ + 


k 


) + 45 




4 






•(f) 


15 + 





+ 



(6.28) 



The expressions for (g) start from the Ohm's law term 2£/L = AiruAD / L (we remind that 
g is measured in units of e 2 /h = e 2 /2nh and includes factor 2 due to the spin, while v is 
the density of states per spin projection). The other terms constitute the weak- localization 
corrections. The leading terms in var(g) are well-known values of the universal conductance 
fluctuations in the case of quasi-lD geometry. 

The opposite condition L ^> £ defines the region of strong localization. In this case, the 
asymptotic behavior of (g n ), n = 1, 2, is as follows: 



GOE 
GUE 
GSE 



(g n ) 
(g n ) 
(g n ) 



2 -3/2-n 7r 7/2 (e/L) 3/2 ( 
2 3-n 7r 3/2^ /L) 3/2 e 



-L/4£ 



2 15/2-n 7r 3/2^ /L) 3/2 e -L/8 € 



(6.29) 



Let us recall that £ is defined here as £ = 2tti>AD independently of the symmetry. The 
formulas ( 6.29|) demonstrate therefore well-known dependence of the localization length on 
the symmetry of the ensemble, L; oc oc (3. Let us stress, however, that in the GUE' case 



8N 



(broken time-reversal and spin-rotation symmetries), the results for which can be obtained 
via the relation Q6.25Q the localization length is the same as in GSE ||101| , |229| . This is 
because the transition GSE— >GUE' not only changes (3 from 4 to 2, but also breaks the 
Kramers degeneracy, increasing by factor of 2 the number of coupled channels. These two 
effects compensate each other. More generally, localization length is proportional to /3/s, 
where s is the degeneracy factor. 

DMPK equations. As has been already mentioned, calculation of higher moments (g n ) 
with n > 2 and thus of the whole distribution function V(g) has not been achieved in 
the supersymmetry approach because of technical difficulties (necessity to increase the size 
of the Q-matrix with increasing n). The conductance distribution in the localized regime 
L ^> £ can be approximately calculated from the Dorokhov-Mello-Pereyra-Kumar (DMPK) 
approach. Within this approach, pioneered by Dorokhov [ |23(J| | and developed by Mello, 
Pereyra and Kumar ||231|| , one derives a Fokker-Planck (diffusion) equation for the distribu- 
tion of transmission eigenvalues T n (the eigenvalues of the transmission-matrix product ttf). 
This equation is conveniently written in terms of the distribution function P(X\, A2, • • • ; L), 
where A, = (1 — Tj)/Tj and has the form 



N 



(6.30) 



BP 2 ^ d W1 . . T d T _ ln 
9L 7 ^ BX n BX n 

N N 

j=u n i a 

i=l j=i+l 

where I is the mean free path, N is the number of transverse modes, and 7 = (3N + 2-/3. 
Overview of the results obtained within the DMPK approach can be found in the review 
article of Beenakker ||88|| . It has been shown recently |)6| that in the limit N ^> 1 the DMPK 
approach is equivalent to the supersymmetric cx-model considered above. Let us note that the 
DMPK approach is restricted to calculation of transport properties of the quasi-lD conductor 
(which can be expressed through the transmission eigenvalues Tj), while the supersymmetry 
method allows to study all sorts of quantities which can be expressed through the Green's 
functions (e.g. statistics of levels, eigenfunction amplitudes, local density of states etc, - see 
other sections of this article). 

In the localized regime L Nl the distribution function P(xi, . . . , x^; L) with Aj = 



sinh Xi (so that Tj = 1/ cosh Xj) takes the form [2"3"0", |232|, |233 



P(x 1 , ...,Xi) 



2nL 



N/2 N 



H exp 



71=1 



2L\ Xn e. 



, 1< ii <C s 2 < • ■ • < 



(6.31) 

where £ n = 7// (1 + f3n — (3) are the inverse Lyapunov exponents. In the limit L — > 00 all x n 
tend to non-random values L/^ n , which is a manifestation of the Oseledec theorem (matrix 



generalization of the "law of large numbers") |[234j , [235| , p36|| . For L/Nl S> 1 the conductance 



N 



g = G/(e 2 /h) = 2j2T n 



(6.32) 



71=1 



is dominated by x±, e.g. g ~ 2/ cosh x\ ~ 8e J , which implies the Gaussian distribution 
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Figure 10: The average and the variance of the conductance logarithm as a function of 
the sample length (for the symplectic symmetry class). The box size N and the number of 
channels are equal to 10. Each data point corresponds to an average over 100 realizations 
of disorder. The dashed line corresponds to the formula ( |6.33| ). 
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(6.33) 



with the average (lng) = — and the variance var(<?) = -j = —2 (In g). We have already 
encountered the same log-normal distribution in Sec. |3.2.4| , when we calculated the distri- 
bution of the product of the wave function intensities in two points located close to the 
opposite edges of the sample. The result Q6.33j ) is fully confirmed by numerical simulations 
as illustrated in Fig. |I0"| . 

The log-normal form ( |6.33| ) of the conductance distribution holds for g < 1 only (since 
one transmission eigenvalue cannot produce a conductance larger than unity); for g > 1 the 
distribution P(g) decays very fast. The moments (g n ) with n > 1 are thus determined by 
the probability to have g ~ 1, which yields (with exponential accuracy) 



(g n ) ~ exp 



2 7 / 



(6.34) 



in full agreement with the asymptotics of the first two moments, Eqs. ( |6.29|) , found from the 
supersymmetric approach. [Eqs. ( |6.29| ) contain also preexponential factors, for derivation 
of which the approximation ( |6.31|) , ( |6.33| ) of the solution of the DMPK equations is not 
sufficient.! 



Finally, we note that in the strictly ID case Abrikosov ||147|| derived the exact result for 
the conductance distribution function P(g) for arbitrary value of L/l. In the limit L/l ^> 1 
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it approaches the log- normal distribution (also found by Melnikov [ |145| , |146|| ), which can 
be obtained from Eq. (|6.33|) by setting A^ = 1 (i.e. 7 = 2). In fact, essentially equivalent 
results were obtained much earlier in the context of the classical wave propagation, see 



Refs. 237, 238, 239 



7 Statistics of wave intensity in optics 

As was already mentioned in the introduction, our statistical considerations are applicable 
not only to properties of eigenfunctions and energy levels of quantum particles (electrons), 
but also to those of intensities and eigenfrequencies of classical waves. This is related to 
similarity of the stationary Schrodinger equation and the wave equation. Propagation of the 
classical field i/) u with frequency u in an inhomogeneous medium is described by the wave 
equation 

[V 2 + fc 2 (l+Mr))R(r) = 0, (7.1) 

supplemented by appropriate sources. Here k = u/c, with c the speed of propagation in the 
average medium, and /z(r) describes fluctuations of the refraction index. The field ipu can 
describe a component of the electromagnetic or acoustic wave. The impurity diagrammatic 
technique |[240| , [24 1| , |242| , [243|| and the cx-model approach |[244j , |245|| can be developed in full 



analogy with the case of the Schrodinger equation in random potential. 

Let us consider an open system with a permanently radiating source (for example, a 
point-like source would correspond to addition of the term oc 5(r — r ) in the r.h.s. of 
Eq. (|7.1| ). The problem of fluctuations of the wave intensity V ; S( r ) i n such a situation has 



a very long history. Almost a century ago Rayleigh proposed a distribution which bears his 
name: 

P (I) = exp(-J) , (7.2) 

where / is the intensity normalized to its average value, I = 1/(1). A simple statistical 
argument leading to Eq. ( |7.2| ) is based on representing ^(r) as a sum of many random 
contributions (plane waves with random amplitudes and phases). This is essentially the 
same argument that was used by Berry to describe fluctuations of eigenfunctions ipi in 
chaotic billiards and which leads to the RMT statistics of |^( r )| ( see Section |3|). Let us 
mention, however, a difference between the two cases (emphasized by Pnini and Shapiro 
246|| ). In the case of an open system, i[) u is a sum of traveling waves, while for the closed 



system ipi is represented as a sum of standing waves. As a result, the Rayleigh statistics has 
the same form as the statistics of eigenfunction amplitudes in a closed system with broken 
time-reversal invariance (unitary class), where the eigenfunctions are complex. 

Diagrammatic derivation of Eq. ( |7.2| ) is very simple ||247|| (see below); for the case of a 



smooth randomness an essentially equivalent derivation ||248|| using path-integral arguments 
was given. However, similarly to distributions of other quantities studied above (eigenfunc- 
tion amplitude, local DOS, relaxation time etc.), distribution of optical intensities show 
deviations from the Rayleigh law, which will be studied below. More specifically, we will 
consider fluctuations of the intensity I(r, ro) at a point r induced by a point-like source at 
ro, with the both points r and ro located in the bulk of the sample. We will assume quasi-lD 
geometry of the sample with length L much larger than the transverse dimension W (see 



Fig. 11 



Let us note that there was a considerable activity recently in studying statistics of the 
transmission coefficients T a b of a disordered waveguide ||249| , |250| , [72] , [73], [74]]. In this for- 



mulation of the problem, a source and a detector of the radiation are located outside the 
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Figure 11: Geometry of the problem. Points r = (xo,yo,^o) and r = (x,y,z) are the 
positions of the source and of the observation point (detector) respectively. From [ [TT|j . 



sample. The source produces a plane wave injected in an incoming channel a, and the in- 
tensity in an outgoing channel b is measured. The transmission coefficients are related to 
the transmission matrix t (already mentioned in Section |6|) as T ao = \t ao \ 2 . One can also 
define the transmittance T a = T ao summed over the outgoing channels. Finally, the total 
transmittance T = ^2 ab T ab is the optical analogue of the conductance g (fluctuations of 
which were studied in the Sec. ||). Combining the diagrammatic approach with the results 
on distribution of transmission eigenvalues in the metallic regime (following from the DMPK 



equations), Nieuwenhuizen and van Rossum |72j calculated the distribution functions P(T ( 
and P(T ab ) for g 3> 1. The results have the following form: 

dx 

2~Kl 

dx 



exp[x,s a - $ con 



»] 



Pis 



—100 

00 dv 



ab 



2m 



exp 



h xv — <P, 



x 



x 



gin 2 



a/1 + x/g + y/xjg^j 



(7.3) 
(7.4) 
(7.5) 



where s a = T a / (T a ) and s ao = T ab /(T ab ). For large g 1 (i.e. in the "metallic" regime), the 
distribution P(s a ), Eq. (|7.3|), has a Gaussian form, 



P(s, 




exp 



- T («- 



(7.6) 



with log-normal tails at small s a < 1 and an exponential asymptotics at large s a ^> 1 
of the form P(s a ) ~ exp(— ps a ). The distribution of the channel-to-channel transmission 
coefficients s a fe is a close relative of the distribution of the point-to-point transmitted intensity 
J(r, r ) which was discussed above and will be studied in detail below. For not too large s ab 
it has the Rayleigh form, P(s ab ) ~ e~ Sab , with the leading perturbative correction given by 



P(s ab ) 



-Sab 



l + —(s 2 ab -4s ab + 2) 

3g 



Sab < y/g 



(7.7) 



(this form of the correction to the Rayleigh law for the intensity distribution was found in a 
number of papers, see Refs. ||251| , [252| , |248| , |249|| ). In the intermediate region ^fg < s a b < 9 
Eq.O yields 

P(s ab ) ~ exp -s ab ^1 + ^s a6 + . . A (7.8) 

(correction of this type was also obtained earlier, Refs. ||252| , |253|| ). For large s a b > g the 
distribution acquires a stretched-exponential form 



P(s ab ) ~ exp (-2y/gs ab ) 



(7.9) 
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We will return to Eq. ( \l.9\ ) below, when comparing the results for the distribution of /(i"o, r) 
with P(s a b). In the localized regime both distributions P(s a ) and P(s a b) are determined 
by the single (largest) transmission eigenvalue and acquire the same log-normal form as the 
conductance distribution, Eq. (|6.33| ). We note also that van Langen, Brouwer and Beenakker 
74| were able to calculate P(s a ) and P(s a b) in the case of the unitary symmetry, (3 = 2, in the 



whole range of the parameter L/£ (from weak to strong localization). The distributions P(s a ) 
and P(s a b) were studied experimentally by Garcia and Genack ||254| , |255|| and by Stoychev 



and Genack ||256|| ; their findings are in good agreement with the theoretical results. 

Now we return to the statistics of the intensity J(r, r D ). The field at the point r is 
given by the (retarded) Green's function G B (r,r Q ) and the radiation intensity is J(r, r Q ) = 
\Gr(t, r Q )| 2 . The average intensity (/(r, r Q )) is given diagrammatically by a diffuson T(r ly r 2 ), 
attached to two external vertices. The vertices are short-range objects and can be approx- 
imated by a 5-function times (if Air), so that (J(r,r )) = (£/A7t) 2 T(t, r Q ). For the quasi 
one-dimensional geometry, the expression for the diffuson reads 

where I is the elastic mean free path, A is the cross-section of the tube, z-axis is directed 
along the sample, z< = min(;z, z Q ) and z > = max(2, z Q ). We assume that \z — z$\ ^> W. 

The intensity distribution P(I) is obtained, in the diagrammatic approach, by calculating 
the moments (I n ) of the intensity. In the leading approximation ||247| , one should draw n 
retarded and n advanced Green's functions and insert ladders between pairs {G_r, Ga} in all 
possible ways. This leads to (I n ) = n\(I) n and, thus, to Eq. (|7.2j ). 

Corrections to the Rayleigh result come from diagrams with intersecting ladders, which 
describe interaction between diffusons. The leading correction is due to pairwise interactions. 
The diagram in Fig. ^ represents a pair of "colliding" diffusons. The algebraic expression 
for this diagram is 

c(r ' rJ = 2 (s)7(n A 

x T(r,ri)T(r, r 2 )T(r 3 , r Q )T(r 4 , r Q ) 

+ 2 (Vi ■ V 2 ) + 2(V 3 ■ V 4 )] f[ 6 (P " r I , (7-11) 

where k is the wave number and V t acts on r t . The factor (i/An) 4 comes from the 4 external 
vertices of the diagram, the T's represent the two incoming and two outgoing diffusons and 
the expression in the curly brackets corresponds to the internal (interaction) vertex ( "Hikami 
box") [[25 7| , 258] . Finally, the factor 2 accounts for the two possibilities of inserting a pair 



of ladders between the outgoing Green's functions. Integrating by parts and employing the 
quasi one- dimensional geometry of the problem, we obtain (for z Q < z): 

C(z,z )~2(I(z,z )) 2 (l + ^j , (7.12) 

where (I(z,z )) = j^^xHT' ls ^ e avera g e intensity, 

L 3 

7 = 29 L\3z + z ) - 2Lz(z + z ) + 2zl(z - z ) > 1 ' (? ' 13) 
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Figure 12: Diagram for a pair of interacting diffusons. The external vertices contribute the 
factor (£/4n) A . The shaded region denotes the internal interaction vertex, see Eq. ( |7.11| ). 
From fh]. 
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and g = kl£A/3iiL ^> 1 is the dimensionless conductance of the tube. For simplicity, we will 
assume that the source and the detector are located relatively close to each other, so that 
\z — z a \ < L, in which case Eq. ( |7.13|) reduces to 7 = gL 2 /2z(L — z). (All the results are 
found to be qualitatively the same in the generic situation z ~ z — z ~ L — z~ L.) 

In order to calculate (I n ) one has to compute a combinatorial factor which counts 
the number iVj of diagrams with % pairs of interacting diffusons. This number is iVj = 

t (n\/i\)(n/2) 2 \ For not too large n <C 7 1 / 2 it is sufficient to keep the 



[ni 



] 2 /[2 2l i!(n-2z)!] 
1 contribution, yielding the leading perturbative correction [71 



(I n ) = n\ 

or in terms of the distribution function, 



1 + —n(n - 1) 

37 



(7.14) 



P(I) 



1 + — (J 2 -4/ + 2) 

37 



K 7 



1/2 



(7.15) 



For larger n (or, equivalently, /) we have to sum up the series over i, yielding 



in 
(/>» 



[n/2] 



2n2 
37" 



n\ exp(2n 2 /37) 



(7.16) 



Although i cannot exceed n/2, the sum in Eq. ( |7.16| ) can be extended to 00, if the value of n is 
restricted by the condition n C7. Eq. ( |7.16| ) represents the leading exponential correction 
to the Rayleigh distribution. Let us discuss now effect of higher order "interactions" of 
diffusons. Diagrams with 3 intersecting diffusons will contribute a correction of n 3 /^ 2 in the 
exponent of Eq. (|7.16| ), which is small compared to the leading correction in the whole region 

Likewise, diagrams with 4 intersecting 



n < 7, but becomes larger than unity for n > 7 2 / 3 



diffusons produce a n 77 s correction, etc. Restoring the distribution P(I), we find |71 



P(I) ~ exp 



3 7 



r 



(7.17) 



Eqs. ( |7. 14|) , ( [T.15| ) , ( |7.17|) are analogous to those found in the case of transmission coefficient 
statistics, ( fT7T) , (|7.8|) , with the only difference that the parameter g is now replaced by 7/2 
(if we consider the limit zq —>■ 0, z — > L, then 7 — * 2g, so that both results are consistent). 
As has been already mentioned, deviations of this form from the Rayleigh distribution of 
intensities were found earlier by various authors ||251| , |252| , [248| , [253j , [249|| ; a value of the 
parameter governing the strength of deviations (here I/7) depends, however, on the geometry 
of the problem. It should be realized that Eq. ( [7.17]) is applicable only for / 7 ~ g and, 
thus, does not determine the far asymptotics of P{I)- The latter is unaccessible by the 
perturbative diagram technique and is handled below by the supersymmetry method. 

For technical simplicity, we will assume now that the time reversal symmetry is broken 
by some magnetooptical effects (unitary ensemble). The moments of the intensity at point 
r due to the source at r Q are given by 



in 
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16tt 5 
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(7. It 



95 



where S[Q] is the zero-frequency a-model action, 



4 

which reduces in the quasi- ID geometry to 



S[Q] = -T^R I ^rStr(Vg) 2 , (7.19) 



S[Q] = -(gL/8) J dz Str (dQ/dz) 2 . 



Assuming again that the two points r and r Q are sufficiently close to each other, \z — z \ <C L 
and taking into account slow variation of the Q-field along the sample, we can replace the 
product Qi2,bb(z)Q 2 i,bb(zo) by Qi2,bb(z)Q 2 i,bb(z). We get then the following result for the 
distribution of the dimensionless intensity y = (167t 2 /A;q)/: 

P(y) = J dQ6(y + Q 12 ,bbQ2i,bb)Y(Q) , (7.20) 

where Y(Q) is a function of a single supermatrix Q defined by Eq. Q3.8|) . Using the fact that 
the function Y(Q) depends only on the eigenvalues 1 < Ai < oo, — 1 < A2 < 1, we find 



d 


d 2 
^ dy 2 


dy 


Ai 




Ai 





p (y) = \ — + y—) I dX ^ 

K 2 )S(y+l-Xl) . (7.21) 



The function Y(Xi, A2) can be evaluated at g ^> 1 via the saddle-point method (see |34]] and 



Section |4. 3| ) with the result 

K(A 1 ,A 2 )~exp{-|[^ + ^]} . (7.22) 

where Ai = cosh^i, A 2 = cos# 2 (0 < Q\ < 00, < 9 2 < n). In fact, the dependence of Y 
on 82 is not important, within the exponential accuracy, because it simply gives a prefactor 
after the integration in Eq. ( |7.21|) . Therefore, the distribution function P(y) is given by 

P(y) ~ Y(X 1 = x/T+^j, A 2 = 1) ~ exp (- 7 0?/2) , (7.23) 

where 9\ = ln(yT + y + y/y). Finally, after normalizing y to its average value (y) = 2/7, we 
obtain Q: 

,. 1 / , ,,<:.. '.,r/. 1 1 (7.24) 



P(I) ~ exp <J — 



In" J 1 + 27/7 + J 2I/7 



For / < 7, Eq. ( |7.24| ) reproduces the perturbative expansion ( |7.17| ), while for 7 ^> 7 it 
implies the log-normal asymptotic behavior of the distribution P(I): 

lnP(J) ~ -(7/8)ln 2 (87/ 7 ) . (7.25) 

The log-normal "tail" ( [7.25| ) should be contrasted with the stretched-exponential asymp- 
totic behavior of the distribution of transmission coefficients, Eq. (|7.9| ). As was found in 
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|7T| , these two results match each other in the following way. When the points z and zq ap- 
proach the sample edges, zo = L — z <C L, an intermediate regime of a stretched-exponential 
behavior emerges: 



lnP(J) 



Sz 



16(f) 2 | 



9 « 1 « 9 (t 



/ > # — 



(7.26) 



Thus, when the source and the detector move toward the sample edges, the range of validity 
of the stretched-exponential behavior becomes broader, while the log-normal "tail" gets 
pushed further away. In contrast, when the source and the detector are located deep in 
the bulk, zq ~ L — z ~ L, the stretched-exponential regime disappears, and the Rayleigh 
distribution crosses over directly to the log- normal one at I ~ g. 

Let us now describe the physical mechanisms standing behind these different forms of 
P(I) jnj. The Green's function G R (r Q ,r) can be expanded in eigenfunctions of a non- 
Hermitean (due to open boundaries) "Hamiltonian" as 



G*(r , r) = J2^(r )Mr)(k 2 -E { + z 7i 



Since the level widths 7$ are typically of order of the Thouless energy E c ~ D/L 2 , there is 
typically ~ g levels contributing appreciably to the sum. In view of the random phases of 
the wave functions, this leads to a Gaussian distribution of G r (tq,t) with zero mean, and 
thus to the Rayleigh distribution of 7(ro,r) = |G (ro, r)| 2 , with the moments (I n ) = n\. 
The stretched-exponential behavior results from such disorder realizations where one of the 
states ipi has large amplitudes in the both points r and r. Considering both ipi(r ) and ipi(r) 
as independent random variables with Gaussian distribution and taking into account that 
only one (out of g) term contributes in this case to the sum for G R , we find (I n ) ~ n\n\/g n , 
corresponding to the above stretched-exponential form of P(P). Finally, the log-normal 
asymptotic behavior corresponds to those disorder realizations where G R is dominated by 
an anomalously localized state, which has an atypically small width 7$ (the same mechanism 
determines the log-normal asymptotics of the distribution of local density of states, see 
Section |4.3|). 



8 Statistics of energy levels and eigenfunctions in a 
ballistic system with surface scattering 

In the preceding part of this article we considered statistical properties of spectra of dis- 
ordered diffusive systems. Using the supersymmetric a-model approach, we were able to 
demonstrate the relevance of the random matrix theory (RMT) and to calculate deviations 
from its predictions both for the level and eigenfunction statistics. Generalization of these 
results to the case of a chaotic ballistic system (i.e. quantum billiard) has become a topic of 
great research interest. For ballistic disordered systems the cx-model has been proposed ||75|| , 
with the Liouville operator replacing the diffusion operator in the action. It has also been 
conjectured that the same cx-model in the limit of vanishing disorder describes statistical 
properties of spectra of individual classically chaotic system. This conjecture was further 
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developed in ||259| , 75 , f7Tf where the a-model was obtained by means of energy averaging, and 
the Liouville operator was replaced by its regularization — the Perron-Frobenius operator. 
However, straightforward application of the results of Refs. fTO, 111], B5|, E7J to the case of 



an individual chaotic system is complicated by the fact that the eigenvalues of the Perron- 
Frobenius operator are unknown, while its eigenfunctions are extremely singular. For this 
reason the cr-model approach has so far failed to provide explicit results for any particular 
ballistic system. 

In this section, we consider a ballistic system with surface disorder leading to diffusive 
scattering of a particle in each collision with the boundary. This models behavior of a 
quantum particle in a box with a rough boundary which is irregular on the scale of the wave 
length. Since the particle loses memory of its direction of motion after a single collision, this 
model describes a limit of an "extremely chaotic" ballistic system, with typical relaxation 
time being of order of the flight time. (This should be contrasted with the case of a relatively 
slight distortion of an integrable billiard | |196| , |26U| , |261[ |.) One might naively think that all 
results for such a model could be obtained by setting I pa L in a system with bulk disorder. 
In fact, the level statistics in a system with bulk disorder and arbitrary relation between 
mean free path I and system size L was studied in ||168| , |169| , [262|| . However, the results 
presented below show that systems with bulk and surface disorder are not equivalent. 

To simplify the calculations, we will assume a circular geometry of the billiard. A similar 
problem was studied numerically in Ref. [|263| , p64j] for a square geometry. We consider 
only the case of unitary symmetry (broken time- reversal invariance); generalization to the 
orthogonal case is straightforward. We follow Ref. [78 1 in the presentation below. The level 
statistics for the same problem was independently studied in Ref. [7U|. Very recently, the 
same approach was used ||265|| to calculate the persistent current in a ring with diffusive 
scattering. 

Our starting point is the sigma-model for ballistic disordered systems 
effective action for this model has the form 



75, 31 . The 
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Here a 8 x 8 supermatrix field g is defined on the energy shell of the phase space, i.e. g depends 
on the coordinate r and direction of the momentum n. The momentum dependence of the 
field g distinguishes the ballistic a-model from the diffusive case where the supermatrix field 
Q depends on r only. The angular brackets denote averaging over n: (0(n)) = J dnO(n) 
with the normalization J dn = 1. Like in the case of the diffusive a-model, the matrix 
g is constrained by the condition g 2 (r,n) = 1, and can be represented as g = UAU^ 1 , 
with A = diag(l, 1, 1, 1, —1, —1, —1, —1). Since we are interested in the clean limit with no 
disorder in the bulk, the second term in the action (|8.1|) containing the elastic mean free 
time r is zero everywhere except at the boundary where it modifies the boundary condition 
(see below). 

As was explained in Section ^|, the statistical properties of energy levels are governed by 
the structure of the action in the vicinity of the homogeneous configuration of the g-field, 
g(r, n) = A. Writing U = 1 — W/2 + . . ., we find the action in the leading order in W, 



F [W] = -™ J drdnSti 



W 21 [K - iuj) W 12 



\.2) 



where the indices 1,2 refer to the "advanced-retarded" decomposition of W, and K is the 
Liouville operator, K = tipnV. This "linearized" action has the same form as that of a 
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diffusive system, with the diffusion operator being replaced by the Liouville operator. This 
enables us to use the results derived for the diffusive case by substituting the eigenvalues 
and eigenfunctions of the operator K for those of the diffusion operator. 

The operator K should be supplemented by a boundary condition, which depends on 
the form of the surface roughness. As a model approximation we consider purely diffuse 
scattering |[266| , 267]] for which the distribution function ip(r, n) of the outgoing particles is 



constant and is fixed by flux conservation^: 

V?(r,n)=7r/ (Nn>(r,n')4' (Nn) < 0. (8.3) 

i(Nn')>0 

Here the point r lies at the surface, and N is an outward normal to the surface. This 
boundary condition should be satisfied by the eigenfunctions of K. 

The eigenvalues A of the operator K corresponding to angular momentum I obey the 
equation 

1 r 

J { (£) = -1 + - / d6 sin exp [2iW + 2£ sin 6} = 0, (8.4) 

2 Jo 

where £ = RX/vf, and R is the radius of the circle. For each value of Z = 0, ±1, ±2, . . . 
Eq. ( |3.4|) has a set of solutions with = = Q_ k , which can be labeled with 

k = 0, ±1, ±2, . . . (even I) or k = ±1/2, ±3/2, . . . (odd I). For I = k = we have £ o = 0, 
corresponding to the zero mode ip(r, n) = const. All other eigenvalues have positive real 
part Re^fc > and govern the relaxation of the corresponding classical system to the 
homogeneous distribution in the phase space. 

The asymptotic form of the solutions of Eq. ( |8.4| ) for large \k\ and/or \l\ can be obtained 
by using the saddle-point method, 

f 0.66/ ±0.141n/± 0.55mA;, < k < / . . 

Kkl ~ \ (lnfc)/4 + 7ri(fc + l/8), 0</<£; ' [ } 

Note that for k = all eigenvalues are real, while for high values of k they lie close to the 
imaginary axis and do not depend on I (see Fig. [T3| ). 

Level statistics, low frequencies. As was explained in Section |2.2j [see Eq. (|2.27Q ], in 
the range of relatively low frequencies (which for our problem means uj Vp/R, see below) 
the level correlation function R(s = uj/A) has the form 



,2 , 



[TTS) 2 



\irv F J 



R(s) - 1 = 5{s) - SI . n ™ + A ( — ) sin 2 vrs. (8.6) 



The first two terms correspond to the zero-mode approximation and are given by RMT, 
while the last one is the non-universal correction to the RMT results. The information 
about the operator K enters through the dimensionless constant A = ' iki'i where the 
prime indicates that the eigenvalue £oo = is excluded. The value of A, as well as the 
high-frequency behavior of R(s) (see below), can be extracted from the spectral function [H] 

l k 



12 Exact form of the boundary condition depends on the underlying microscopic model. In particular, the 
diffuse scattering can be modelled by surrounding the cavity by a disordered layer with a bulk mean free 



path I and a thickness d 3> I- The corresponding boundary condition |26§|, 269] differs from Eq. (3.3) by 



a parameterless function of order unity. For a re view of the boundary conditions corresponding to various 



microscopic realizations of the rough surface see [ 270 
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Figure 13: First 11 x 11 (0 < k, I < 11) eigenvalues of the Liouville operator K in units of 
Vp/R, as given by Eq. (B^). From Ref. 



According to the Cauchy theorem, Si can be represented as an integral in the complex plane, 



\v F J 2-ni J c (z - iuR/v F ) 2 J t ( z ) 

where the contour C encloses all zeroes of the function Ji(z). Evaluating the residue at 
z = iuR/vF, we find 

S l {u) = -{R/v F f ^ hxJ^z). (8.8) 

dZ z=iu)R/v F 

Considering the limit u — > and subtracting the contribution of A o = 0, we get 

A = -19/27 - 1757r 2 /1152 + 64/(9tt 2 ) w -1.48. (8.9) 

In contrast to the diffusive case, this constant is negative: the level repulsion is enhanced 
with respect to result for RMT. Eq.( fj.6p is valid as long as the correction is small compared 
to the RMT result, i.e. provided u is below the inverse time of flight, vp/R. 

Level statistics, high frequencies. In the range u ^> A the connected part of the level 
correlation function R(u) — 1 can be decomposed into the smooth Altshuler-Shklovskii (AS) 
part R&s(uj) = (A 2 /27r 2 )Re S(u) and the part R osc which oscillates on the scale of the 
level spacing (see Sec. fj.2[ ). Evaluating the asymptotic behavior of S[(u>) from Eq.(|8.8|), we 
find in the high-frequency regime when uo ^> v f/R: 

n , , /Ai?\ 2 / V F \V2 / ujR 7r\ . 
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The oscillating part of the level correlation function R osc (s) for frequencies to ^> A is 
given by Eq. (|2.34 ) with the spectral determinant 



D(s) = s- 2 11 (l-iaA/A^-^l + isA/Aw)- 1 . 

Since A~ 2 d 2 In D(s)/ds 2 = —2ReS(u), we can restore D(s) from Eqs.(|8.7|), ( pT8|) up to a 
factor of the form exp(ci + c 2 s), with cj and C2 being arbitrary constants. These constants 
are fixed by the requirement that Eq.( |2.34| ) in the range A C w C should reproduce 

the low-frequency behavior (|3.6|). As a result, we obtain 



D(s) = f-V - TT i • (8.H) 

Here iV = (vp/RA) 2 = (ppR/2) 2 is the number of electrons below the Fermi level. For high 
frequencies to 3> vpjR this yields the following expression for the oscillating part of the level 
correlation function: 

*- (u > = is (v) cos W (8 - 12) 

It is remarkable that the amplitude of the oscillating part does not depend on frequency. 
This is in contrast to the diffusive case, where in the AS regime (u above the Thouless 
energy) the oscillating part R osc (u) is exponentially small, see Eq. (|2.39|) . 



The level number variance. The smooth part of the level correlation function can be 
best illustrated by plotting the variance of the number of levels in an energy interval of width 
E = sA, 

E 2 (s) = J (s-\s\)R(s)ds, (8.13) 

A direct calculation gives for s N 1 ^ 2 

tt 2 S 2 (s) = 1 + 7 + ln(27rs) + „4s 2 /(2iV) (8.14) 

and for s > N 1/2 

16ArV2 



7r 2 S 2 (s) = 1+7 + lm 



7T 



2 



7T 2 f2N l / 2 \ 1 ' 2 [AS 7T 



. ,cos —-^r - - . (8.15) 

Here 7 ~ 0.577 is Euler's constant, and A is defined by Eq. (|3.9|). The first three terms at 



the rhs of Eq. (|8.14|) represent the RMT contribution (curve 1 in Fig. [T^ 



As seen from Fig. [TJ|, the two asymptotics ( |8.14| ) and ( |8.15| ) perfectly match in the 



intermediate regime, s ~ N 1 / 2 . Taken together, they provide a complete description of 
S 2 (s). According to Eq. ( ^.15|) , the level number variance saturates at the value E 2 _ 



7r _2 (l + 7 + ln(16iV 1//2 /7r 2 )), in contrast to the behavior found for diffusive systems or 
ballistic systems with weak bulk disorder | |168| , |262| . The saturation occurs at energies 
s ~ iV 1 / 2 , or in conventional units E ~ VpjR. This saturation of S 2 (s), as well as its 
oscillations on the scale set by short periodic orbits, is expected for a generic chaotic billiard 



80 , 81] . It is also in good agreement with the results for S 2 (s) found numerically for a 



tight-binding model with moderately strong disorder on boundary sites ||263| , |264| 
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Figure 14: Level number variance ^{E) as a function of energy; s = E/A. Curve 1 shows 
the RMT result, while curves 2 and 3 correspond to asymptotic regimes of low (|8.14) and 



high ( |B.15| ) frequencies. The saturation value E 2 °^ is given in the text. From Ref. [[78 



Eigenfunction statistics. Now we turn to correlations of the amplitudes of an eigen- 
function in two different points defined by Eq. ( p.67|) . As was discussed in Sec. |3.3.3| , these 
correlations are governed by the ballistic propagator IIs(ri, r 2 ) (see Eq. ( |3.98| ) and the text 
preceding it). Direct calculation gives: 



n B (ri,r 2 ) 



U 1 (r 1 ,r 2 ) + L T 2 (r 1 ,r 2 ) 



(8.16) 



n^n.ra) = ILfin-Tj-V- 1 J dring^ri-ra) (8.17) 
-V- 1 J dr 2 4 0) (ri - r 2 ) + V~ 2 J dr'^uf^ - r' 2 ); 

where ng (r) = 1/(ttpf\t\), and (r,8) are the polar coordinates. This formula has a clear 
interpretation. The function Ub can be represented as a sum over all paths leading from ri 
to r 2 , with possible surface scattering in between. In particular, Hi corresponds to direct 
trajectories from r x to r 2 with no reflection from the surface, while the contribution n 2 is due 
to the surface scattering. The first term in the numerator 4k 2 — 1 comes from trajectories with 
only one surface reflection, while the second sums up contributions from multiple reflections. 

Let us summarize the main results of this section. We have used the ballistic cr-model 
approach to study statistical properties of levels and eigenf unctions in a billiard with diffusive 
surface scattering, which exemplifies a ballistic system in the regime of strong chaos. It was 
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found that the level repulsion and the spectral rigidity are enhanced compared to RMT. 
In particular, the level number variance saturates at the scale of the inverse time of flight, 
in agreement with Berry's prediction for a generic chaotic system |80, [Sl|. As another 



manifestation of the strong spectral rigidity, the oscillating part of the level correlation 
function does not vanish at large level separation. We calculated also the ballistic analog of 
the diffusion propagator in this model, which governs correlations of eigenfunction amplitudes 
in different spatial points. 

While we focused our attention on the statistics of levels and wave functions in a closed 
ballistic sample, the surface nature of scattering will also modify statistical properties of 
the transport characteristics for an open system. In this connection, we mention the recent 
papers ||271 , |272|| where quantum localization and fluctuations of the transmission coefficients 
and of the conductance were studied for a quasi- ID wire with a rough surface. A number of 
important differences compared to the case of bulk disorder was found. 

Note that the motion in a quasi- ID wire with surface scattering is closely related to the 
PRBM ensemble of Sec. [5731 According to (|8.3| ), the probability density for a particle to 
leave the surface after a scattering event with an angle S C 1 with respect to the surface 
is V(6) ~ 9. Since the distance of the ballistic flight is r ~ O^ 1 for small 6, this yields 
V(r) ~ r~ 3 for r larger than the transverse size of the wire. Thus, we get a power-law 
"tail" at large r of the form (|5.43|) with a = 3/2. According to Sec. [57|, this is precisely 



the marginal value separating the regions of conventional KQ 1 (q, 0) oc q 2 and unconventional 
Kq 1 ^, 0) oc q 2a ~ x behavior of the diffusion propagator K (q,u). It is clear that at a = 3/2 
the propagator acquires a logarithmic correction, iT ( T 1 (g, 0) oc q 2 In q, see Ref. ||271|| for details. 



9 Electron-electron interaction in disordered mesoscopic 
systems 

In the preceding sections we have considered statistical properties of energy levels and eigen- 
functions of a single particle in a disorder potential. However, if an electronic system is 
considered, the Coulomb interaction between the electrons has to be taken into account. 
The influence of the electron-electron interaction on transport properties of disordered sys- 
tems has been intensively studied during the last two decades, in particular in connection 
with such phenomena as weak localization [|273|| and universal conductance fluctuations [|274 



The electron-electron interaction sets the length scale 1$ (phase breaking length) below which 
the electron wave function preserves its coherence. Also, interplay of the interaction and dis- 
order leads to a singular (at Fermi energy or at zero temperature) correction to the density 
of states and to the conductivity [ 273| . 



More recently, another kind of problems has attracted the research interest: to what ex- 
tent does the electron-electron interaction influence the properties of the electron spectrum 
in a disordered dot? This ineterest is largely motivated by a recent progress in nanotech- 
nology which allows to observe experimentally discrete electronic levels in semiconductor 
quantum dots [53, |275| and in small metallic grains |276| , |277| . In fact, there are two types of 



the quantum dot spectra studied experimentally via measuring their I-V characteristics: (i) 
excitation spectrum, when excited levels are probed in a dot with given number of electrons 
by increasing the source-drain voltage, and (ii) addition spectrum, when electrons are added 
one by one by changing the gate voltage. 

As concerns the excitation spectrum, it was found in ||278| , |99|j that the quasiparticle 



levels with energies below the Thouless energy E c (counted from the Fermi energy) have a 
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width smaller than the mean level spacing and thus form a well-defined discrete spectrum. 
To show this, let us calculate the r.m.s. value of the matrix element Vijki of the screened 
Coulomb interaction which describes a decay of the quasiparticle state \i) to a three-particle 
(more precisely, two particles + one hole) state \jkl). Using Eqs. (|3.78|) , ( |3.88|) for a diffusive 



dot, one finds |99| ((V^l 2 ) ~ (A/g) 2 . On the other hand, the density of states of the three- 
particle states is u 3 (E) ~ E 2 /A 3 , so that the Golden Rule width of the one-particle states 
is T(E) = 2ir(\V\ 2 )v 3 (E) ~ E 2 /Ag 2 . Comparing T(E) with A, we get the above threshold 
E ~ gA = E c . The number of the discrete excited levels is thus of order of the intrinsic 
dimensionless conductance g of the dot (not to be confused with the tunneling conductance in 
the Coulomb blockade regime determined by the contacts). These findings are in agreement 
with experiment p3 . 

The fine structure of these excited single particle levels was studied theoretically in 
80| (see also [325, 326 ]) where the Hamiltonian of the many-body interacting 
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problem was considered as a tight-binding model in the Fock space, with matrix elements 
of the Coulomb interaction playing a role of the hopping terms. It was found that only for 
the energy above E* ~ s/gA is the Golden Rule applicable and the levels have the regular 
Lorentzian shape. For smaller energies the quasiparticle excitation consist essentially of a sin- 
gle peak with a small admixture of other many-particle exact eigenstates. This corresponds 
to the Anderson localization in the Fock space. 

Though the properties of the single-particle excitations discussed above are of most in- 
terest, one can also discuss statistical properties of exact many-body levels. For excitation 
energy £> Aa typical state consists of ~ ^JE/ A ^> 1 quasiparticles, each of them having 
an energy ~ \J EA. Comparing the density of those states to which any of such states is 
coupled with the matrix element of the interaction, one estimates the border for the complete 
derealization of many-body states as E ch ~ Ag 2/3 [fM], |7| (see also ||282| ). For energies 
E <C E c h the many-body level statistics is Poissonnian, while for the energies above E c h it 
should acquire the Wigner-Dyson form. Numerical studies of statistical properties of spectra 
of many-body fermionic systems ||283 , [284j , |285| , [286| , |287|| indeed show such a crossover from 
the Poisson to RMT statistics with increasing excitation energy. 

Statistical properties of the addition spectrum have been studied by several experimental 
groups recently |83|. |288| . |36l , [289| , |290fl . We will address the theoretical aspects of this problem 



(following Refs. @, |3§) in Sec. gj]. 

Let us briefly mention another topic which has attracted a great deal of research interest 
recently. This is the problem of localization of two interacting particles raised by Shep- 
elyansky ||291|| (see also an early paper by Dorokhov ||292|| ). It was suggested in |[291|| that 
the effective localization length £2 of two interacting particles in a ID disordered sample 
can be much larger than the localization length £1 of a non-interacting particle. This paper 
stimulated a considerable analytical and numerical activity. In the quasi-lD case, the prob- 
lem can be mapped onto an ensemble of random banded matrices with strongly fluctuating 
diagonal elements, which was studied numerically in [ |293| and analytically (via the super- 
symmetric cr-model approach) in ||294| , |295|| . The results confirmed the original conjecture of 
Shepelyansky, £ 2 oc £f. In the case of a strictly ID system, the situation is somewhat more 
complicated, since the corresponding random matrix ensemble was found |JTJ to be of the 
power-law random banded matrix type (see Sec. |5.3j ) in the critical regime a = 1. Still, the 
localization length enhancement was predicted, £2 oc £f , but with a non-universal exponent 
H (1 < \i < 2) depending on the interaction strength. Generalization of the Shepelyansky' s 
idea to the vicinity of the Anderson transition point in higher-dimensional systems was pro- 
posed by Imry ||296|| . A detailed review of recent activity in this direction can be found in 
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9.1 Coulomb blockade: fluctuations in the addition spectra of 
quantum dots 

The Coulomb blockade effect shows up the most distinctly in the addition spectrum experi- 
ment, where one measures the conductance of an almost isolated quantum dot as a function 
of the gate voltage. Detailed description of the experimental setup can be found e.g. in 
82], [29711 . At low temperature T one observes a sequence of almost equidistant peaks sep- 



arated by "valleys" where the conductance is very small. The positions of the peaks are 
determined by the condition that the energies of the dot with N and N + 1 electrons are 
equal. We will assume below that the following inequality is met: T <C kT <C A, where 
r = + T r is the sum of the tunneling rates to the leads. The peak width is then set by 
the temperature, while the height is given by ||298 | 



e 2 7r r^rv 

feax = ~h2kTT l + T r - (9 ' 1} 

The peak heights show strong fluctuations induced by RMT-like fluctuations of eigenfunction 
amplitudes, as predicted in [13|, Q and observed experimentally in 0, j95fl . In the valleys 
between the peaks, it is energetically costly (of order of the charging energy e 2 /C, where 
C is the dot capacitance) to add an electron to (or to remove from) the dot, and the 



conductance is determined by virtual processes (so-called elastic cotunneling ||299|| ) and is 
strongly suppressed. 



The issue of the statistics of the peak spacings was addressed for the first time in |83 



Basing on numerical data for very small systems, the authors of |83| concluded that the 



r.m.s. deviation of the peak spacing Sn is proportional to the charging energy e 2 /C, with a 
coefficient « 0.15. We will show below (our consideration will closely follow Refs. [Q, |136|| ) 



that the fluctuations are in fact much smaller, of order of the mean level spacing A. Let us 
note that in an analogous problem for classical particles the fluctuation magnitude rms(S'Ar) 
would indeed be proportional to the mean value (Sn) IPOOl |301|| . The physical reason for 



smaller fluctuations in the quantum case is in the delocalized nature of the electronic wave 
functions, which are spread roughly uniformly over the system. These theoretical conclusions 
were confirmed recently by thorough experimental studies, as discussed in the end of the 
section. 

The simplest theoretical model which may be used to study distribution of the peak 
spacings is as follows. One considers a dot as a fixed size diffusive mesoscopic sample and 
assumes that changing a gate voltage by an amount SV g simply reduces to a uniform change 
of the potential inside the dot by a constant ^SVg, with certain numerical coefficient 7 
("lever arm"). Such a model was used for numerical simulations of the addition spectra in 
Refs. [|^, |302j| . We start by considering the statistics of peak spacings within this model |84 



we will later return to the approximations involved and relax some of them. We will neglect 
the spin degree of freedom of electrons first; inclusion of the spin will be also discussed in 
the end of the section. 

The distance between the two consecutive conductance peaks is given by (Fig. [T5]) 



Sn — (£n+2 ~ £n+i) — (£n+i — &n) 

(9.2) 
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Figure 15: Energy levels before and after the electron #(iV + 1) has been added to the dot. 



From Ref. 84 . 



where £n is the ground state of a sample with N electrons. In the second line of Eq.( 
we rewrote Sn in terms of the Hartree-Fock single electron energy levels, with [i\ denoting 
the energy of the state #j in the dot containing % electrons. It is convenient to decompose 
Sjv in the following way 

- Wn+1 ~ Vn ) + l/% ~~ MiV ) 

= E l + E 2 (9.3) 

The quantity E 2 is the distance between the two levels of the same one-particle (Hartree- 
Fock) Hamiltonian (describing a dot with N electrons) and is expected to obey RMT; 
in particular {E 2 ) = A and r.m.s.(£/2) = with a numerical coefficient a of order of unity 
[a = 0.52 (0.42) for the orthogonal (resp. unitary) ensemble]. On the other hand, E\ is a 
shift of the level #(N +-2) due to the change of the Hamiltonian — > Hn+i accompanying 
addition of the electron #(iV + 1) to the system. 

The effective interaction U(r, r') between the electrons #(N + 1) and #(N + 2) can be 
found from the RPA-type equation, 

U(r, r') = U (r - r') - J dr lC lr 2 U (r - ri )P ( ri , r 2 )C/(r 2 , r') (9.4) 

where Uo(r) = e 2 /er, e is the dielectric constant, and 

j Po(ri,r 2 ) = z/[5(r 1 -r 2 )-y- 1 ] (9.5) 
is the polarization operator. Its solution has the form |84| 



U(r, r') = U + 5U(r) + 5U(r') + U K (r, r'). (9.6) 
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Here U = e 2 /C is a constant (charging energy), 5U(r) is the change of the self-consistent 
potential due to addition of one electron (i.e. difference in the self-consistent potential in 
the dots with N and N + 1 electrons), and U K (r) is the screened Coulomb interaction. In 
particular, in the experimentally most relevant 2D case (which we will consider below) and 
assuming a circular form of the dot with radius R, we have 

while U K is given in the Fourier space by U K (q) = 2ite 2 /e(q + k s ) with the inverse screening 
length k s = 2ire 2 v/t. 

According to Eq. (|9.6|) , the term E\ can be decomposed into the following three contri- 
butions 

Ei = e 2 /C + jdr{\^ N+1 (r)\ + \^ N+2 (r)\)5U(r) 
j drdr'|^ +1 (r)||^ +2 (rOI^(|r-r'|) 



+ 



= E f) + sji) + E f) (9. 8) 

The first term in Eq. (|9.8|) (the charging energy) determines the average value (E\) and 
thus the average peak spacing (Sn) (since e 2 /C ^> A for a large dot with N 3> 1). This 
is the only contribution to E\ kept by the so-called constant interaction model, which in 
addition neglects fluctuations of the capacitance C. Consequently, fluctuations of Sn within 
the constant interaction model are determined solely by fluctuations of the single-particle 
level spacing E 2 and thus follow RMT: r.m.s. (S^r) = aA. 

The term E\ in Eq. ( |9.3| ) is however an additional source of fluctuations and is thus 
responsible for the enhancement of fluctuations in comparison with RMT. In principle, all 
the three terms E± , E^\ and E^ in Eq. ( |9.8| ) contribute to this enhancement. Fluctuations 
of the first one, E^ = e 2 /C are due to the fact that the capacitance is slightly different from 
its purely geometric value because of a finite value of the screening length. The corresponding 
correction to C can be expressed in terms of the polarization operator P(r,r'). The latter 
is a fluctuating quantity (because of fluctuations of the eigenfunctions in the Fermi sea) and 
contains a random part P r (r,r') leading to the following expression for the random part of 
the charging energy: 

(e 2 /C) r = 2 J drdr'8U{r)P r (r,r')5U{r'). (9.9) 



Evaluating the fluctuations of the polarization operator, we find [138 



/77(0h 48 2, 

vai{hl ) = —v mg 

P 



^ J dx 1 dT^SU{T 1 )U{T li T^)SU{T^ 



l. M x2 



oc - p ln 9 {-) (9.10) 

,(2) 



Now we consider fluctuations of the last term, E\ , in Eq.(|9.8|). Using Eqs. ( |3.74| ), ( |3.86| ) 
the corre 
found to be 



(2) 

for the correlations of eigenfunction amplitudes in two remote points, the variance of Ef J is 



var(Pf ) ) = J dT 1 dr' 1 dr 2 dr' 2 U K {\Y 1 - r' 1 \)U K {\r 2 - r' 2 |)n(ri, r 2 )n(r / 1 , r' 2 ) 
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4A 2 f 
~ ~pr^F2 J dr 1 dr 2 U 2 (r 1 ,r 2 ) 

« Mi) 2 

Finally, fluctuations of the term E± can be also evaluated with a help of Eqs. ( p . 74j ) , ( |3.86| ), 
yielding 

var(4 1} ) = ^ J dr 1 dr 2 5U(r 1 )n(r 1 ,r 2 )5U(r 2 ) 

1 A 2 . . 

oc -_. (9.12) 



It is seen that for g ^> 1 all the contributions Eqs. (|9.10| )-( |9.12| ) are parametrically 
small compared to the RMT fluctuations (which are ~ A). Fluctuations of the term 
related to the change 5U(r) of the self-consistent potential represent parametrically leading 
contribution to the enhancement of the peak spacing fluctuations with respect to RMT. 

Let us now discuss approximations made in the course of the above derivation: 

i) The dot was supposed to be diffusive in the calculation. For a ballistic dot one should 
replace IT(r,r / ) by its ballistic counterpart. This would mean that the parameter g is 
replaced by ~ iV 1 / 2 ~ L/\p, where N is the number of electrons in the dot and L 
the characteristic linear dimension. The numerical coefficient would depend, however, 
on "how strongly chaotic" is the dot. The role of the eigenf unctions fluctuations and 
correlations ( "scars" ) in enhancement of the peak spacing fluctuations was studied in 
[ p03|l via numerical simulations of a dot with N « 100 electrons. 

ii) It was assumed that changing the gate voltage results in a spatially uniform change of 
the potential in the sample. This has led us to the expression Eq. (|9.7|) for the change 
of the self-consistent potential 5U(r) accompanying the addition of one electron to the 
dot. This result would correspond to a gate located far enough from the sample. In 
a more realistic situation, when the gate is relatively narrow and located close to the 
sample, the potential change SU(r) (as well as the additional electron density) will be 
concentrated on the side of the dot facing the gate. The change of the potential 5U(r) 
corresponds then to a slight deformation of the dot with adding each electron to it. To 
estimate 5U (r) in this case, we can consider a model problem of a point-like charged 
object (modeling the gate) located a distance d from the edge of the dot. Assuming 
that the dot size is larger than d, we can approximate the dot [while calculating 5U (r)] 
by a half plane. In one period of the Coulomb blockade oscillations the gate charge 
changes by e (electron charge). The charge distribution induced in the dot is 

6p(x,y) = ~\ - ; 1 (9.13) 

where the closest to the gate point of the dot is chosen as the coordinate origin and 
the x axis is directed along the dot edge. The corresponding change of the potential 
is 5U{r) = (e 2 z//e) _1 5p(r). Substituting this into the first line of Eq. (|9.12j ), we come 

to the same result 

r.m.s.^f)) oc A/y/fig, 
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in the diffusive limit d > I, and to 

T.m.s.(E[ 1] ) ~ A/ v /f3k F d 

in the case d < I. 

Effect of the quantum dot deformation on the peak spacing fluctuations has been 
recently considered in Ref. |304|| . The authors of [|304j1 characterized the strength of 



the deformation by a phenomenological dimensionless parameter x and assumed that 
it can be large (x > 1), strongly affecting the spacing distribution. However, the above 
estimates indicate that for a typical geometry this parameter is much less than unity, 
x ~ 1/ y/g or x ~ 1/ y/kpd. 

It was assumed that the dot energy and the measured gate voltage are related through 
a constant (or smoothly varying) coefficient 7. This "lever arm" 7 depends, however on 
the dot-gate capacitance, which is also a fluctuating quantity. If the gate size and the 
distance to the gate is of the order of the size of the dot, these fluctuations should be of 
the same order as fluctuations of the dot self-capacitance and thus lead to additional 
fluctuations parametrically small compared to A, see Eq. ( |9 . 1 0[ ) . 



) The calculation was done within the random phase approximation, which assumes 
that the ratio of the interparticle Coulomb energy to the kinetic energy is small, 
r s = \/2e 2 /ehvF C 1 (wf is the Fermi velocity). However, most of the experimen- 
tal realizations of semiconductor quantum dots correspond to r s ~ 1. Since this value 
is still considerably lower than the Wigner crystallization threshold, the calculations 
should be still valid, up to a numerical factor a{r s ) [depending on r s only and such 
that a(r s < 1) = 1]. 

) We considered the model of spinless electrons up to now. Let us briefly discuss the role 
of the spin degree of freedom. Within the constant interaction model, it would lead to 



a bimodal distribution ||302|| of peak spacings 

V(S N ) = \ [s(S N - e 2 /C) + -L PwD ( Sn ~^ /C )] , (9.14) 

where Pwd{s) is the Wigner-Dyson distribution and A denotes the level spacing in the 
absence of spin degeneracy. The value of the coefficient a in the relation r.m.s.(Sjv) = 
aA is then increased (compared to the spinless case) and is equal to 1.24 (1.16) for the 
orthogonal (resp. unitary) ensemble. Taking into account fluctuations of eigenf unctions 

(and thus of E\) however modifies the form of the distribution. The value of the term 

(2) 

E\ representing the interaction between two electrons is larger in the case when 
4>n+2 an d ipN+i correspond to two spin-degenerate states (i.e. have the same spatial 
dependence of the wave function), since 

(J drdr'|^ 2 (r)||^ 2 (r')|C/«(|r-r'|)^ - Q drdr'\tf(r)\\^(r')\U^ ~ r'l)) 

= Jy~2 J dvdv'k d {\v - r'\)U K (\r - r'|) ~ A (9.15) 

for r s ~ 1 (the coefficient depends on r s , see [|4]]). Therefore, filling a state ip^ pushes 
up the level ipii (with respect to other eigenstates) by an amount of order of A. This 
removes a bimodal structure of the distribution of peak spacings and slightly modifies 
the value of the coefficient a. 
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Basing on the above analysis, we can make the following general statement. Imagine that 
we fix r s ~ 1 (i.e. fix the electron density and thus the Fermi wave length) and the system 
geometry, and then start to increase the linear dimension L of the system. Then, while the 
average value of the peak spacing Sn scales as (Sn) ~ e 2 /C oc 1/L, its fluctuations will scale 
differently: r.m.s.(Sjv) ~ A oc 1/L 2 . This conclusion is also corroborated by diagrammatic 
calculations of Ref. |^5| . 

As was mentioned in the beginning of this section, this result should be contrasted with 



that for classical particles [ P00| , |301|| , where the fluctuations are proportional to the mean 
value (Sn)- The parametrically smaller fluctuations in the quantum case are due to the 
delocalized nature of the electronic wave functions, which are spread roughly uniformly over 
the system. 

The above prediction was confirmed by recent experiments [^, |87], |290|| , where a thorough 
study of the peak spacing spacing statistics was carried out. It was found that the low- 
temperature value of r.m.s.(SV), as well the typical temperature scale for its change are 
approximately given by the mean level spacing A (while in units of E c the magnitude of 
fluctuations was as small as 1-4%). 

Several recent papers studied statistical properties of the peak spacing numerically. Stopa 
[ 303 ] used the density functional theory and found the fluctuations of the addition energies 
to be approximately 0.7A, in agreement with the above results. Disappearance of the bi- 
modal character of the peak spacing distribution with increasing r s was observed recently 
in numerical simulations by Berkovits [|305|| via exact diagonalization. However, the system 
size and the number of particles in the exact diagonalization studies are very small, so that 
it is difficult to draw any quantitative conclusions concerning the dots with large number of 
electrons from the results of ||305|| . 



10 Summary and outlook 

In this article we have reviewed the recent progress in the study of statistical properties 
of disordered electronic systems. We have discussed statistics of energy levels and eigen- 
function amplitudes, as well as of several related quantities (local density of states, escape 
time, conductance, etc.) In most of the article the supersymmetric cr-model approach was 
used, as a unique and powerful tool allowing one to calculate various distribution and cor- 
relation functions. Within this approach, we have employed a number of complementary 
methods of treating the cr-model: exact solution (in particular, transfer-matrix method in 
ID), perturbation theory, renormalization group, saddle-point method. The main emphasis 
has been put onto system-specific deviations from the universal predictions of the random 
matrix theory. The results presented constitute a detailed and, in many respects, complete 
description of fluctuations of spectra and wave functions of disordered systems. 

Still, there is a number of directions in this field in which a more complete understanding 
is needed, so that the corresponding research remains active at present. Let us point out 
some of them. 

• Statistical properties of energy levels and wave functions of ballistic systems whose 
classical counterparts are chaotic require more theoretical studies. It remains to be 
understood what are the conditions of applicability of the ballistic cr-model derived in 



76, [77]. Apparently, a certain amount of disorder present in a system (or, in other 



words, some ensemble averaging) is needed to justify the derivation in these papers, but 



the precise conditions have not been quantified yet (see Ref. ||306|1 where this problem is 
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discussed in the context of correlations of quantum maps). Also, a discrepancy between 
the results of the cx-model and of the semiclassical approach in treating repetitions of 



periodic orbits [ |307| , 308 ] is to be resolved. It remains to be seen whether the ballistic 
cx-model approach can be developed to predict the weak-localization corrections (recent 
attempts [ |309| , |310|| did not lead to any definite conclusions), the asymptotic "tails" of 



the distributions (like those discussed in Section f| for diffusive systems) etc. 

Interplay of disorder and electron-electron interaction in mesoscopic systems continues 
attracting considerable research interest. In Sec. |]of this article we discussed only the 
issue of fluctuations in an (almost) closed system. In the experiment, coupling of the 
dot to the outside world is controlled by the gate(s) and may be varied. There arise 
a rich variety of regimes depending on the coupling strength, magnetic field, and the 
strength of the Coulomb interaction (parameter r s ). Fluctuations in open quantum 
dots have been studied recently, both theoretically ||311| , |312| , |313|| and experimentally 



1 |315| , |289| , |29U| |, but the issue is not fully understood yet. In Ref. | ]316| | it was found 
that application of a strong magnetic field affects dramatically the addition spectrum 
of a quantum dot, leading to strong fluctuations and to bunching of the Coulomb 
blockade peaks. A consistent explanation of these experimental results is still missing. 

Very recently, several groups [[317| , |318| , |319| , |320| , j321|] simultaneously performed self- 



consistent Hartree-Fock calculations of addition spectra of quantum dots. Though the 
method allows one to study considerably larger systems as compared to the exact diag- 
onalization method, the obtained results were not sufficient to extract unambiguously 
a scaling dependence of the spacing distribution function on the parameters of the 
problem. A drawback of the Hartree-Fock method with the bare Coulomb potential 
is that the exchange interaction does not get screened, whereas the screening was cru- 
cially important for the theoretical consideration in Section |9.1| . More work in this 
direction may be expected in the nearest future. 

The physics of the Anderson metal-insulator transition, including fluctuations at the 
critical point, remains an actively studied field. Parameters of modern computers 
allow one to evaluate numerically critical indices and various distribution functions 
with a high accuracy. It would be very interesting to study numerically the critical 
point in higher dimensions (see, in this respect, recent paper [ [L78|| ) as well as in an 



effectively infinite-dimensional tree-like sparse random matrix model and in a long- 
range ID model (power- law random banded matrix ensemble considered in Sec. |5.3|) . 
It remains to be seen what is the status of the conjecture of Ref. EDI relating the 



spectral compressibility to the eigenfunction multifractality. Also, the form of the 
conductance distribution at the critical point remains an open problem. Ref. [p8| 
predicted a log- normal form of the distribution function at g ^> (g) ~ 1, similarly to 
the distribution functions of local density of states and of relaxation times. However, 
in contrast to the latter two quantities, an anomalously large conductance cannot be 
explained by a single anomalously localized state (since a single state cannot produce 
the conductance larger than 1). It seems much more probable that the conductance 
distribution falls off much more fast, in a Gaussian (or similar) fashion at large g, 
but no corresponding calculation has been available so far. Technically, the problem 
is that calculating higher moments of the conductance requires increasing of the size 
of the Q-matrix in the cx-model, which complicates tremendously a non-perturbative 
treatment of the distribution function. For recent numerical studies of the problem see 
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Refs. |27], |328|, |329 
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A Abbreviations 



AKL 


Altshuler, Kravtsov, and Lerner, Ref. 28 


ALS 


anomalously localized state 


DMPK 


Dorokhov-Mello-Pereira-Kumar (equations) 


DOS 


density of states 


GOE 


Gaussian orthogonal ensemble 


GSE 


Gaussian symplectic ensemble 


GUE 


Gaussian unitary ensemble 


IPR 


inverse participation ratio 


LDOS 


local density of states 


LN 


logarithmically- normal 


PRBM 


power-law random banded matrix 


RBM 


random banded matrix 


RG 


renormalization group 


RPA 


random phase approximation 


RMT 


random matrix theory 


WD 


Wigner-Dyson (level statistics) 


ID, 2D, 3D 


one-dimensional, two-dimensional, three-dimensional 
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